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This  technical  note  is  the  final  report  covering  the  period  1  Oct  8U  to  30 
Sep  81  in  response  to  NASA  A83011B  sponsored  by  Nasa  Ames  Research  Center 
and  administered  by  Mr.  Tom  Gregory.  Capt.  G.  Sisson  was  the  principal 
^investigator  and  was  aided  by  staff  and  cadets  at  the  Air  Force  Academy. 
This  technical  note  consists  of  two  separate  papers.  The  first  covers  the 
compressible  calibration  of  the  seven-hole  probe  and  is  a  reprint  of  an 
article  in  the  Aeronautics  Digest  -  Fall/Winter  1980,  USAFA-TR-81-4 .  The 
second  paper  describes  the  associated  software  for  the  data  acquisition 
system  and  contains  program  listings.  Together,  the  two  papers  describe  the 
entire  data  acquisition  system  developed  to  measure  flow  field  properties 
quickly  and  economically  in  wakes  using  the  seven-hole  probes.  , - 
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CALIBRATION  OF  SEVEN-HOLE  PROBES 
SUITABLE  FOR  HIGH  ANGLES  IN  SUBSONIC  COMPRESSIBLE  FLOWS 

A. A.  Gerner  and  C.L.  Maurer* 

v^/  Abstract 

This  paper  illustrates,  by  example,  a  method  for  calibrating  seven-hole  probes  to 
measure  local  total  and  static  pressures  and  relative  flow  angles  of  up  to  70  degrees 
in  subsonic  compressible  flows.  To  conserve  air  in  our  blowdown  wind  tunnel,  we  used 
the  method  of  Latin  Squares  to  statistically  sample  a  large  and  otherwise  unmanageable 
data  set,  thereby  reducing  to  a  minimum  the  number  of  data  points  required  to  construct 
a  polynomial  fit  to  the  data.  The  three-variable  third  order  polynomials  found  to  repre¬ 
sent  the  probe  calibration  permit  all  the  desired  output  quantities  to  be  found  explicit¬ 
ly  from  pressures  measured  on  the  probe  in  an  unknown  flow  field.  This  method  determines 
the  flow  angle  to  within  ±2  degrees  with  95  percent  certainty. 

’A 

I.  Introduction 

Many  present  and  future  aircraft  designs  are  beginning  to  employ  such  devices  as 
leading-edge  strakes,  forward-swept  wings,  and  canards.  These  devices  have  demon¬ 
strated  the  potential  for  enhancing  aircraft  maneuverability  and  control  by  producing 
strong  vortices.  Some  modern  aircraft,  such  as  the  Concorde,  actually  rely  on  vortices 
and  complex  flows  to  create  lift  However,  in  some  instances,  primary  and  secondary 
vortices  can  interact  unfavorably,  causing  separation  and  loss  of  lift  on  portions  of 
a  wing.  To  be  sure,  these  vortex-laden  flow  fields  are  quite  intricate  and  difficult 
to  analyze.  Flow  visualization  techniques  offer  a  way  to  gain  insight  into  vortex  inter¬ 
actions,  but  suffer  from  their  inability  to  provide  quantitative  information.  To  over¬ 
come  this  limitation,  small  probes  can  be  inserted  directly  into  the  flow  stream  to  gather 
meaningful  pressure  information.  Historically,  non-nulling  five-hole  probes  have  been 
used  to  determine  local  total  and  static  pressures  at  a  particular  point  in  a  flow,  as 
well  as  flow  directions  up  to  40  degrees  relative  to  their  axis.  Nevertheless,  it  is 
not  inconceivable  for  the  local  flow  angles  of  strong  vortices  to  exceed  60  degrees. 

For  this  reason,  the  Air  Force  Academy,  under  a  grant  from  NASA-Ames,  has  developed  a 
unique  seven-hole  probe.  In  addition  to  local  total  and  static  pressure  measurements, 
these  probes  have  demonstrated  the  ability  to  determine  flow  angles  up  to  80  degrees 
relative  to  their  axis.  When  combined  with  a  computerized  data  acquisition  system,  they 
are  capable  of  taking  data  at  a  rate  of  nearly  two  data  points  per  second,  much  faster 
than  nulling  probes  which  require  considerable  time  to  balance  probe  tip  pressures  be¬ 
fore  each  pressure  measurement  can  be  taken. 

In  addition,  these  probes  are  very  small  (about  one-tenth  of  an  inch  in  diameter), 
so  they  do  not  significantly  disturb  the  flow  they  are  measuring.  But  because  of  this 
small  size,  they  suffer  from  iuherent  manufacturing  defects.  As  a  result,  each  probe 
must  be  calibrated  before  it  can  be  a  useful  measuring  device.  Gallington  describes 
such  a  calibration  procedure  (Ref.  1)  which  is  both  fast  and  effective,  but  is  only  valid 
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for  incompressible  flows.  This  constraint  effectively  restricts  the  probes  to  surveying 
flows  with  free-stream  Mach  numbers  of  0.3  and  below  (Ref.  2).  However,  the  mere  presence 
of  a  body  in  an  airflow  causes  the  flow  to  accelerate  in  certain  regions,  resulting  in 
local  velocities  greater  chan  the  free-stream  velocity.  In  addition,  local  flow  velocities 
in  vortex  fields  are  likely  to  exceed  those  of  the  free  stream  too.  As  a  result,  even 
chough  the  free-stream  conditions  are  slow  enough  to  justify  an  incompressible  flow  as¬ 
sumption,  the  local  flow  conditions  might  in  fact  be  compressible.  The  purpose  of  this 
report  Chen  is  to  develop  a  power  series  calibration  scheme  which  accurately  determines 
the  actual  flow  conditions  from  pressures  measured  on  the  probe  in  compressible  flows. 

Because  of  the  similarities  in  compressible  and  incompressible  theory,  we  begin 
our  discussion  by  developing  fully  the  incompressible  calibration  theory.  Then,  by  anal¬ 
ogy,  we  expand  this  theory  to  obtain  the  desired  form  of  the  calibration  suitable  for 
compressible  flow.  Next,  we  describe  the  apparatus  and  procedures  we  used  to  calibrate 
a  probe  in  subsonic  compressible  flow,  and  finally,  we  discuss  the  results  of  that  cali¬ 
bration. 

II.  Incompressible  Flow  Theory 

To  simplify  our  discussion  on  the  calibration  of  seven-hole  probes  in  compressible 
flow,  we  will  begin  with  a  background  in  the  theory  of  incompressible  flow  calibrations. 
Much  of  the  information  presented  in  the  following  sections  is  drawn  from  the  work  of 
Gallington  (Ref.  1).  Readers  wishing  further  information  on  the  theory  of  incompressible 
probe  calibrations  are  highly  encouraged  to  consult  this  reference. 

Axis  System  for  Low  Flow  Angles 

To  begin  our  discussion  of  calibration  theory,  we  restrict  our  treatment  to  low 
flow  angles;  typically,  those  values  for  which  the  angle  between  the  velocity  vector 
of  the  flow  stream  and  the  probe's  axis  are  less  than  30  degrees.  The  more  familiar 
reference  system  for  low  flow  angles  measures  velocity  vectors  in  terms  of  the  angle 
of  attack,  a,  and  the  angle  of  sideslip,  B.  However,  in  choosing  our  reference  system, 
we  adopted  the  tangential  reference  system  illustrated  in  Figure  1.  In  this  system,  aT 
is  taken  to  be  the  projection  on  the  vertical  plane  of  the  angle  between  the  velocity 
vector  and  the  probe's  axis.  And  to  preserve  symmetry,  8^  is  defined  as  the  projection 
on  the  horizontal  plane  of  the  angle  between  the  probe's  axis  and  the  relative  wind. 

For  this  reason,  the  tangential  reference  system  differs  slightly  from  the  u-B  system, 
but  will  be  used  to  evaluate  the  low  angle  flow  properties. 

Pressure  Coefficients  for  Low  Flow  Angles 

For  low  flow  angles  it  is  desirable  to  define  dimensionless  pressure  coefficients 
which  utilize  all  seven  measured  probe  pressures  and  are  sensitive  to  changes  in  flow 
angularity  with  respect  to  the  probe's  x-axis.  From  Figure  2,  one  such  pressure  coeffi- 
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Figure  1.  Low  Angle  Reference  System 


cient  sensitive  to  changes  in  angle  of  attack  in  the  x-z  planes  is  defined  as: 

P..  -  Pi 

c  - - r — 

01  P7  -  Pi-s 


(1) 


where  the  numerator  measures  changes  in  flow  angularity  based  on  the  differences  in  oppo¬ 
site  port  pressures,  and  the  denominator  nondimens ionalizes  the  term  with  the  apparent 
dynamic  pressure.  This  pseudo-dynamic  pressure  is  obtained  from  the  difference  between 
the  central  port  pressure,  P7 ,  which  approximates  the  total  pressure  at  low  angles,  and 
the  average  of  the  six  surrounding  pressures,  Pi_* ,  which  collectively  approximates  the 
static  pressure.  From  the  definition  of  this  pressure  coefficient,  it  is  easy  to  see 
that  two  other  possibilities  also  exist:  one  which  measures  the  pressure  differential 
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Figure  2.  Port  Numbering  Convention  and  Principal  Axes 
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between  ports  three  and  six,  and  the  other,  which  measures  the  pressure  differential 
between  ports  two  and  five.  The  complete  set  of  these  pressure  coefficients  include: 
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But  before  these  coefficients  can  be  of  any  use  to  us,  they  must  be  resolved  into  the 
a,j,  -  8t  reference  system.  This  is  done  by  weighing  the  contribution  of  each  coefficient 
in  Eqns.  (2)  along  the  respective  axis,  which  results  in  the  following  two  equations: 
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The  first  equation  defining  Ca  contains  all  three  coefficients  of  Eqns.  (2).  In  particu¬ 
lar,  Cai  has  the  greatest  significance,  which  only  makes  sense  since  it  lies  directly 
along  the  axis  of  interest.  The  equation  defining  Cg  only  takes  into  account  the  last 
two  coefficients  of  Eqns.  (2),  assigning  to  each  an  equal  significance.  The  fact  that 
Cai  is  not  included  in  this  equation  again  makes  sense,  since  it  is  directly  aligned 
with  the  a^-direction,  ideally  making  it  insensitive  to  changes  in  the  perpendicular  in¬ 
direction.  In  summary,  the  procedure  of  obtaining  C0  and  Cg  requires  two  tasks.  First, 
determine  Cai ,  Ca2 ,  and  Ca3  from  the  seven  measured  pressures  using  Eqns.  (2)  and  then 
substitute  these  intermediate  quantities  into  Eqns.  (3)  for  the  desired  coefficients. 

Having  defined  the  two  angular  pressure  coefficients,  it  is  now  appropriate  to  dis¬ 
cuss  the  remaining  low  angle  pressure  coefficients,  CQ  and  Cq,  defined  as: 
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CQ  is  the  apparent  total  pressure  coefficient  and  functions  as  a  correction  factor 
to  convert  actual  pressures  measured  by  the  probe  to  accurate  values  of  local  total  pres¬ 
sure.  From  the  numerator,  it  is  seen  that  the  coefficient  measures  the  difference  between 
the  pseudo-total  pressure  measured  by  the  probe,  P7,  and  the  actual  total  pressure.  Just 
as  with  Eqn.  (1),  the  coefficient  is  nondlmensionalized  by  the  denominator,  which  is 
a  measure  of  the  apparent  dynamic  pressure  as  previously  described. 

Cq  functions  much  like  C0,  but  instead  of  correcting  probe  pressures  to  total  pres¬ 
sure,  Cq  relates  these  pressures  to  the  actual  dynamic  pressure.  In  this  coefficient, 
the  numerator  represents  the  probe's  approximation  of  dynamic  pressure  while  the  denomi¬ 
nator  consists  of  the  actual  dynamic  pressure. 
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Axis  System  for  High  Flow  Angles 

Up  to  this  point,  our  discussion  has  been  limited  to  a  description  of  the  pressure 
coefficients  used  for  flow  angles  below  30  degrees.  Yet,  the  real  advantage  of  using 
a  seven-hole  probe  lies  in  its  ability  to  determine  flow  angles  as  high  as  80  degrees 
to  the  probe's  x-axis. 

A  reference  system  better  suited  to  high  angle  measurement  than  the  tangential  sys¬ 
tem  is  the  polar  reference  system,  which  measures  flow  angularities  in  terms  of  8  and  <J> 
and  is  shown  in  Figure  3.  In  this  system,  6,  the  pitch  angle,  is  the  angle  the  velocity 
vector  makes  with  respect  to  the  probe's  x-axis;  and  #,  the  roll  angle,  describes  the 
azmuthal  orientation  of  the  velocity  vector  in  the  y-2  plane,  measured  counterclockwise 
from  the  negative  z-axis  as  viewed  from  the  front.  Although  a  singularity  exists  direct¬ 
ly  along  the  x-axis,  this  does  not  represent  any  problems  in  high  angle  measurement  and 
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Figure  3.  High  Angle  Reference  System 

is  avoided  entirely  by  switchihg  to  the  tangential  reference  system  at  low  angles. 

Pressure  Coefficients  for  High  Flow  Angles 

At  low  angles  of  attack,  all  seven  of  the  measured  probe  pressures  are  used  to  form 
the  pressure  coefficients.  However,  at  high  angles  of  attack,  as  illustrated  in  Figure 
the  flow  tends  to  detach  over  the  downstream  portions  of  the  probe.  Pressure  ports 
lying  in  this  separated  region  are  insensitive  to  changes  in  flow  angularity;  consequent¬ 
ly,  it  is  not  feasible  to  use  their  pressures  in  a  meaningful  coefficient.  As  a  result, 
the  pressure  coefficients  for  high  angle  measurements  must  be  defined  so  that  they  include 
only  the  pressures  from  ports  in  attached  flow. 

Typically,  the  separation  points  of  a  cylinder  in  turbulent  flow  are  over  100  degrees 
from  the  frontal  stagnation  point  (Ref.  3).  And  for  a  conical  body,  such  as  the  probe  tip, 


61 


US AFA-TR-8 1 -4 


Figure  4.  Flaw  Over  Probe  at  High  Angle  of  Attack 

the  flow  is  likely  to  remain  actached  longer.  In  addition,  the  u-velocity  component 
is  also  likely  to  extend  the  separation  points.  Thus,  for  the  condition  depicted  in 
Figure  4,  pressure  ports  three,  four,  five,  and  seven  lie  in  reliably  attached  flow; 
port  one  is  in  separated  flow;  and  the  disposition  of  ports  two  and  six  is  uncertain. 
Using  only  those  ports  in  attached  flow,  a  coefficient  sensitive  to  the  angle  of  pitch 
is  defined  as: 


P  *.  -  P  ; 


P-  - 


Pi  +  Ps 


(5) 


Following  the  same  rationale  as  in  the  low  angle  case,  the  numerator  measures  changes 
in  0  based  on  the  differences  in  opposing  port  pressures.  In  this  example,  P7,  the  smal¬ 
ler  of  the  two  pressures,  is  subtracted  from  P4,  Here  again,  the  coefficient  is  nondimen- 
sionallzed  by  dividing  through  with  the  apparent  dynamic  pressure.  This  pseudo-dynamic 
pressure  is  determined  from  the  difference  between  the  peripheral  port  pressure,  P4, 
which  at  high  angles  approximates  the  total  pressure,  and  the  average  of  P3  and  P5,  which 
when  taken  together  are  relatively  independent  to  changes  in  roll  and  approximate  the 
static  pressure. 

Using  a  similar  argument,  a  coefficient  which  changes  In  proportion  to  roll  angle 
is  appropriately  defined: 
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The  numerator  of  Eqn.  (6)  is  sensitive  to  changes  in  ♦,  in  that  as  the  velocity  vector 
rolls  in  either  direction,  the  windward  pressure  rises  and  the  leeward  pressure  falls. 

In  this  way,  the  difference  between  the  two  pressures  varies  significantly  for  variations 
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in  roll,  yet  the  average  of  their  sums  remains  relatively  constant.  Once  more,  the  coef¬ 
ficient  is  nondimensionalized  with  the  same  denominator  as  in  the  previous  case. 

Obviously,  the  above  two  coefficients  are  only  valid  for  a  narrow  range  in  roll 
angle  about  port  four.  That  is,  as  we  rotate  the  velocity  vector  to  either  side  of  port 
four,  the  region  of  separated  flow  approaches  either  port  five  or  three.  Therefore, 
to  insure  that  all  pressures  are  taken  from  ports  in  attached  flow,  we  restrict  Eqns.  (5) 
and  (6)  to  roughly  a  60-degree  pie-shaped  sector  centered  on  port  four.  In  this  way, 
six  pie-shaped  sectors  are  summarily  defined  for  high  angle  measurement,  such  that  each 
has  its  own  set  of  coefficients  based  on  the  pressures  in  attached  flow.  The  remaining 
angular  pressure  coefficients  are  defined  with  the  same  method  used  to  develop  Eqns. 

(5)  and  (6),  resulting  in  the  following  set  of  equations: 
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Similarly,  the  C0  and  Cq  coefficients  are  developed  with  the  same  rationale  used 
to  derive  their  low  angle  counterparts;  the  only  difference  resides  in  the  choice  of 
the  pressures  which  are  roughly  equivalent  to  total  and  static  pressures.  These  pressures, 
of  cou»"  ,  vary  in  relation  to  the  sector  a  particular  coefficient  describes.  The  cam¬ 
let  set  of  these  coefficients  include: 
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Division  of  Angular  Space 

Having  defined  a  host  of  coefficients  for  low  and  high  angles  and  for  various  sec¬ 
tors  around  the  probe,  a  question  arises  as  to  when  a  particular  set  of  coefficients 
should  be  used.  To  arbitrarily  assign  angular  cut-offs,  based  on  probe  symmetries,  would 
be  naive  since  actual  data  might  suggest  better  division  lines.  With  this  in  mind,  a 
better  scheme  for  locating  the  sector  division  lines  is  based  on  the  isobars  depicted 
in  Figure  5.  This  method  defines  seven  sectors,  the  central  low  angle  sector  and  the 
six  high  angle  periphery  sectors,  and  allocates  data  points  to  a  given  sector  based  on 
the  highest  port  pressure  measured  on  the  probe. 


Polynomial  Power  Series  Expansion 

Once  the  data  points  are  allocated  to  the  proper  sector  with  its  corresponding  pres¬ 
sure  coefficients,  a  fourth  order  polynomial  expansion  is  used  to  solve  for  the  desired 
quantities.  In  two  variables  (i.e.,  using  the  two  angular  pressure  coefficients)  this 
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SECTOR  1 


Figure  5.  Division  of  Angular  Space 
expansion  takes  on  the  following  form: 
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Si 


where  A  is  either  a^,  8^,,  CG  or  Cq  for  low  angles  and  0,  >  C0  or  Cq  for  high  angles, 

with  the  subscript  denoting  the  ith  such  quantity.  The  K's  are  the  calibration  coeffi¬ 
cients,  with  the  superscripts  denoting  the  quantity  to  which  a  particular  set  of 
K's  belong,  and  the  subscripts  identifying  the  coefficient  of  a  particular  term  in 
the  power  series  expansion.  Note  that  in  the  high  angle  case,  the  Ca's  and  Cg's  are 
replaced  by  Cg's  and  C^'s  respectively.  In  matrix  notation  for  "n"  data  points  of  a 
particular  sector,  a  set  of  Eqn.  (9)'s  are  represented  as: 
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But  to  simplify  further  discussion,  Eqn.  (10)  is  rewritten  as: 

[A]  -  [C] [K]  (11) 

In  this  form,  the  n  x  1  A-matrix  contains  n  values  of  the  flow  parameters  of  interest, 
the  n  x  15  C-matrix  contains  the  expanded  angular  pressure  coefficients  for  each  of  the 
n  flow  parameters,  and  the  15  x  1  K-matrix  contains  the  calibration  coefficients  for 
the  particular  flow  parameter  of  interest. 

Determining  the  Calibration  Coefficients 

During  the  calibration  process,  the  quantities  within  the  A-matrix  are  set  up  by 
the  experimenter  in  terms  of  the  known  tunnel  conditions,  and  the  terms  within  the  C- 
matrix  are  determined  from  the  measured  probe  pressures.  A  calibration  procedure,  there¬ 
fore,  involves  the  calculation  of  the  unknown  K-matrix.  This  calculation  is  performed 
by  rearranging  Eqn.  (11)  to  solve  for  the  unknown  calibration  coefficients.  With  matrix 
algebra,  this  is  performed  according  to  the  procedure  outlined  by  Netter  and  Vasserman 
(Ref.  4): 

First  multiply  each  side  of  Eqn.  (11)  by  the  transpose  of  the  C-matrix: 

[C]T  [A]  -  [C]T  [C] [K]  -  [CTCJ[K]  (12) 

Next  multiply  each  side  by  the  inverse  of  the  recently  created  CTC-matrix: 

[CTC]"1[C]T[A]  -  [CTC]_1[CTC][K]  (13) 

Realizing  that  the  product  of  a  matrix  and  its  inverse  results  in  the  identity  matrix, 
Eqn.  (13)  simplifies  to  yield  a  solution  for  the  unknown  K-matrix  in  terms  of  the  known 
C-  and  A-matrlces: 
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(K]  -  [CTC]-1[C]T{A]  <U> 

This  technique  determines  the  calibration  coefficients  by  a  least  squares  curve  fit  to 
the  experimental  data. 


Determining  the  Desired  Flow  Properties 

Once  the  calibration  coefficients  are  determined,  the  calibration  process  is  complete 
The  probe  is  ready  to  be  inserted  in  an  unknown  flow  field  and  the  desired  flow  properties 
determined.  Once  in  the  flow  field,  the  probe's  measured  pressure  readings  allow  us 
to  determine  the  angular  pressure  coefficients;  these  coefficients  are  then  manipulated 
to  fill  the  C-matrix  of  Eqn.  (10).  Since  the  K-matrix  is  already  known,  the  desired 
flow  properties  in  the  A-matrix  are  then  determined  explicitly.  For  a  particular  condi¬ 
tion,  the  solutions  for  the  desired  flow  properties  take  on  the  following  functional  forms 


Inner  Sector  (low  flow  angles) 

aT  °T  aT 

aT  -  f(Ca,  Cg)  =  Kx  +  K2  Ca  +  K3  ACg  +  . 

8t  BT 
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c  c  c 
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C  C  C 

C  -  f(C  ,  Cft)  -  Kx  q  +  K2  qC  +  Kj  %  +  . 

q  a  p  a  p 

Outer  Sectors  (high  flow  angles) 
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f(Cen.  c*n)-  Kx  n  +  k2  nc9n  +  K,  %n  +  . 

Cq 
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(16) 


Once  the  local  total  and  dynamic  pressure  coefficients  are  specified,  it  is  possible  to 
determine  the  local  total  and  dynamic  pressures.  This  is  accomplished  by  rearranging  Eqns 
(4)  in  the  case  of  the  inner  sector,  and  Eqns.  (8)  in  the  case  of  the  outer  sectors, 
thereby  solving  for  the  desired  pressures.  This  calculation  involves  the  polynomial 
result  for  CQ  or  C^  and  the  seven  measured  probe  pressures.  As  ah  example,  the  inner 
sector  equations  for  local  total  and  dynamic  pressures  are  derived  as  follows: 
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Recalling  that 
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(4) 


the  local  total  and  dynamic  pressures  are  then  solved  by  manipulating  Eqns.  (4);  thus: 


P  T  *  P?  -  C  <P7  -  Pi-6) 

oL  o 


(17) 


The  same  procedure  is  extended  to  determine  the  local  total  and  dynamic  pressures  in 
the  outer  sectors. 

Ill.  Extension  to  Compressible  Flow 

Up  to  this  point  we  have  discussed  the  form  of  a  calibration  procedure  which  is 
only  valid  for  incompressible  flow.  This  is  so  because  the  calibration  coefficients 
depend  directly  upon  the  angular  pressure  coefficients,  which  as  pressure  coefficients 
are  themselves  dependent  upon  Mach  number  (Ref.  5).  As  such,  seven-hole  probes  with 
the  present  method  of  calibration  are  limited  to  surveying  flows  within  the  incompres¬ 
sible  regime.  This  restriction  is  lifted  easily  enough  with  the  inclusion  of  an  addi¬ 
tional  pressure  coefficient  representative  of  compressibility.  Consequently,  the  number 
of  variables  in  the  polynomial  curve  fit  is  increased  from  two  (for  incompressible  flow) 
to  three  (for  compressible  flow). 

Determining  a  Coefficient  of  Compressibility 

Just  like  the  angular  pressure  coefficients  (C  and  C„  or  C.  and  C.),  the  compres- 

a  p  h  9 

sibility  coefficient  (hereafter  denoted  C^)  must  be  nondimens ional  and  determined  strict¬ 
ly  from  pressures  measured  on  the  probe.  In  addition,  Cjj  must  have  the  feature  that 
it  approaches  zero  at  very  low  Mach  numbers.  That  is,  at  very  low  Mach  numbers  all  terms 
bearing  a  compressibility  coefficient  should  be  insignificant  in  the  power  series  expan¬ 
sion.  This,  in  effect,  brings  us  back  to  an  incompressible  power  series  expansion  in 
two  variables,  where  the  significant  terms  only  contain  the  angular  pressure  coefficients. 
A  further  constraint  requires  the  compressibility  coefficient  to  approach  a  finite  value 
in  the  hypersonic  limit.  In  other  words,  for  very  high  Mach  numbers,  large  changes  in 
Mach  number  should  have  a  neglible  effect  on  the  compreasibility  coefficient.  This  re¬ 
flects  a  limitation  shared  by  all  pressure  probe  methods  in  that  Mach  number  becomes 
Indeterminate  at  hypersonic  speeds  (Ref.  1).  These  requirements  are  satisfied  by  model¬ 
ing  the  compressibility  coefficient  after  the  dynamic  to  total  pressure  ratio  depicted 
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in  Figure  6.  This  figure  is  idealized  for  isentropic  flow;  as  such,  it  is  not  useful 
for  flow  speeds  much  beyond  Mach  one.  However,  since  our  calibration  is  limited  to  sub¬ 
sonic  speeds,  the  isentropic  idealization  is  valid  through  the  sonic  limit. 


MACH  NUMBER  (M) 


Figure  6.  Compressibility  Effects  as  a  Function  of  Mach  Number 

To  develop  a  compressibility  coefficient  from  probe  measured  pressures,  we  need 
to  represent  the  total  and  dynamic  pressures  from  probe  measured  pressures.  For  the 
low  angle  case,  a  pseudo-total  pressure  is  approximated  by  P7  with  a  pseudo-dynamic  pres¬ 
sure  approximated  by  P7  -  Pj-6.  Thus,  for  the  inner  sector,  a  compressibility  coefficient 
in  terms  of  the  dynamic  to  total  pressure  ratio  is  modeled  by: 


°M7 


(18) 


The  compressibility  coefficients  for  each  of  the  outer  sectors  are  modeled  in  a  similar 
manner;  these  outer  sector  coefficients  include: 


Ps  +  Pz 
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Selection  of  Data  Points 

Typical  incompressible  probe  calibrations  take  approximately  80  data  points  in  two 
variables  (Ca  and  Cg)  for  each  of  the  seven  sectors  (Ref.  1).  This  results  in  a  total 
of  about  560  data  points  for  a  complete  calibration.  To  extend  this  present  scheme  into 
yet  another  dimension  (Mach  number)  would  create  a  data  set  of  intractable  proportions. 
Consequently,  it  is  necessary  to  represent  the  data  set  with  a  sample  of  more  manageable 
proportions.  In  addition,  this  sample  must  be  chosen  such  that  the  density  of  chosen 
data  points  throughout  the  data  set  is  homogenous.  In  other  words,  the  sample  must  be 
an  accurate  representation  of  the  data  set,  otherwise  the  calibration  routine  will  not 
offer  consistent  accuracy  throughout  the  range  of  data. 

A  method  of  ensuring  a  homogenous,  yet  random,  sample  of  a  three-dimensional  parame¬ 
ter  space  is  suggested  by  Cochran  and  Cox  (Ref.  6).  The  technique  is  known  as  the  method 
of  Latin  Squares,  an  example  of  which  is  shown  in  Figure  7.  The  plan  depicted  in  this 
figure  is  a  3  x  3  square  which  actually  represents  a  three-dimensional  parameter  space; 
one  variable  along  the  vertical  axis,  a  second  variable  along  the  horizontal  axis,  and 
the  third  variable  denoted  by  the  letters  along  the  axis  going  into  the  page.  This 
square  is  better  visualized  in  Figure  8,  where  it  is  shown  in  three  dimensions  instead 
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Figure  7.  3x3  Latin  Square  Plan 
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Figure  8.  3x3  LaCin  Square  Shovm  in  Three  Dimensions 

of  the  two-dimensional  rendition  of  Figure  7.  The  3x3x3  cube  represents  the  entire 
parameter  space,  and  each  sub-cube  represents  a  discrete  data  point  within  the  entire 
data  set.  In  Figure  7,  the  A's  denote  the  data  points  to  be  selected  in  the  A-plane 
(i.e.,  the  plane  in  which  the  third  variable,  Mach  number,  is  held  constant),  the  B's 
denote  the  points  selected  in  the  B-plane,  and  so  on.  In  this  way,  the  entire  data 
set  is  sampled  by  the  points  which  appear  as  the  shaded  cubes  in  Figure  8.  The  unique 
feature  of  this  sampling  technique  is  that  no  matter  from  which  axis  direction  the  cube 
is  viewed,  the  entire  frontal  space  will  appear  covered  by  data  points,  one  point  deep; 
that  is,  no  two  points  will  appear  to  overlap  each  other.  Furthermore,  in  any  given 
plane  where  one  variable  Is  held  constant,  every  value  of  the  remaining  two  variables 
is  sampled  exactly  once.  Consequently,  the  method  of  Latin  Squares  guarantees  a  homo¬ 
genous  density  of  data  points  comprising  the  sample  of  the  representative  data  set.  Of 
course,  the  technique  of  Latin  Squares  is  not  limited  to  a  3  x  3  plan,  but  can  be  ex¬ 
panded  to  a  12  x  12  plan  if  desired;  however,  the  more  common  plans  range  between  the 
5x5  and  8x8  squares  (Ref.  6).  This  is  true  because  as  plan  size  increases,  the  ratio 
of  points  sampled  to  total  points  in  the  data  set  decreases.  To  illustrate,  the  3x3 
square  of  Figure  7  samples  a  27-point  data  set  with  nine  points,  whereas  a  10  x  10  square 
uses  100  points  to  sample  a  1000-point  set.  Accordingly,  the  3x3  Latin  Square  samples 
33  percent  of  the  entire  data  set,  while  the  10  x  10  Latin  Square  samples  only  10  percent 
of  the  complete  set.  Nevertheless,  a  3  x  3  square  describes  a  given  parameter  apace 
with  only  27  data  points,  while  a  10  x  10  square  divides  the  same  parameter  space  into 
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1000  discrete  da ca  points.  Thus,  a  larger  Latin  Square  represents  a  given  parameter 
space  with  greater  resolution,  but  saiup  [es  the  resulting  data  set  with  a  smaller  percen¬ 
tage  of  points.  The  experimenter,  therefore,  is  faced  with  a  compromise  in  choosing 
the  Latin  Square  best  suited  to  his  needs. 

For  reasons  to  be  discussed  later,  we  selected  a  6  x  6  Latin  Square  (shown  in  Figure 
9)  for  each  sector.  This  gives  us  a  daca  set  with  six  different  values  for  each  of  the 
three  variables  0,  $,  and  M,  for  a  total  of  216  daca  points  per  sector.  Using  the  sampling 
technique  of  Latin  Squares,  each  216-point  data  set  is  represented  by  a  sample  of  36 
points.  As  such,  a  total  of  252  data  points  will  be  experimentally  tested  in  all  seven 
sectors  for  a  complete  calibration  of  the  probe. 
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Figure  9.  6x6  Latin  Square  Plan 

Polynomial  Power  Series  Expansion  in  Three  Variables 

In  the  incompressible  calibration,  a  fourth  order  polynomial  expansion  in  two  vari¬ 
ables  was  used  for  a  tota'  of  15  terms  with  15  corresponding  calibration  coefficients. 

By  adding  a  third  variable  to  the  calibiation,  the  number  of  terms  in  the  fourth  order 
expansion  jumps  from  15  to  35  with  35  corresponding,  calibration  coefficients.  In  order 
to  obtain  a  valid  estimation  of  standard  deviation,  a  surplus  of  about  20  data  points 
over  the  number  of  calibra-ion  coefficients  is  necessary  (Ref.  1).  This  sets  the  required 
number  of  data  points  for  a  given  sector  at  55;  however,  our  sample  size  has  already 
been  constrained  to  36  points  due  to  the  selection  of  a  6  x  6  Latin  Square.  It  is  possi¬ 
ble  to  take  more  than  one  36-point  Latin  Square  sampling  within  the  216-point  data  set, 
but  in  the  Interest  of  keeping  the  total  number  of  points  required  for  a  complete  cali¬ 
bration  down  to  a  manageable  number,  so  that  time  spent  in  the  wind  tunnel  can  be  mini¬ 
mized,  we  elect  to  limit  each  of  the  seven  samples  to  36  data  points.  As  a  result,  a 
fourth  order  curve  fit  is  no  longer  feasible;  consequently,  the  polynomial  expansion  is 
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reduced  to  the  next  lowest  order. 

A  third  order  polynomial  expansion  in  three  variables  requires  20  calibration  coef¬ 
ficients.  This  leaves  us  with  «  16  data  point  surplus,  which  is  sufficient  to  calculate 
a  valid  standard  deviation.  Using  the  same  format  at  Eqn.  (9),  a  third  order  expansion 
in  three  variables  takes  on  the  following  form: 


A±  -  Ki  + 


Order  of  Teraa 
0th 


K2Cai  +  K,C8l  +  K4CMl  +  1st 

KsC^  +  K^C8l  +  +  KtcaiCei  +  K&ai  CMi  +  KtoC8lCMi  +  2nd 

KilC^  +  K^cgj  +  KiaC^  +  ^ 

Ki«CaiC8i  +  K^C8iCMi  +  ^8CaiC^  +  K^Cg^  +  K^oC^Cg^ 


(20) 


where  A  is  either  a^,,  8T,  COI  or  Cq  for  the  inner  sector  and  8,  <f>,  C0,  or  Cq  tor  the 
outer  sectors,  with  the  subscript  denoting  the  ith  such  quantity.  The  K's  are  the  cali¬ 
bration  coefficients,  where  the  superscript  denotes  the  quantity  to  which  a  particular 
set  of  K's  belong,  and  the  subscripts  identify  the  coefficient  of  a  particular  term  in 
the  power  series  expansion.  In  the  case  of  an  outer  sector,  the  Ca's  and  Cg's  are  re¬ 
placed  by  C$'s  and  Ccji's,  respectively. 


Determination  of  Mach  Number 

As  stated  previously,  this  experiment  limits  its  scope  to  surveying  Mach  numbers 
slightly  below  sonic  flow  on  down  to  incompressible  flow.  Consequently,  there  are  no 
shocks  ahead  of  the  probe  and  the  isentropic  flow  relation  applies: 


P 

o 

Rearranging  Eqn.  (21)  and  deriving  Mach  number  in  terms  of  the  dynamic  to  total  pressure 
yields : 
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And  since  we  know  C0,  Cq,  and  the  seven  probe  pressures,  we  can  explicitly  determine 
the  local  Mach  number  from  the  dynamic  to  total  pressure  ratio.  For  the  inner  sector, 
this  ratio  is  determined  from  Eqns.  (4)  and  the  seven  probe  pressures: 
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Analaguus  filiations  for  the  outer  sectors  ate  developed  in  the  same  way,  except  that 
Eqns.  (8  an  1  the  appropriate  probe  piessures  are  used  instead.  During  calibration, 
tne  actual  values  of  Eqn.  id)  and  its  outer  sector  counterparts  are  known  from  the  mea¬ 
sured  tunnel  conditions.  Solving  Eqn.  (21)  in  terms  of  the  dynamic  to  total  pressure 
ratio  as  a  function  of  Mach  number  gives  us: 
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Thus,  if  tie  in',  and  static  pressures  cannot  be  measured  directly,  we  can  still  deter¬ 
mine  the  dvn  !■  total  pressure  racio  if  only  the  Mach  number  is  known.  By  noting 

the  dlffereu  .tween  Eqns.  (1~)  and  (24),  we  can  estimate  how  accurately  our  polynomial 
curve  fit  deconn.oes  the  Mach  number  for  the  inner  sector;  similar  arguments  are  also 
extended  to  the  outer  sectors. 


Estimating  Accura cy  of  the  Curve  Fit 

Providing  there  are  approximately  twenty  more  data  points  than  the  number  of  calibra¬ 
tion  coefficients.  It  is  statistically  feasible  to  calculate  a  global  estimate  of  the 
accuracy  of  the  curve  fit  for  each  of  the  flow  parameters.  This  is  done  by  computing 
the  standard  deviation  of  the  difference  between  the  experimental  data  and  the  polynomial 
prediction  of  that  daca.  For  flow  angles  and  Mach  number,  the  following  relation  applies 
(Ref.  7): 


d(A) 
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where  a  is  the  standard  deviation;  n,  the  total  number  of  data  points;  and  A,  the  desired 
flow  parameter. 

Even  though  the  total  and  dynamic  pressures  are  determined  from  C0  and  Cq.  the  ac¬ 
curacies  of  these  pressures  are  not  representative  of  the  accuracies  obtained  for  CQ 
and  Cq.  To  estimate  the  uncertainty  in  determining  the  total  and  dynamic  pressures  (for 
the  inner  sector,  for  example)  requires  the  defining  equations,  Eqns.  (17),  and  applying 
to  them  the  method  of  Kline  and  McClincock  (Ref.  8).  Then,  by  taking  the  standard  devia¬ 
tion  of  these  uncertainties  and  nondimeusionalizing  them  with  the  dynamic  pressure,  we 
arrive  at  the  following  estimates  for  the  accuracies  of  these  pressures: 
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where  Cq  is  the  average  value  of  this  coefficient  for  a  given  sector  denoted  by  "n". 

Providing  the  errors  between  the  actual  values  and  the  polynomial  predictions  of 
those  values  are  normally  distributed,  there  is  a  68.3  percent  probability  that  the  poly¬ 
nomial  prediction  will  fall  within  one  standard  deviation  of  the  actual  value.  And  at 
1.96a,  this  certainty  is  increased  to  95  percent. 


IV.  Apparatus  and  Calibration  Procedure 
Probe  Geometry 

The  seven-hole  probe  was  constructed  at  the  Air  Force  Academy  by  packing  seven 
properly-sized  stainless  steel  tubes  into  a  larger  stainless  steel  tube  as  shown  in  Figure 
10.  Although  the  inner  seven  tubes  have  an  outside  diameter  of  only  0.028  inches  with 
a  .005-inch  wall  thickness,  accurate  alignment  is  reasonably  insured  since  the  tubes 
can  only  pack  in  one  unique  way.  Once  assembled  in  this  arrangement,  the  tubes  are  then 
soldered  together  and  machined  to  provide  the  25-degree  half  angle  at  the  tip. 


.018  IN. 

DIAMETER 


Figure  10.  Probe  Geometry 


Probe  Mounting  Hardware 

The  test  facility  we  used  has  a  variable  position  sector  capable  of  traversing  25 
degrees  in  either  direction  for  a  total  sweep  of  50  degrees.  Yet,  for  a  complete  cali¬ 
bration  of  the  probe,  a  range  of  0  through  80  degrees  is  necessary.  Consequently,  we 
constructed  two  stings  (see  Figure  11):  a  15-degree  bent  sting  for  low  angle  measurements 
from  -10  to  A0  degrees  angle'  of  pitch,  and  a  55-degree  bent  sting  for  high  angle  measure¬ 
ments  from  30  to  80  degrees  angle  of  pitch.  In  addition,  both  stings  were  designed  to 
permit  the  tip  of  the  probe  to  pivot  about  a  fixed  point  in  the  center  of  the  tunnel. 
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Figure  11.  Scing  Geometry 


This  insured  uniform  flow  over  the  tip,  despite  changes  in  the  angle  of  pitch. 

Referring  to  Figure  12,  each  sting  has  36  holes  drilled  in  the  front  face.  These 
holes  are  evenly  spaced  at  10-degree  intervals  and  circumscribe  a  complete  circle  in 
roll.  In  this  way,  roll  angle  is  accurately  set  by  engaging  the  alignment  pin  on  the 
probe  holder  with  the  alignment  hole  on  the  sting's  face. 


Data  Point  Selection 

Since  there  are  a  total  of  36  discrete  roll  angles,  we  are  allowed  to  test  six  dif¬ 
ferent  roll  angles  in  each  outer  sector.  This  is  the  primary  reason  why  we  chose  the 
6x6  Latin  Square.  Although  the  technique  of  Latin  Squares  allows  us  to  conveniently 
sample  large  data  sets,  there  are  some  drawbacks  associated  with  seven-hole  probe  cali¬ 
bration.  Primarily,  we  can  no  longer  test  the  entire  range  of  data  and  then  allocate 
the  data  to  a  given  sector  based  on  the  highest  port  pressure.  Instead,  we  must  deter¬ 
mine  beforehand  which  points  will  be  tested  and  the  sector  they  belong  to.  As  a  result, 
we  are  forced  to  draw  the  angular  cutoffs  for  each  sector  without  knowing  where  they 
actually  lie. 

In  terms  of  angle  of  pitch,  experience  with  incompressible  calibrations  has  shown 
that  30  degrees  is  the  smallest  angle  at  which  almost  all  points  still  fall  in  the  outer 
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INCIDENCE  SETTING 


Figure  12.  Probe,  Holder,  and  Sting  Assembly 

sectors.  Consequently,  we  selected  30  degrees  as  the  lower  limit  of  high  angle  measure¬ 
ment.  And  to  match  the  number  of  roll  measurements  with  an  equal  number  of  pitch  settings, 
we  measured  pitch  from  30  to  80  degrees  in  10-degree  increments  for  the  outer  sectors. 
Specifying  the  roll  angles  for  each  sector,  however,  presented  some  difficulties.  Assum¬ 
ing  a  perfectly  constructed  probe,  the  outer  sector  boundaries,  formed  by  the  isobars 
between  two  adjacent  peripheral  ports,  are  coincident  with  a  roll  angle  on  which  data 
is  taken.  This  predicament  is  illustrated  in  Figure  13a  for  sector  four.  In  this  exam¬ 
ple,  it  is  uncertain  whether  the  points  along  the  =  330  degrees  line  belong  in  sector 
four  or  five,  or  whether  the  points  on  the  =  30  degrees  line  belong  in  sector  three 
or  four.  To  resolve  this  problem,  we  rotated  the  probe  5  degrees  and  then  permanently 
fixed  it  to  the  holder  in  this  position.  From  Figure  13b  we  see  that  the  isobars  have 
shifted  off  the  points  and  that  each  sector  has  six  clearly  defined  roll  angles.  Realis¬ 
tically,  the  probe  is  not  constructed  perfectly,  causing  the  isobars  to  deviate  from 
their  ideal  positions.  But  even  so,  the  isobars  can  deviate  up  to  5  degrees  in  either 
direction  before  the  points  of  one  sector  fall  in  another  sector.  At  this  point,  one 
might  ask  if  the  probe  can  be  accurately  calibrated  with  the  given  offset.  In  actuality, 
the  offset  makes  no  difference.  Since  only  the  probe  was  rotated,  leaving  the  index 
mark  of  the  probe  unchanged,  the  calibration  process  will  interpret  the  pressures  in 
terms  of  the  roll  angle  set  by  the  index  mark.  An  alternative  to  the  offset  technique 
would  use  the  roll  angles  between  those  presently  measured  (i.e.,  5,  15,  25  .  .  .  355 
degrees);  however,  this  is  an  option  we  chose  not  to  take. 

Now  that  we've  determined  the  angular  coordinates  of  the  data  points  to  be  taken  in 
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A.  Probe  and  Holder 
Before  5-llegree  Offset 

Figure  13.  Probe  and  Holder  Before 

the  outer  sectors,  it's  time  we  considered  Che 


B.  Frohe  and  Holder 
After  5-Degree  Offset 

and  After  Offset  for  Sector  Four 


inner  sector.  Since  wo  cannot  evenly 
represent  the  inner  sector  with  a  square  matrix  in  -  and  r,  we  simply  covered  the  parame¬ 
ter  space  with  an  even  densicy  of  36  data  points.  The  values  of  pitch  angle  were  8, 

16,  and  24  degrees,  with  roll  angles  selected  to  insure  a  uniform  distribution.  The  6x6 
Latin  Square  matrix  was  then  systematically  filled  with  these  ordered  pairs.  Although 
data  were  taken  in  terms  of  pitch  and  roll,  these  angles  weie  converted  to  angles  of  attack 
and  sideslip  prior  to  calibration  calculations. 


Figure  14.  Distribution  of  Data  Points 
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Since  the  angular  coordinates  of  all  data  points  are  known,  the  last  task  is  to 
determine  the  Mach  number  at  which  each  is  sampled.  The  Mach  numbers  for  each  of  the 
six  letters  in  the  Latin  Square  of  Figure  9  are:  A,  M  =  0.37;  B,  M  =  0.45;  C,  M  =  0.53; 

D,  M  =  0.66;  E,  M  =  0.77;  F,  M  =  0.91.  The  resulting  distribution  of  data  points,  and 
the  Mach  number  each  was  tested  at,  can  be  inferred  from  Figure  14,  which  depicts  the 
inner  sector  and  a  typical  outer  sector. 

Test  Procedure  and  Software 

The  experiment  was  conducted  in  the  1  x  1  foot  blowdown  wind  tunnel  at  the  United 
States  Air  Force  Academy.  All  pressure  measurements  were  made  with  a  Scanivalve  Corpor¬ 
ation  Model  T  Scanivalve  and  pressure  transducer,  calibrated  by  a  Wallace  and  Tiernan 
sonar  manometer  digital  U-tube.  Data  were  collected  and  reduced  by  a  Digital  PDP  11/45 
computer  and  LPS-11  Laboratory  Peripheral  System. 

Before  any  data  were  taken,  the  test  conditions  had  to  be  set  for  each  run.  Angle 
of  pitch  was  set  first  with  a  Gunner's  quadrant.  The  quadrant  was  clamped  to  the  sting's 
incidence  setting  block  (refer  to  Figure  12)  and  set  to  the  desired  angle  of  pitch  minus 
the  angle  of  the  elbow  on  the  sting  (either  15  or  55  degrees).  Next,  the  sector  was 
hydraulically  positioned  to  level  the  quadrant.  Since  the  sector  normally  has  a  tendency 
tc  drift  between  runs,  we  froze  the  mechanism  by  driving  it  against  blocks.  Shims  were 
used  to  fine  tune  the  adjustments.  Using  this  technique,  the  pitch  was  set  to  within 
±0.0112  degrees  of  the  desired  angle,  and  without  drift.  Once  set,  the  pitch  was  not 
changed  until  all  points  at  that  angle  were  taken.  The  roll  angle  was  set  by  pulling 
the  c;)ring-  loaded  probe  forward  to  clear  the  alignment  pin,  then  rotating  it  to  the  de¬ 
sired  setting  and  engaging  it  into  the  new  position.  For  a  complete  360-degree  revolu¬ 
tion,  the  probe  was  first  rotated  counterclockwise  (as  viewed  from  the  front)  from  0 
to  180  degrees  in  10-degree  increments.  Next,  it  was  rotated  clockwise  to  the  190-degree 

position.  And  finally,  the  probe  was  rotated  counterclockwise  from  190  to  350  degrees 

in  10-degree  increments.  This  subjected  the  tubes  to  a  maximum  twist  of  180  degrees. 

After  setting  the  roll  angle,  the  tunnel  was  closed  by  a  hydraulic  ram  and  bolted  shut. 
Mach  number  was  then  set  by  manually  adjusting  the  tunnel's  inner  geometry  with  a  series 
of  hand  cranks  on  the  outside  of  the  tunnel. 

Prior  to  each  run,  the  above  variables  were  verified  and  input  into  the  computer 
by  a  program  titled  TRISHP.  Once  the  input  conditions  were  set,  air  was  blown,  and  after 
a  few  seconds  to  acquire  steady-state  conditions,  data  were  recorded.  Following  this, 
the  tunnel  was  reopened  and  the  process  repeated  for  a  new  data  point.  Typically,  each 

run  took  about  ten  minutes  to  set  up.  Consequently,  tests  ran  for  over  a  week  to  acquire 

all  252  data  points. 

After  all  data  were  taken,  a  second  program  titled  TRICAL  performed  the  matrix  opera¬ 
tions,  determined  the  calibration  coefficients,  and  estimated  the  accuracy  of  the  poly¬ 
nomial  expansions  in  fitting  the  known  data.  Having  done  this,  the  probe  calibration 


79 


USAFA-TR-81-4 


is  complete, 

V.  Discussion 

The  primary  purpose  of  this  experiment  is  to  create  a  power  series  curve  fit  which 
accurately  determines  the  actual  flow  conditions  from  pressures  measured  on  the  probe. 

The  best  way  to  evaluate  the  performance  of  the  curve  fit  is  to  analyze  the  standard 
deviations  between  the  experimental  data  and  the  polynomial  determination  of  those  data. 
These  standard  deviations  are  presented  in  Table  1.  in  addition,  the  standard  deviations 
obtained  from  past  incompressible  calibrations  are  included  in  Table  2  for  comparison 
(Ref.  1). 

From  Table  1,  the  standard  deviations  of  and  8^,  are  0.78  degrees  and  0.72  degrees 
respectively.  Although  these  values  appear  to  be  quite  good,  they  are  significantly 

Table  1 

STANDARD  DEVIATIONS 
COMPRESSIBLE  FLOW  CALIBRATION 


INNER  SECTOR 

OUTER  SECTORS 

S  ECTOR^-''''''^ 
EXPRESSION 

m 

AVERAGE 

1  -  6 

°(aT) 

0.78° 

0(6) 

4.27° 

d(6x) 

0.72° 

0(f>) 

0.57° 

°(PoL> 

2.57, 

0(POL> 

5.7% 

PoL  "  P“>L 

PoL  '  P»L 

1 

o(p  -  p  ) 

v  oL  “L 

1.4% 

0(POL  "  P-L> 

.  c.  2% 

P  “  P 
oL 

P  -  P 
oL  °°L 

c(M) 

0.006 

O(M) 

0.061 

higher  than  those  obtained  from  the  incompressible  calibrations.  Specifically,  the  in¬ 
compressible  flow  calibrations  produced  standard  deviations  of  0.42  degrees  for  and 
0.36  degrees  for  8T.  Similar  variations  are  also  noted  in  the  percent  errors  of  total 
and  dynamic  pressures.  Yet,  when  we  recall  that  the  Incompressible  calibrations  used 
a  fourth  order  power  series  with  80  data  points  per  sector  and  the  compressible  calibra¬ 
tion  used  a  third  order  power  series  with  only  36  data  points  per  sector,  the  slight 
differences  are  quite  understandable.  Still  another  measure  of  curve-fit  accuracy 
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in  the  inner  sector  is  apparent  from  Figures  15  and  16.  Both  of  these  figures  depict 
isolines  of  Ca  and  Cg  plotted  against  and  with  held  constant.  Figure  15  pre' 
sents  this  information  for  a  low  Mach  number,  with  Figure  16  representative  of  a  high 
Mach  number.  In  both  cases,  the  iso  lines  of  the  coefficients  are  nearly  straight  and 
relatively  orthogonal  to  each  other,  implying  a  linear  dependence  on  their  respective 
angle  and  independence  to  the  other  angle.  These  properties,  however,  begin  to  break 
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Figure  16.  Isolines  of  Ca  and  Ca,  versus  and  8^  for  Low  Angles  at  High  Mach  Number 


down  as  we  exceed  the  limits  of  our  experimental  data.  But  since  the  boundaries  of  the 
experimental  data  coincide  with  the  inner-outer  sector  interface  (see  Figure  17),  data 
points  lying  outside  the  inner  sector's  experimental  boundary  will  fall  into  the  outer 
sectors.  Consequently,  tne  accuracy  in  determining  flow  properties  is  uniform  through¬ 
out  the  inner  sector. 

In  the  outer  sectors,  the  average  standard  deviation  of  the  error  in  calculating  ♦ 
is  0.57  degrees  for  the  compressible  calibration  and  1.17  degrees  for  previous  incompres 
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Figure  17.  Ideal  and  Actual  Sector  Boundaries  Based  on 
the  Highest  Probe  Measured  Pressures 

sible  calibrations.  Judging  from  Figure  17,  one  might  expect  the  average  standard  devi¬ 
ation  of  <(>  for  the  compressible  calibration  to  be  substantially  greater  than  that  for 
the  incompressible  calibrations,  because  the  actual  boundaries  of  sectors  one  through 
four  do  not  coincide  with  the  ideal  boundaries  for  which  the  calibration  was  made.  Al¬ 
though  this  phenomenon  is  entirely  consistent  with  the  manufacturing  anomalies  associated 
with  probe  construction,  the  actual  boundaries  cannot  be  determined  in  advance  and,  there¬ 
fore,  cannot  be  taken  into  account  in  a  calibration  scheme  using  the  method  of  Latin 
Squares.  Nevertheless,  despite  the  actual  locations  of  the  sector  boundaries,  the  stand¬ 
ard  deviation  in  $  for  the  compressible  calibration  agrees  favorably  with  its  incompres¬ 
sible  counterpart.  The  reason  for  this  is  shown  in  Figure  18.  Even  though  Figure  18 
is  based  on  data  taken  from  incompressible  calibrations,  it  is  representative  of  seven- 
hole  probes  in  general.  That  is,  the  coefficient  of  roll  continues  to  behave  linearly, 
or  in  a  way  that  can  be  represented  easily  by  a  polynomial,  slightly  beyond  the  isobaric 
sector  boundaries.  As  such,  the  effect  of  using  data  within  10  degrees  in  roll  from 
the  selected  boundaries,  as  Figure  18  suggests,  does  not  have  an  adverse  effect  on  fitting 
the  data.  Consequently,  the  disagreement  between  the  actual  and  ideal  sector  boundaries 
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Figure  18.  Incompressible  Data  Showing  versus  :  in  Sector  1 
for  tr  =  84  Degrees 

of  Figure  17  does  not  significantly  Influence  the  determination  of  roll  angle. 

Despite  the  small  error  in  calculating  the  roll  angle,  the  standard  deviation  of 
the  error  in  calculating  pitch  angle  at  4.27  degrees  is  unacceptably  high,  especially 
when  compared  to  the  0.84-degree  standard  error  of  previous  incompressible  calibrations. 

The  actual  variation  between  the  experimental  data  and  its  polynomial  prediction  is  pre¬ 
sented  in  Figure  19,  which  depicts  a  typical  outer  sector.  From  this  illustration,  we  can 
see  that  the  greatest  error  in  determining  the  pitch  angle  occurs  at  high  angles  of  attack. 
This  occurs  as  a  result  of  the  polynomial's  inability  to  fit  the  actual  data.  Typically, 
the  shape  of  the  Cn  versus  a  curve  looks  much  like  the  lift  curve  of  a  stalling  airfoil. 

Figure  20  illustrates  such  a  curve  for  sector  one,  based  on  the  extensive  data  available 

from  incompressible  probe  calibrations.  The  reason  why  the  curve  hooks  over  as  it  does 
is  evident  after  examining  the  two  pressures  comprising  the  numerator  of  the  coefficient 

of  pitch.  According  to  Figure  20,  the  center  port  pressure  decreases  with  increasing 

angle  of  pitch.  Beyond  some  point,  a  suction  develops  at  this  port,  causing  the  pressure 
there  to  dip  below  the  free  stream  static  pressure.  But  near  80  degrees  angle  of  pitch, 
the  suction  breaks  and  the  pressure  begins  to  increase.  As  this  occurs,  the  slope  of 
the  Cp7  curve  approaches  the  slope  of  the  Cp,  curve.  Once  the  two  slopes  are  equal, 
the  rate  of  change  of  the  numerator  is  zero,  causing  the  slope  of  the  coefficient  of 
pitch  also  to  be  zero.  No  calibration  may  be  made  beyond  this  point,  because  each  value 
of  Cq  then  corresponds  to  two  values  of  9.  And  since  9  is  a  function  of  Cp ,  the  poly- 
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Figure  19.  Three-Dimensional  Data  Set  Showing  Correlation  Between 
Experimental  and  Polynomial  Data  for  Sector  5 
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Figure  20.  Mechanism  for  the  Breakdown  of  Linearity  in  the  Coefficient 
of  Pitch  at  High  Angles  of  Attack 
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nomial  can  only  calculate  one  unique  value  of  0  for  a  given  value  of  Cg.  Therefore, 
based  on  the  Information  presented  in  Figure  20,  we  would  expect  to  be  capable  of  cali¬ 
brating  the  probe  out  to  80  degrees  angle  of  pitch.  This  was  indeed  the  case  for  the 
past  incompressible  calibrations,  but  not  the  case  in  our  compressible  calibration,  as 
Figure  19  vividly  points  out. 

The  failure  of  the  calibration  at  high  angles  of  pitch  is  most  likely  the  result 
of  probe  geometry.  Specifically,  the  half  angle  of  the  probe  used  in  our  compressible 
calibration  is  25  degrees  (see  Figure  10)  as  opposed  to  the  30-degree  half  angle  of  the 
probes  used  in  the  incompressible  calibrations.  This  steeper  half  angle  was  incorpor¬ 
ated  to  permit  a  closer  approximation  of  total  pressure  (as  measured  by  the  peripheral 
ports)  at  high  angles,  in  the  hope  of  extending  the  range  of  calibration  beyond  80  degrees 
of  pitch.  However,  this  reasoning  overlooked  the  effect  increasing  the  half  angle  would 
have  on  the  central  port.  That  is,  a  steeper  half  angle  requires  the  flow  passing  around 
the  probe  tip  to  turn  through  a  greater  angle,  causing  the  suction  over  the  central  port 
to  break  at  a  lower  angle  of  pitch.  Consequently,  the  slope  of  the  coefficient  of  pitch 
levels  off  earlier,  limiting  the  calibration  to  a  pitch  angle  below  80  degrees.  But 
even  as  this  theoretical  upper  limit  of  the  calibration  is  approached,  the  breakdown 
in  linearity  of  the  coefficient  of  pitch  with  increasing  angle  of  pitch  forces  the  poly¬ 
nomial  to  work  harder  to  fit  the  data.  Since  at  high  angles  of  attack  small  changes  in 
the  coefficient  of  pitch  result  in  large  changes  in  the  angle  of  pitch,  a  small  mismatch 
between  the  actual  data  and  the  polynomial  curve  fit  translates  into  a  large  error  as 
Figure  21  illustrates.  Thus,  as  a  result  of  the  difficulty  in  fitting  the  actual  data 
near  the  theoretical  limit  of  calibration  as  established  by  the  probe’s  geometry,  the 
calibration  appears  to  be  limited  to  an  angle  of  pitch  below  80  degrees. 

To  confirm  these  suspicions,  we  ran  a  second  data  reduction,  excluding  all  80-degree 
pitch  data.  The  standard  deviations  of  the  truncated  data  set  are  displayed  in  Table  3 
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Figure  21.  Limitations  of  Third  Order  Curve  Fit 
at  High  Angles  of  Attack 
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Table  3 

STANDARD  DEVIATIONS 
COMPRESSIBLE  FLOW  CALIBRATION 


INNER  SECTOR 

— 

OUTER  SECTORS 
(Q  «  80°  Data  Truncated) 

S  ECTOR^-- 

^ — ''expression 

7 

SECTOR  - 

EXPRESSION 

AVERAGE 
!  -  6 

0.78° 

0(4) 

0.79° 

o(bT) 

0.72° 

0(4,) 

0.39° 

°<PoL> 

2.5% 

°<PoL> 

1.1% 

PoL  -  P~L 

PoL  '  P«L 

U(PoL  -  P»L> 

1.4% 

°(PoL  '  P»L> 

4.1% 

P  -  P 

oL  ~L 

P  ~  P 

oL  *>L 

O(M) 

0.006 

o(M) 

0.022 

and  reflect  a  significant  decrease  in  the  error  associated  in  calculating  the  experimen¬ 
tal  quantities.  For  example,  the  standard  deviation  in  calculating  pitch  angle  went 
down  from  4.27  degrees  to  0.79  degrees.  The  correlation  between  the  data  points  of  sec¬ 
tor  five  also  improved  dramatically  and  is  depicted  in  Figure  22.  However,  one  must 
keep  in  mind  that  since  the  80-degree  points  were  removed  from  the  data  set,  the  data 
set  is  no  longer  as  accurately  represented.  That  is,  in  addition  to  losing  the  80-degree 
data,  we  also  lost  an  equal  number  of  associated  roll  angle  and  Mach  number  data.  Never¬ 
theless,  the  greatly  reduced  standard  deviations  support  our  contention  that  for  our 
probe  in  the  vicinity  of  80  degrees,  the  third  order  curve  fit  is  incapable  of  accurately 
fitting  the  data. 

VI.  Conclusions  and  Recommendations 

The  polynomial  expansion  in  three  variables  accurately  extends  the  calibration  of 
seven-hole  probes  intc  the  compressible  regime.  Based  on  the  reasonably  close  correla¬ 
tions  between  the  standard  deviations  of  past  incompressible  calibrations  and  the  stand¬ 
ard  deviations  of  the  compressible  calibrations,  the  method  of  Latin  Squares  furnishes 
a  sample  space  which  accurately  represents  a  much  larger  three-dimensional  parameter 
space.  Finally,  the  third  order  curve  fit  accurately  represented  the  parameter  space 
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Figure  22.  Three-Dimensional  Data  Set  Showing  Correlation  Between 
Experimental  and  Polynomial  Data  for  the  Truncated  Data  Set  of  Sector  5 
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out  to  an  angle  of  pitch  of  70  degrees,  but  fell  completely  apart  when  required  to  fit  the 
data  extending  to  80  degrees  angle  of  attack. 

It  is  recommended  that  an  additional  Latin  Square  be  run  in  an  outside  quadrant. 

A  set  of  calibration  coefficients  can  then  be  determined  using  that  Latin  Square  of  data 
only;  the  same  should  be  done  with  the  original  square.  The  results  can  then  be  examined 
to  see  how  well  the  calibration  coefficients  of  one  square  determine  the  values  of  the 
other  square  and  vice  versa.  Close  correlations  should  establish  beyond  a  reasonable 
doubt  the  validity  of  a  Latin  Square  in  representing  a  parameter  space.  Next,  both  Latin 
Squares  should  be  combined  to  determine  another  set  of  calibration  coefficients,  which 
in  turn  can  be  used  to  see  if  more  data  points  significantly  improve  the  degree  of  fit. 
Lastly,  using  the  combined  set  of  data  points,  extend  the  calibration  to  a  fourth  order 
polynomial  and  examine  it  to  see  if  the  addition  of  the  fourth  order  terms  significantly 
improve  the  degree  of  fit  up  to  80  degrees  angle  of  attack. 
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Symbols 

the  ith  value  of  a  particular  data  point,  where  A  Is  either 
a^,  Bt<  c0>  or  Cq  for  low  angles  and  -3 ,  <J>,  C0n,  or  Cq^  for 
high  angles 

the  1th  value  of  the  experimentally  known  value  of  A 

the  ith  value  of  the  polynomial  predicted  value  of  A 

coefficient  representative  of  compressibility  effects 

apparent  total  pressure  coefficient  for  low  angles 

apparent  total  pressure  coefficient  for  high  angles 

apparent  dynamic  pressure  coefficients  for  low  angles 

average  value  of  the  apparent  dynamic  pressure  coefficients 
for  a  given  sector 

apparent  dynamic  pressure  coefficients  for  high  angles 
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n=l-20 


coefficient  of  pressure  (arrived  at  by  nondimensionalizing 
an  individual  port  pressure) 

angle  of  attack  pressure  coefficient  for  low  angles 

intermediate  pressure  coefficients  used  to  determine  Ca 
and  Cg 

angle  of  sideslip  pressure  coefficient  for  low  angles 

pitch  angle  pressure  coefficient  for  high  angles 

roll  angle  pressure  coefficient  for  high  angles 

calibration  coefficient,  where  A  denotes  the  parameter  to 
which  a  particular  set  of  K's  belong,  and  n  denotes  the 
particular  term  of  the  calibration  coefficient  in  the  poly¬ 
nomial  power  series  expansion 

Mach  number  of  free  stream 

local  Mach  number 

pressure  at  port  "n" 

total  pressure  of  free  stream 

local  total  pressure 

static  pressure  of  free  stream 

local  static  pressure 

average  of  pressures  1  through  6 

local  velocity  components  with  respect  to  probe 

local  velocity  with  respect  to  probe 

angle  of  attack 

angle  between  probe  axis  and  velocity  vector  projected  on 
vertical  plane  through  probe's  x-axis 

angle  of  sideslip 

angle  between  probe  axis  and  velocity  vector  projected  on 
horizontal  plane  through  probe's  x-axis 

total  angle  between  velocity  vector  and  probe's  x-axis 

ratio  of  specific  heats 

angle  between  the  plane  containing  the  velocity  vector  and 
the  probe's  x-axis  and  the  x-z  plane  measured  positive  clock¬ 
wise  from  port  number  four  as  viewed  from  the  front 


standard  deviation 
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SEVEN-HOLE  PROBE  DATA  ACQUISITION  SYSTEM 


G.  Sisson*,  T.  Bolick**,  F.  Jonas** 

Abstract 

This  paper  describes  the  computer  software  (programs)  used  in 
conjunction  with  the  seven-hole  probe.  The  complete  data  acquisition  system 
includes  the  seven-hole  probe,  a  computer-driven  traverse  mechanism, 
computerized  data  acquisition  equipment,  and  computer  graphic  displays. 

I .  Introduction 

The  USAF  Academy  and  NASA-Ames  Research  Center  have  conducted  a  joint 
research  effort  over  the  past  three  years  which  has  resulted  in  a  new  means 
of  surveying  aerodynamic  flow  fields.  A  newly  invented  seven-hole  probe  was 
used  to  measure  pressures  and  velocities  in  the  wake  of  wind  tunnel  models 
as  part  of  this  effort  (Ref.  1-3).  With  this  flow  field  information  the 
aircraft  designer  can  improve  aircraft  designs,  increasing  range 
performance,  handling  characteristics,  and  overall  mission  performance. 

Flow  field  or  wake  surveys  are  not  new.  They  have  been  performed  for 
years  with  different  types  of  instrumentation  such  as  hot-wire  anemometers, 
pitot  tubes,  four-  and  five-hole  probes,  and  laser  doppler  velocimeters  and 
anemometers.  Each  type  of  instrumentation  has  inherent  advantages  and 
disadvantages.  The  seven-hole  probe  is  a  new  form  of  instrumentation  that 
capitalizes  on  high  data  rates,  ease  of  use,  low  costs,  and  an  excellent 
range  of  applicability. 

The  first  seven-hole  probe  surveys  (Ref.  1)  resulted  in  the 
measurement  of  local  total,  static,  and  dynamic  pressures  to  within  5 
percent  accuracy  (95  percent  confidence  level)  and  local  flow  angles  to 
within  +2  degrees  (95  percent  confidence  level)  for  incompressible  flows. 

*Hewlett-Packard ,  Ft.  Collins,  Colorado 
**Captain,  USAF,  Department  of  Aeronautics,  DFAN 


l 


‘J.'M'IOO  1'l.ltmi  t  C'Y  l*ltn  n.l.-vwwvtr*  > 


These  results  were  very  impressive.  Nevertheless,  certain  improvements  were 
possible . 

The  extension  of  the  probe  calibration  to  compressible  subsonic  flows 
(Ref.  3)  was  one  such  improvement.  This  is  very  significant  for  two 
reasons.  One,  most  aircraft  spend  the  majority  of  their  flight  time  in  this 
regime;  and  two,  even  in  an  incompressible  flow  field,  wakes  can  be 
produced  which  exhibit  compressibility. 

The  computer  programs  used  in  the  survey  technique  (Ref.  2)  placed 
many  limitations  on  the  user;  thus  another  area  for  improvement  was 
identified.  The  final  data  was  not  available  while  the  wind  tunnel  was 
still  running,  making  data  validation  difficult.  The  data  was  always 
collected  at  uniformly  spaced  points  in  the  flow,  which  resulted  in 
abnormally  large  run  times  or  unusually  sparce  data  sets.  Thus,  real  time 
or  on-line  data  analysis  was  not  possible,  resulting  in  longer  test 
programs  to  identify  pertinent  flow  field  phenomena.  In  addition,  the 
computer  inputs  were  very  cryptic  and  the  programs  were  not  reusable  due  to 
poor  programming  techniques  and  lack  of  program  documentation. 

This  report  discusses  improvements  to  these  computer  programs  that 
alleviate  these  limitations.  The  complete  program  docua entation  follows  in 
the  appendices  for  the  seven-hole  probe  data  acquisition  system. 

II .  Discussion 

Data  validation  was  a  problem  with  the  old  software-  Four  computer 
programs  were  executed  in  sequence  to  generate  a  plot  which  would  tell  the 
user  if  all  the  electrical  and  pneumatic  connections  for  the  experiment 
were  correct.  A  schematic  of  this  procedure  is  shown  in  Figure  1.  The 
second  program,  TOPWING,  was  the  only  one  that  used  the  wind  tunnel 


PROGRAMS 


DATA  FILES 


CCEF7 


SCALAR 


Figure  1.  Schematic  of  Old  Software 

apparatus.  Because  wind  tunnel  schedules  are  always  tight,  it  was  very 
tempting  to  run  this  program  over  and  over  until  all  necessary  wind  tunnel 
data  was  collected.  Only  then  could  the  data  be  reduced,  plotted,  and 
analyzed.  On  more  than  one  occasion  the  data  proved  to  be  erroneous  or 
invalid  because  of  changes  in  the  electr-ical  and  pneumatic  connections 
during  the  test. 

The  solution  to  this  problem  was  to  make  a  single  computer  program 
perform  the  same  functions  of  the  previous  four  programs.  This  single 
program  could  then  be  made  to  include  data  reduction  and  display  while  the 
test  was  being  conducted,  thus  allowing  on-line  data  analysis.  Any  data 
anomalies  could  Chen  be  detected  immediately,  allowing  the  user  either  to 
retake  the  data  or  inspect  the  instrumentation.  However,  a  single  program 
proved  to  be  impossible,  as  it  required  too  much  space  in  the  computer  (RAM 
or  core  memory). 
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The  new  software  consists  of  three  programs  that  use  a  common  data 
file  format.  These  three  programs  are  linked  together  through  system  level 
control  language  command  files,  thus  abbreviating  core  storage 
requirements.  The  result  is  that  the  user  sees  only  one  task  instead  of 


® - "I  START  CMPSHP 


Oafina  Vofuma 
To  8a  Sampiad 


Mar*  Any  Point* 
That  Wflf  8* 
Ovariayad  A* 
Inviaujla 


[Take  100  Points  of 
I  Data  m  Oastred 
VoHjma  &  Add  It 
j  To  Pll#— Out  File 
is  OATAOUT 


6NO  CMPSHP 


Daiata  Old  OATAIN 
&  Malta  OATAOUT 
Naw  OATAIN 


f 


© 


START  MENU 


Get/Chenge  Menu 
Deecriptton  lor 
Plotting  Program 


ENO  MENU 
t  START  GRAFS 


Draw  tha  Graph 
as  Reouastad 


END  GRAFS 


* 


Claan  Up 
Tamporary  Usa 
Filaa 


Figure  2. 


Compressible  Seven-Hole  Prooe  DAS  Schematic 


4 


h«  ^iouuc«*<j  on  *  nont  copier 


several  independent  tasks.  This  is  illustrated  in  Figures  2  and  3,  which 
show  the  flow  for  data  acquisition  tasks  and  for  the  separate  plotting 
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Figure  3.  Plots  Program  Schematic 


5 


tvc-pi uouceu  uiia  ua  uovermiiem  copier, 


tasks  respectively.  The  complete  program  listings  are  in  Appendix  A.  The 
first  program,  CMPSHP,  performs  the  functions  previously  done  by  VOLCAL, 
TOPWING,  and  C0EF7  .  It  allows  the  user  to  specify  the  volume  of  space  to  be 
sampled  using  a  three-degree-of -freedom  traversing  mechanism  (Ref.  2).  The 
data  is  then  acquired  in  the  defined  space  and  reduced  to  pressure 
coefficients,  angles,  and  velocity  ratios.  The  second  program,  MENU,  is  a 
highly  interactive  program,  which  allows  the  user  to  define  the  desired 
graphic  display  of  the  data.  The  third  program,  GRAFS,  uses  the  description 
provided  by  MENU  and  produces  the  desired  graphic  display  of  the  data.  The 
data  can  be  displayed  as  contour  and/or  axonometric  plots  of  the  scalar 
data,  and/or  vector  plots  of  the  cross-flow  velocity  data.  These  last  two 
programs,  GRAFS  and  MENU,  perform  the  same  logial  function  as  the  old 
SCALAR  program,  but  with  an  enhanced  user  interface  and  greater 
flexibility. 

Using  the  old  software,  the  seven-hole  probe  was  positioned  at  points 
in  a  rectangular  grid  which  had  nearly  uniform  spacing  between  points,  as 
shown  in  Figure  4.  This  was  the  most  efficient  spacing  without  a  priori 
knowledge  of  the  location  of  the  pressure  and  velocity  gradients  in  tne 
flow.  Nevertheless,  it  meant  collecting  a  large  amount  of  data  in  small 
grids  for  both  the  uniform  flow  and  the  gradient  regions.  This  condition 


Figure  4.  Uniformly  Spaced  Grid 
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was  far  from  optimum  from  a  data  collection  point  of  view.  What  was 
required  was  a  non-unif orraly  spaced  grid  which  concentrated  points  in 
regions  of  interest,  as  shown  in  Figure  5.  Thus,  the  software  was  modified 
to  display  this  non-unif ormly  spaced  data  as  a  series  of  grids.  In 
addition,  these  grids  were  allowed  to  overlap  so  that  two  or  more  might 
cover  the  same  position  in  the  plots.  This  allowed  more  flexibility  in 
collecting  the  data  and  still  did  not  require  an  a  priori 

knowledge  of  the  location  of  major  flow  field  features.  Data  acquisition  is 
thus  accomplished  by  first  examining  a  widely-spaced  grid.  After  examining 
the  data,  a  newer,  denser  grid  is  then  specified  by  the  user  to  cover  only 
those  regions  where  more  detail  is  required.  This  procedure  can  be  repeated 
as  many  times  as  is  required  to  get  the  necessary  resolution.  This  is  what 


Figure  5.  Non-Uniformly  Spaced  Grid 

is  meant  by  interactive  data  collection. 

Although  this  technique  solved  the  problem  of  non-unif ormly  spaced 
data,  it  caused  another:  which  data  should  be  used  at  a  point  where  more 
than  one  grid  is  defined  there?  The  obvious  and  perhaps  best  solution  woulu 
be  to  average  the  data.  The  new  software,  however,  does  not  perform  any 
averaging  in  the  interest  of  saving  time  and  space.  Instead,  the  points  in 
the  previously  defined  grids  are  marked  as  invisible  by  the  newly  defined 
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grid,  as  shown  in  Figure  6,  and  illustrated  in  Appendix  B. 

With  the  previously  discussed  modifications,  the  new  software  provided 
for  interactive  data  collection  and  plots  of  the  data  available  within 
seconds  after  the  data  collection  was  completed.  It  was  a  simple, 
inexpensive,  and  highly  accurate  means  of  flow  visualization,  which 
provided  quantitative  results  for  many  aerodynamic  flow  proper  > .  The 
software  still  suffered  one  major  flaw,  however,  in  that  it  was  not 


Figure  6.  Multiple  Grid  Display  with  Invisible  Points 

user-or iented .  The  keyboard  input  section  of  the  display  program  was  a 
series  of  questions.  Running  this  program  frequently  required  answering 
these  same  questions  many  times.  This  created  a  situation  for  possible 
errors  and  turned  out  to  be  a  source  of  frustration  for  the  user.  Improving 
the  input  section  of  the  display  program  was  solved  by  moving  all  t he 
question  and  answer  software  to  the  MENU  program.  The  user  responses  were 
then  written  out  to  a  scratch  file,  which  was  read  in  when  the  MENU  program 
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was  next  executed.  This  data  is  then  displayed  on  the  computer  terminal 
screen,  and  the  user  is  free  to  change  any  parameter  in  any  sequence.  When 
the  parameters  displayed  are  what  the  user  wants,  a  null  input  (no  change 
desired)  signals  the  program  to  plot  the  data.  If  the  plot  required  is  to 
be  the  same  as  the  last  one,  then  only  a  null  input  is  required.  This  is 
much  faster  and  less  error-prone  than  the  old  approach.  In  addition,  this 
allows  the  user  to  check  his  input  parameters  constantly.  Thus,  moving  all 
the  user  input  software  to  the  MENU  program  has  the  adued  benefit  of  making 
the  remaining  display  software  a  fixed  reuseable  unit,  a  sort  of  macro 
subroutine.  This  user  interface  is  shown  in  Appendix  B,  exactly  as 
presented  to  the  user. 

Finally,  the  software  is  sufficiently  general  enough  to  serve  as  a 
framework  for  use  with  other  flow  field  measurement  devices  such  as 
hot-wire  anemometers  or  the  four-  and  five-hole  probes.  Very  little 
modification  would  be  required  to  the  existing  software  for  these  devices. 

In  addition,  a  standard  data  file  format,  as  shown  in  Appendix  C,  has  been 
defined  to  standardize  the  acquisition  and  recording  of  experimental  data 
for  use  in  these  programs. 

III.  Conclusions 

In  conclusion,  the  highly  successful  seven-hole  probe  flow  survey  tech¬ 
nique  has  been  improved  through  software  modifications  to  allow: 

a.  interactive  data  collection  in  regions  of  interest 

b.  user-directed  specified  non-uniformly  spaced  data  grids 

c.  the  time  required  to  collect  and  analyze  data  to  be  reduced 

In  essense,  the  process  has  become  a  computer-aided  flow  visualization  system, 
user-oriented  enough  to  allow  one  to  perform  wake  surveys  without  expert 
intervention. 
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APPENDIX  A 


This  appendix  contains  complete  and  documented  listings  for  the  fol¬ 
lowing  programs  used  in  the  seven-hole  probe  data  acquistion  system: 


a.  CMPSHP 

b.  MENU 

c.  GSUBS 

d.  COMCAL 

e.  FILEIO 

f.  GRAFS 
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03  03  O)  to 

CO  CO  03  0) 

U3  03  03  03 

03  03  03  03 

03  03  03  CO 

03  03  03  CO 

03  03  03  03 
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ft  ft 
ft  ft 
ft  ft 
ft  ft 
ft  ft 
ft  ft 

ftftfta.ftftftftfta.fta. 
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X  x  xx  : 
X  X  x  X  : 


XZXZXXXSX 

zxzxxxxxz 
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trrrsr tr ttrr rr 
zscssctssscscc 
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liiiiiillSi 


UUUUUUUUUu 
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C 

c 

C  COMPRESSIBLE  FLOW  SEVEN  HOLE  PROBE  DATA  ACQUISITION  PROGRAM 

C 

c 

CMHHtHHHHIHHtHHHtMHHHHHHftMHHtMHtHHtHtHH 

C 

C  THIS  PROGRAM  IS  DESIGNED  TO  PERFORM  DATA  ACQUISITION  AND 

C  REDUCTION  OF  DATA  TO  BE  TAKEN  WITH  A  SEVEN  HOLE  PROBE. 

C 

C  THE  INPUTS  AND  OUTPUTS  OF  THE  PROGRAM  ARE  AS  FOLLOWS: 

C 

C  FILE  ’DATAIN.DAT*  IS  THE  INPUT  DATA  FILE.  IT  CONTAINS  ANY 

C  DATA  THAT  MAY  HAVE  ALREADY  BEEN  TAKEN  IN  THE  AXIAL  LOCATION 

C  PLANE  TO  WHICH  WE  WISH  TO  ADD  MORE  DATA.  THE  DATA  IS  STORED  IN 

C  THE  STANDARD  DFAN  DATA  FILE  FORMAT.  EACH  RECORD  HAS  20 

C  DIMENSIONS  AS  DEFINED  BELOW.  THERE  ARE  100  DATA  POINTS  TO 

C  EACH  PLANE  OF  DATA  CONSISTING  OF  10  LINES  OF  10  POINTS  EACH. 

C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r, 

c 

c: 


FILE  ’DATAOUT.DAT’  IS  THE  OUTPUT  DATA  FILE.  IT  UILl  CONTAIN 
EVERYTHING  THAT  WAS  IN  ’DATAIN.DAT'  AND  THE  DATA  THAT  WILL 
BE  TAKEN  DURING  THIS  RUN.  HOWEVER  t  ALL  DATA  POINTS  THAT  WERE 
IN  ’DATAIN.DAT*  THAT  WILL  BE  OBSCURED  (COVERED  UP/OVERLAYED ) 

BY  LOCATIONS  IN  THE  DATA  TAKEN  DURING  THIS  RUN  WILL  NOW  BE 
MARKED  AS  INVISIBLE.  THE  DATA  IS  STORED  IN  THE  SAME  FORMAT  AS 
FOR  ’DATAIN.DAT’. 

FILE  ’TARE.DAT*  IS  BOTH  AN  INPUT  AND  AN  OUTPUT  FILE*  DEPENDING 
ON  THE  RESPONSE  OF  THE  USER.  IT  WILL  INITIALLY  BE  A  NULL  FILE 
INDICATING  THAT  TARE  DATA  (AMPLIFIER  ZERO  CONDITIONS)  MUST  BE 
TAKEN.  THIS  DATA  WILL  BE  OUTPUT  TO  A  NON-NULL  ’TARE.DAT’  FILE 
WHICH  UILL  BE  USED  THE  NEXT  TIME  AROUND.  OLD  TARE  DATA  MAY 
BE  IGNORED  AND  NEW  TARE  DATA  TAKEN  AT  THE  REQUEST  OF  THE  USER. 
FILE  ’TARE.DAT*  HAS  ONLY  ONE  RECORD  WITH  12  ENTRIES  r  ONE  FOR 
EACH  INPUT/OUTPUT  DATA  VALUE. 


C  FILE  ’CPROBE.DAT’  IS  THE  DATA  FILE  THAT  CONTAINS  THE  CALIBRATION 

C  COEFFICIENTS  FOR  THE  C0MFRE3S I BLE  SEVEN  HOLE  PROBE.  IS  IS  AN 

C  INPUT  FILE  FOR  THIS  PROGRAM. 

C 

C  EACH  DATA  POINT  IN  ’DATAIN.DAT’  AND  ’DATAOUT.DAT’  HAVE  THE 

C  FOLLOWING  20  DIMENSIONS  IN  THE  ORDER  SHOWN: 

C 

C  ENTRY  DESCRIPTION 

C 

C  1  Y  LOCATION  -  INCHES 

C  2  Z  LOCATION  -  INCHES 

C  3  VISIBLE  ( 0  =  YES  ?  1=N0> 

C  4  ALPHA  -  DEGREES 

C  5  BETA  -  DEGREES 


6 

9 

o 

10 
11 
1 :: 
1.2 
14 


ALPHAT  -  DEGREES 


BETAT  -  DEGREES 


THETA  -  DEGREES 
PHI  -  DEGREES 
CTOTAL 
CSTATIC 


ua  » luvernmt"  *ipi<T 


c 

15 

CZERO 

c 

16 

C  D  Y  N 

c 

17 

CM 

r 

13 

RM  ACH 

2 

19 

U  VELOCITY  COMPONENT 

c 

20 

V  VELOCITY  COMPONENT 

cttnt ***********************************************************  i  ***  i 

r» 

C  DATA  DECLARATIONS  SECTION 

C 

C* *********************************************** ************* #****4 *  * 

DIMENSION  DAT ( 20  * 100) ,IDLIST(2C) , XX ( 4 ) , THAT ( 20 ) ,  Y  Y  (  4  )  , 

1  RSCRC50) • ZERO < 12) , CAL( 12)  , STRAP ( 12)  * 

2  K ALPHA < 7, 20)  , KBETA ( 7 , 20 >  r KQ ( 7 , 20 ) , KZERO ( 7 , 20 ) , KM ( 7 , 20 ) 

REAL  KALPHA»KBETA»KQ»KZERO,KM 

LOGICAL *1  ST  I TLE ( 20 ) , L3CR( 10,20),DLABEL(10,20) f LT I TLE ( SO ) 

LOGICAL  EOF 

C*** ******************  ********* **************************** **** ********* 

c 

C  DATA  PRESETS  AREA 

C 

c***** ***************************************************************** 


NAMES/LABELS  FOR  EACH  DIMENSION 
DATA  DLABEL/ ' Y ' , '  ' , '  ' , '  ' , '  ' 


'  T  '  ,  '  A  ' 
'P' » 'H' 
'  T '  r'  A ' 
'Z't'T' 
'  I '  »  '  ' 


'  T  '  »  '  A  ‘ 


'  E  '  ,  '  R  ' 
'  Y '  , 'N' 

'  '  t  '  ' 

'  C  '  ,  '  H  ' 
'  V  '  ,  '  E  ' 


LIST  OF  ALL  VARIABLES  TO  EXTRACT  FROM  INPUT  DATA  FILE 
GET  THEM  ALL 

DATA  IDLI3T/1 ,2,3, 4,5,o,7,8,9, 10, 11 » 12, 13, 14, 15, 16 , 17 , 13 , 19 , 20/ 
DATA  YES/'Y'/ 

CALIBRATION  COEFFICIENTS  FOR  THE  AMPLIFIERS  AND  A/D  CONVERTORS 
THE  ORDER  IS  AS  FOLLOWS?  HORIZONTAL  (Z)  LOCATION  ,  PI,  PC-  PS, 
P4 »  PS,  P6 ,  P  TOTAL,  P  STATIC,  P7,  VERTICAL  (Y)  LOCATION,  TOTAL 
TEMP,  DESIRED  VERTICAL  LOCATION,  DESIRED  HOR I ZONT AL  ;_0CA7I0N 


DATA  CAL/ 10 . , 1 .2872, 1 .2229,1 .2136,1 .2007 , 1 . 2756 ,1.23*’, 

1 . 1487 , 1 . 0955, 1.2218,10.,!./ 

STRAPPING  VOLTAGE  FOR  ALL  AMPLIFIERS  -  ALL  AT  +/-  1  VOL T 
DATA  STRAP/ 1 2*1 . / 

NUMBER  OF  DIMENSIONS  IN  A  DATA  POINT  IS  20 
DATA  ND/20/ 

DATA  EOF/. FALSE./ 


w.*>  tiovof  MiuPiil  copier 


C******************************************************* ******  ********* 

c 

C  FORMAT  STATEMENTS 

C 

C **********  ********  v *****************************************  ******** 

C  FORMAT  STATEMENT  FOR  FIRST  HEADER  RECORD  OF  STD  FILE  FORMAT 

11  FORMAT (515) 

12  FORMAT  (  '  DO  YOU  WANT  TO  RETAKE  TAR  E/AM  F'LIFIER  ZERO  DATA7  EY/NI') 

13  FORMAT (Al) 

14  FORMAT (  '  THE  PRESENT  DATA  TITLE  I S  t ' /2X  »  40 A  1 / 

1  '  DO  YOU  WISH  TO  CHANGE  IT7') 

15  FORMAT (  '  ENTER  DATA  TITLE  (40  CHARS  MAX)') 

16  FORMAT ( 40A1  ) 

17  FORMAT (  '  START  TUNNEL  AND  ENTER  RETURN  WHEN  READY') 

13  FORMAT  (  '  ENTER  ATMOSPHERIC  PRESSURE  (  I N - H (3  >  '  ) 

19  FORMAT ( F 1 0 . 0 ) 

20  FORMAT (  '  COORDINATE  PAIRS  ARE  ( VERT , HOR IZ >  '  > 

21  FQRMATC2F10.0) 

22  FORMAT (  '  ENTER  LOCATION  OF  LOWER  LEFT  AND  ENTER  RETURN') 

23  FORMAT (  '  ENTER  LOCATION  OF  UPPER  LEFT  AND  ENTER  RETURN'  ) 

24  FORMAT (  '  ENTER  LOCATION  OF  LOWER  RIGHT  AND  ENTER  RETURN') 

25  FORMAT (  '  ENTER  LOCATION  OF  UPPER  RIGHT  AND  ENTER  RETURN') 

26  FORMAT (  '  LOWER  LEFT  I ' , 2F 1 0 . 3 , / ,  '  UPPER  LEFT)  ' ,  2  F  1 0 . 3  »  /  > 

1  '  LOWER  RIGHT!  '  f2F10.3f/r  '  UPPER  RIGHT.'  '  »  2  F 1  0 . 3  »  /  » 

2  '  ARE  THESE  LOCATIONS  CORRECT7') 

C************* ****** *************************  ******** ***************** 

c 

C  MAIN  PROGRAM  CODE  STARTS  HERE 

C 

C ******* ************************************************************** 

C  GET  THE  CALIBRATION  COEFFICIENTS  FILE  ASSIGNED  TO  THE  PROGRAM 

CALL  ASSIGN ( 2, ' C152f23CPR0BE.DAT' ,17) 

CALL  FDBSET (  2  , 'READONLY' , , ,7) 

DEFINE  FILE  2(7,  256  >  U  » IREC ) 

C  GET  ALL  OF  THE  CALIBRATION  COEFFICIENTS 

1REC-1 

DO  1000  1=1,7 

RE  AIK  2  '  IREC  )  ( l\ALPHA<  I ,  J  )  ,  J  =  1 , 20  )  , 

1  (KBETAd,  J)  ,  J=l,-20)  , 

2  (KQ< I , J) , J=i , 20) , 

3  ( KZERO ( I ,J) ,J=1 ,20) , 

4  'KM< I , J) • J=1 ,20) 

1000  CONTINUE 

C  WE  ARE  DONE  WITH  THE  CALIBRATION  COEFFICIENTS  DATA  FILE 

C  CLOSE  AND  RELEASE  IT 

CALL  CLOSE ( 2 ) 

C  GET  THE  LPS-11  LABORATORY  PERIPHERIAL  SYSTEM  ASSIGNED  TO  UNIT  1 

CALL  ASLSLN(l) 

C  GET  THE  INPUT  DATA  FILE  ASSIGNED  AND  ALLOCATED  TO  THE  FROGRAM 

CALL  ASSIGN(2, 'DATAIN.DAT' , 10) 

CALL  FDBSET(2, 'OLD' , 'SHARE' ) 

C  GET  THE  OUTPUT  DATA  FILE  ASSIGNED  AND  ALLOCATED  TO  THE  PROGRAM 

CALL  ASSIGN<3, 'DATAOUT.DAT' , 11) 

CALL  FDBSETC3, 'NEW' , 'SHARE' ) 

C  GET  THE  T ARE/AMPLIFIER  ZERO  DATA  FILE  ASSIGNED  AND  ALLOCATED 

CALL  ASSIGN<4, 'TARE.DAT' ,8) 

CALL  FDBSEK4, 'OLD' , 'SHARE' ) 

C  READ  THE  TARE  DATA  FROM  THE  FILE  -  IF  END  OF  FILE  IS  ENCOUNTERED 

C  THEN  WE  MUST  TAKE  TARE  DATA  -  IF  NOT  THEN  TARE  DATA  IS  OPTIONAL 

C  AT  THE  REQUEST  OF  THE  USER 

READ(4,END=1100)ZER0 

A5 
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^  ■ 


2000 


GIVE 

WRITE 

READ  ( 

IF  TH 

IF  (A 

TAKE 

CONTI 

CALL 

WE  NO 

SINCE 

CALL 

ASSIG 

ACQUI 

CALL 

CALL 

WRITE 

WRITE 

CONTI 

NOW  W 

CALL 

IT  IS 

OR  IF 

IS  A 

CONTI 

ALLOW 


THE  USER  THE  OFTION  TO  RETAKE  THE  TARE  DATA 


SO  CLOSE  AND  RELEASE  IT 
TO  MAKE  A  NEW  ONE 


DATA  FILE  TO  CONTAIN  THE 


5  ? 13 ) ATARE 

E  ANSWER  WAS  NOT  YES  THEN  SKIP  TAKING  TARE  DATA  AGAIN 
TARE. NE. YES)  GO  TO  1200 
TARE/AMPLIFIER  ZERO  DATA  NOW 
NUE 

TARE' ZEROiSTRAP) 

LONGER  WANT  THE  INPUT  TARE  FILE  SO  CLOSE  AND  RELEASE  I 
IT  IS  NOW  OBSOLETE  WE  WILL  WANT  TO  MAKE  A  NEW  ONE 
CLOSE ( 4 ) 

N  AND  CREATE  A  NEW  TARE  DATA  FILE  TO  CONTAIN  THE  JUST 
RED  AMPLIFIER  ZERO  DATA 
ASSIGN (4, ' TARE . DAT ' ? 8 ) 

FDBSET ( 4  r 'NEW' ? 'SHARE'  ) 

THE  NEW  DATA  TO  DISK 
<  4 ) ZERO 
NUE 

E  ARE  THRU  WITH  THE  TARE  FILE  COMPLETELY  SO  RELEASE  IT 
CLOSE  <  4 ) 

NOW  TIME  TO  SEE  IF  THE  INPUT  DATA  FILE  WAS  A  NULL  FILE 
THIS  RUN  IS  ADDING  DATA  TO  AN  EXISTING  FILE.  IF  IT 
NEW  RUN  THEN  THE  TITLE  MUST  BE  INPUT  ?  IF  IT  IS  A 
NUATION  RUN  THEN  WE  WILL  TELL  THE  USER  THE  TITLE  AND 
HIM  TO  CHANGE  IT  IF  HE  SO  DESIRES. 


FILE  COMPLETELY  SO  RELEASE 


READ  THE  FIRST  RECORD  -  IF  END  OF  FILE  THEN  GO  TO  GET  TI^LE 
READ (  2  ?  1 1 ?END=1300)MND?NP?NL?NGf NLAT 

GET  THE  REST  OF  THE  HEADER  RECORDS  FROM  THIS  GRID?  INCLUDING 
THE  TITLES. 

CALL  GETHDRI2? MND ? ND ? I DL I ST ? LSCR ? DL ABEL ? ST ITLE ? LT I TLE ? EOF ) 
OUTPUT  THE  TITLE  TO  THE 'USER  -  SEE  IF  HE  WANTS  TO  CHANGE  IT 
WRITE(5?14) <  LTITLE ( J ) ?  J  =  1 ?  40 ) 

R E A D  <  5  ? 1 3 ) AT  I TL 

IF  THE  USER  DOESN'T  WANT  TO  CHANGE  IT  THEN  GO  TO  THEN  NEXT  STEP 

IF  (ATITL.NE.YES)  GO  TO  2000 

CONTINUE 

ASK  FOR  AND  READ  IN  THE  40  CHARACTERS  FOR  THE  TITLE 
WR  ITE ( 5  ? 15 ) 

READ ( 5  ? 1 6 ) (LTITLE(I) ? 1=1? 40) 

MAKE  THE  TIME  AND  DATA  THE  LAST  20  CHARACTERS  OF  THE  60 
CHARACTER  TITLE  ARRAY 
CALL  T IME ( LTITLE ( 41 ) ) 

CALL  DATE ( LTITLE ( 4? ) ) 

CONTINUE 

TELL  THE  USER  TO  START  THE  TUNNEL  AND  HIT  ‘RETURN*  KEY 
WHEN  HE  IS  READY  TO  CONTINUE 
WRI TE  <  5 » 1 7 ) 

CLOSE  THE  INPUT  DATA  FILE  AND  GET  IT  REASSIGNED  SO  WE  ARE 
BACK  TO  THE  START  OF  THE  FILE.  NORMALLY  THIS  WOULD  SIMPLY 
BE  A  REWIND  STATEMENT  BUT  THAT  ISN'T  WORKING  ON  SEQUENTIAL 
DISK  FILES  ON  THIS  SYSTEM  AT  THE  PRESENT  TIME.  THIS  DOES 
THE  SAME  THING. 

CALL  CLOSE ( 2 ) 

CALL  ASSIGN(2? 'DATAIN.DAT' ? 10) 

CALL  FDBSET(2? 'OLD' f 'SHARE' ) 

NOW  WE  WAIT  FOR  THE  USER  TO  SIGNAL  THAT  HE  IS  READY  TO  CONTINUE 
READ'S? 1 3 ) A  G  0 

NOW  WE  ASK  FOR  THE  ATMOSPHERIC  PRESSURE  READING  IN  INCHES  OF 
MERCURY  WHICH  WE  WILL  THEN  CONVERT  TO  PSIA  USING  THE  CONSTANT 
FOR  MERCURY  AT  THIS  ELEVATION  AND  TEMPERATURE. 


INCLUDING 


THEN 


THEN 


TITLE 


'RETURN' 


r 


C 

c 

c 

L 

O 

o 


( 


c 

c 


vV 


c 

i; 

c 

c 

c 

c 

c 

c 

c 

c 


c 


3000 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


3100 

c 

c 

c 


c 


3200 

C 

c 


WRITE! 5 . 1 3  ) 

R  E  A  D  ( 5  •  1  ?  )  P  A  T  M 
PATM  =  F'ATM*  .4392 
CONTINUE 

TELL  THE  USER  THAT  COORDINATES  ARE  TO  BE  ENTERED  AS  FAIRS 
IN  THE  ORDER  ! VER T ICAL . HORIZONTAL /  WHERE  THE  COORDINATE  IS 
A  LOCATION  IN  INCHES  FROM  THE  USER  AND  TRAVERSE  DESIGNATED 
ORIGIN . 

WF:  I TE  (  5  >  20  ) 


ASK  FOR  COORDINATES  FOR  THE  LOWER  LEFT  CORNER  OF  THE  PLANE 
URITE ( 5 . 22 ) 

READ ! 5 f  21  )  XX  •  1 ) . YY ( 1  ) 

ASK  FOR  COORDINATES  FOR  THE  UPPER  LEFT  CORNER  OF  THE  PLANE 
WRITE! 5 . 23) 

READ (5.21 ) XX(2) . YY<2) 

ASK  FOR  COORDINATES  FOR  THE  LOWER  RIGHT  CORNER  OF  THE  PLANE 
URIfE!5.24) 

READ! 5 *21 )XX!3) . YY(3) 


ASK  FOR  COORDINATES  FOR  THE  UPPER  RIGHT  CORNER  OF  THE  PLANE 
UR  I TE  <  5 . 25  ) 

RE AD  <  5  *  2 1 >  XX  C  4 ) .YY(4) 

NOW  THAT  WE  HAVE  THE  COORDINATES  ASK  THE  USER  TO  VERIFY  THAT  HE 
HAS  ENTERED  THEN  CORRECTLY. 

WRITE! 5.  26 ) < ( XX ( I >  *  YY  < I ) ) . 1  =  1 , 4 ) 

READ! 5. 13) AOK 

IF  THE  COORDINATES  ARE  NOT  RIGHT  THEN  GO  BACK  AND  TRY  AGAIN 

IF  ! AOK . NE . YES ) GO  TO  2100 

CONTINUE 

NOW  WE  PROCESS  THE  DATA  THAT  EXISTS  (IF  ANY)  ON  THE  INPUT  DATA 
FILE.  EACH  PLANE/GRID  OF  DATA  IS  PROCESSED  AND  ANY  DATA  POINTS 
THAT  ARE  COVERED  BY  THE  NEWLY  DEFINED  PLANE/GRID  WILL  BE  MARKE 
AS  INVISIBLE .  THEN  TH£~DATA  WILL  BE  OUTPUT  TO  THE  OUTPUT  DATA 
FILE  TO  WHICH  THE  NEW  DATA  WILL  BE  APPENDED  A  LITTLE  LATER. 


GET  A  GRID/PLANE  OF  DATA*  20  DIMENSIONS  PER  POINT.  100  POINTS 
CALL  GETGRD  < 2 . ND . IDLIST . 100 . LSCR . R3CR . LSCR . LSCR . DLABEL. 

DAT. NR. NL.NG.NL AT. EOF) 

IF  THERE  IS  NO  MORE  DATA  TO  PROCESS  THEN  CONTINUE  TO  THE  NEXT 
STEP  IN  THE  PROCESS  -  GO  TO  ACQUIRE  THE  NEW  FLANE  OF  DATA 
IF ( EOF ) GO  TO  3100 

DETERMINE  IF  ANY  POINTS  IN  THIS  PLANE  SHOULD  BE  MARKED  AS 
INVISIBLE  BECAUSE  THEY  ARE  CONTAINED  IN  THE  NEW  PLANE. 

THIS  ROUTINE  MARKS  POINTS  INVISIBLE  AS  REQUIRED. 

CALL  CNTAIN(NP.NL.DAT.XX.YY) 

NOW  THAT  ANY  POINTS  THAT  ARE  CONTAINED  HAVE  BEEN  MARKED  AS 
INVISIBLE.  WE  OUTPUT  THE  PLANE  OF  DATA  TO  THE  OUTPUT  FILE 
AND  SEE  IF  THERE  IS  MORE  DATA  ON  THE  INPUT  FILE  TO  PROCESS. 
CALL  F'UTL  IN  (3.  ND.NP.NL.NG.NL  AT.  ST  I  TLE.LTITLE.DLABEL.DAT) 

GO  TO  3000 
CONTINUE 

NOW  WE  HAVE  PROCESSED  ALL  THE  DATA  FROM  THE  INPUT  FILE  AND  ARE 
ALMOST  READY  TO  TAKE  THE  DATA  FOR  THIS  NEW  PLANE.  FIRST  WE 
DEFINE  OUR  PLANE  AS  10  HORIZONTAL  LINES  OF  10  POINTS  EACH. 

NP=  1 0 
NL=  1 0 

NOW  UE  MARK  ALL  POINTS  IN  THE  NEW  PLANE  AS  VISIBLE. 

DO  3200  1=1.100 
DAT ( 3  « I ) =0  . 

CONTINUE 

NOW  UE  USE  ROUTINE  BLINE  TO  DETERMINE  THE  LOCATION  OF 
EACH  DATA  POINT  IN  THE  PLANE. 


ucpi'uoiiceo  onii  us»  iioverniiie'  uier 


CALL  BLINE  < NP . NL , DAT , XX  >  YY  ) 

n  NOLI  WE  TAKE  THE  DATA  FOR  THIS  PLANE  -  THE  DATA  ACQUISITION 

C  AND  REDUCTION  CALLS  ARE  PART  OF  THE  TAKRED  ROUTINE. 

CALL  TAKRED (NP  »  NL  *  DAT  t ZERO ,  STRAP , CAL, PATH, 
r  3  KALPHA»KBETA»KQ»KZERO»KM> 

C  NOW  THAT  THE  DATA  FOR  THIS  PLANE  HAS  BEEN  TAKEN  AND  REDUCED  t 

C  WE  OUTPUT  IT  TO  THE  OUTPUT  DATA  FILE  WITH  THE  OTHER  DATA 

r  C  ALREADY  THERE  FROM  THE  INPUT  DATA  FILE. 

CALL  PUTLIN(3,ND»NP»NL»NG+1 >  N  L  A  T ,STITLE»LTITLE,D LABEL, DAT) 

C  UE  NO  LONGER  NEED  THE  INPUT  DATA  FILE  SO  RELEASE  IT. 

r  CALL  CLOSE  (  2 ) 

"  C  WE  NO  LONGER  NEED  THE  OUTPUT  DATA  FILE  50  RELEASE  IT  TOO. 

CALL  CLOSE ( 3 ) 

f  C  WE  ARE  THRU  WITH  THE  DATA  ACQUISITION  AND  REDUCTION  OF  A  PLANE 

^  C  OF  DATA  SO  TERMINATE  THE  PROGRAM. 

STOP 

(  END 

ct*****t*t*t***#****tt****ttt***t***********t******************t** 


( 

( 

C 
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L 

c 

c 
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o 
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c 

C  THIS  ROUTINE  IS  USED  TO  TAKE  THE  AMPLIFIER  ZERO  DATA 

C  FOR  ALL  THE  TRANSDUCERS. 

C 

C 

C»*t***#*****t****t**************»**##*****#******t**********tttt** 

SUBROUTINE  TARE ( ZERO , STRAP > 

C  ARRAY  ZERO  WILL  CONTAIN  THE  AMPLIFIER  ZERO  VOLTAGES 

C  ARRAY  STRAP  CONTAINS  THE  STRAPPING  VOLTAGE  FOR  EACH  CHANNEL 

DIMENSION  ZERO < 12) » STRAP (12) , DATA ( 12) 

C  CLEAR  THE  AMPLIFIER  ZERO  ARRAY 

DO  10  1=1 , 12 
ZERO ( I ) =0  . 

10  CONTINUE 

C  CALL  THE  DATA  ACQUISITION  SUBROUTINE  -  IT  ACQUIRES  THE 

C  ACTUAL  VOLTAGES  FOR  EACH  ANALOG  TO  DIGITAL  CONVERSION 

C  CHANNEL  THAT  WE  ARE  INTERESTED  IN. 

CALL  TAK'DAT  (DATA, ZERO, STRAP) 

C  THERE  ARE  CERTAIN  CHANNELS  THAT  DO  NOT  USE  AMPLIFIERS 

C  AND  FOR  WHICH  THE  AMPLIFIER' ZERO  ENTRY  MUST  BE  ZERO. 

C  THESE  ARE  AS  FOLLOWS: 

C  ACTUAL  HORIZONTAL  LOCATION 

C  ACTUAL  VERTICAL  LOCATION 

C  DESIRED  HORIZONTAL  LOCATION 

C  DESIRED  VERTICAL  LOCATION 

C  TEMPERATURE 

C  MOVE  THE  AMPLIFIER  ZERO  VALUES  THAT  WE  NEED  INTO  THE  ARRAY 

DO  20  1*2,10 
ZERO < I ) =DAT A ( I  ) 

20  CONTINUE 

RETURN 
END 

CM 1 1* **********  ***************** *#*###*#**###******#**#**** ********* 

C 

c 

C  THIS  ROUTINE  COMPUTES  THE  DESIRED  LOCATION  OF  ALL  100  DATA 

C  POINTS  BASED  ON  THE  LOCATION  OF  THE  4  CORNERS  DEFINED  IN 

C  ARRAYS  XX  AND  YY.  THE  DATA  POINTS  ARE  FROM  MINIMUM  TO 

C  MAXIMUM  HORIZONTAL  VALUE  FOR  EACH  LINE  AND  EACH  LINE  IS 

C  FOR  A  GREATER  VERTICAL  VALUE  THAN  THE  PREVIOUS  LINE 

C 

r  A8 
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C* ******** t**tt4l**tt»**t*ttt*t*********t#**********t****** ********* 

SUBROUTINE  BLI  NE  •:  N  f  M  »  DAT  f  XX  ,  YY  ) 

C  N  IS  THE  NUMBER  OF  FOINTS  PER  LINE 

C  M  IS  THE  NUMBER  OF  LINES  PER  PLANE 

C  DAT  IS  THE  DATA  ARRAY 

C  XX  IS  THE  4  HORIZONTAL  CORNER  LOCATIONS  ARRAY 

C  YY  IS  THE  4  VERTICAL  CORNER  LOCATIONS  ARRAY 

DIMENSION  DAT (20* 100)  f  XX  v  4 ) » YY ( 4 ) 

C  INITIALIZE  THE  DATA  POINT  INDEX 

K  =  1 

C  COMPUTE  THE  HORIZONTAL  INCREMENT 

UINC=1  .  /  FLOAT 'IN-l  > 

C  COMPUTE  THE  VERTICAL  INCREMENT 

VINC  =  1 . /FLOAT ( M-l  ) 

C  CLEAR  THE  VERTICAL  INCREMENT  VALUE 

V  =  Q  . 

C  CQR  EACH  LINE  IN  THE  PLANE  PROCESS  THE  LOCATIONS 

DO  20  J=  1  »  M 

C  CLEAR  THE  HORIZONTAL  INCREMENT  VALUE 

U  =  0  . 

C  FOR  EACH  POINT  ON  A  LINE  -  COMPUTE  1  HE  LOCATION 

DO  10  1=1 fN 

C  COMPUTE  THE  VERTICAL  LOCATION 

DAK  1 f K>  =  < 1-U) *( 1-V) *XX< 1 ) +( 1-U>*V*XX<2) 

1  +  U*<  l-V)*XX(3)-i-U*V¥XX<4> 

C  COMPUTE  THE  HORIZONTAL  LOCATION 

DAT<  2» K)  =  < 1-U ) *< 1-V) *YY< 1 )  +  < 1-U ) *V*YY ( 2  > 

1  +U*<  l-V)#YY(3)-i-U*V¥YY<4> 

C  INCREMENT  THE  DATA  POINT  INDEX 

K  =  K  + 1 

C  INCREMENT  THE  HORIZONTAL  INCREMENT  VALUE 

U=U+UINC 
10  CONTINUE 

C  INCREMENT  THE  VERTICAL  INCREMENT  VALUE 

V  =  V  +  V I N  C 

20  CONTINUE 

RETURN 

END 

*************************************************** ************ 


THIS  ROUTINE  DETERMINES  IF  EACH  OF  THE  100  DATA  POINTS 
IS  CONTAINED  IN  THE  BOUNDARYS  OF  THE  NEWLY  DEFINED  PLANE. 


********************  *******  ****-***********  *******  **********  **** 

SUBROUTINE  CNTA IN ( N r n f DAT , XX r Y Y > 

N  IS  THE  NUMBER  OF  POINTS  FER  LINE 
M  IS  THE  NUMBER  OF  LINES  PER  PLANE 
DAT  IS  THE  DATA  ARRAY 

XX  IS  THE  HORIZONTAL  CORNER  LOCATIONS  ARRAY 
YY  IS  THE  VERTICAL  CORNER  LOCATION  ARRAY 
DIMENSION  DAT (20? 100) f  XX ( 4 ) »  YY ( 4 ) 

LOGICAL  NAREA 

C  SET  UP  THE  DATA  POINT  INDEX 

K  =  1 

C  FOR  EACH  LINE  IN  THE  PLANE  PROCESS  THE  DATA  FOINTS 

DO  10  J  =  1 »  M 

C  FOR  EACH  POINT  IN  A  LINE  PROCESS  THE  DATA  POINT 

DO  10  1=1 fN 

C  COMPUTE  THE  SLOPE  AND  INTERCEPT  INFORMATION  FOR 

A9 
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RIGHT 


THE  PLANE 


THE  LINE 


THE  PLANE 


:  THE  LINES  THAT  MANE  UP  THE  FOUR  SIDES  OF  THE  PLANE. 

A  1  =  (  :<X  C  1  )  -XX  <  3  )  )  /  (  YY(  1  )  -YY  (  3  >  ) 

A  2  =  (  X  X  (  2  )  -  X  X  (  4  )  >/(YY(2)-YY(4)  ) 

A3= C YY ( 1 ) -YY  (  2  >  V/( XX (  1 >-XX(2 1  ) 

A  4  =  (  YY  (  3  )  -YY  i  4  '•  >  /  <  X  X  •:  3  )  -XX  (  4  )  ) 

B 1 =  XX (  1  ^ -A1*YY< 1 ) 

B2=XX(2)-A2*YY(2> 

D  3  =  Y  Y ( 1 ) -A3*XX<  1  ) 

B4=YY(3)-A4*XX<3) 

SEE  IF  THE  DATA  POINT  IS  TO  THE  LEFT  OF  THE  PLANE 
X  -  A  1  *  D  A  T  (  2  r  K  )  +  B 1 
IF  SO  THEN  IT  IS  NOT  CONTAINED 
IF  ( X.GE.DAT Cl ,K) )G0  TO  10 

SEE  IF  THE  DATA  POINT  IS  TO  THE  RIGHT  OF  THE  PLANE 

X=A2*DAT (  2  f  K  )  +  B  2 

IF  SO  THEN  IT  IS  NOT  CONTAINED 

IF  ( X  .  LE . DAT (  1  f  K ) ) GO  TO  10 

SEE  IF  THE  DATA  POINT  IS  BELOW  THE  PLANE 

X=A3*DAT( 1 fK)+B3 

IF  SO  THEN  IT  IS  NOT  CONTAINED 

IF  ( X  .  GE . BAT ( 2  f  K ) ) GO  TO  10 

SEE  IF  THE  DATA  POINT  IS  ABOVE  THE  LINE 

X=A4*DAT( 1 fK'+B4 

IF  SO  THEN  IT  IS  NOT  CONTAINED 

IF  (X.LE.DAT (2»K) >G0  TO  10 

THE  DATA  POINT  MUST  BE  INSIDE  THE  PLANE  -  MARK  AS  INVISIBLE 
DAT  <  3  »  K ) *  1  • 

10  K  =  K  + 1 

RETURN 
END 

********************************************************************* 


THIS  ROUTINE  CONTROLS  THE  ACTUAL  DATA  ACQUISITION  AND  REDUCTION 
OF  DATA  POINTS  TAKEN  UITH  THE  SEVEN  HOLE  PROBE  IN  COMPRESSIBLE 
FLOW. 


********************************************************************** 

SUBROUTINE  TAKRED ( N f M f DAT t ZERO f STRAP f CAL f P ATM f 
1  KALPHAfKBETAfKQf  KZERO  f KM ) 

N  IS  THE  NUMBER  OF  POINTS  PER  LINE 

M  IS  THE  NUMBER  OF  LINES  PER  PLANE 

DAT  IS  THE  DATA  STORAGE  ARRAY 

ZERO  IS  THE  T ARE/AMPLIFIER  ZERO  DATA  ARRAY 

STRAP  IS  THE  AMPLIFIER  STRAPPING  VOLTAGE  ARRAY 

CAL  IS  THE  CALIBRATION  CONSTANTS  ARRAY 

K ALPHA  t  KBETA  r  KQ t  KZERO  t  KM  ARE  THE  SEVEN  HOLE 

PROBE  CALIBRATION  COEFFICIENTS  FOR  COMPRESSIBLE  FLOW. 

DIMENSION  DATA  < 1 2 ) > ISB  <  2 ) 

DIMENSION  DAT < 20 r 100) f ZERO ( 12) f STRAP (12) f CALC  12) f 
1  KMC  7 f 20) f KALPHA < 7 f  20 ) f KBETA <?f 20) fKZER0<7f20>  fKQ(7f20) 

REAL  KALPHAfKBETAfKZEROfKQfKM 

THE  I  SERF  FUNCTION  COMPUTES  THE  INDEX  OF  A  DATA  POINT  SO  THAT 
THE  TRAVERSE  TRAVEL  FROM  ONE  POINT  TO  THE  NEXT  IS  MIMINIZED • 
ISERF'(MfN)=M0D<M+1f2>*<  10*M  +  l-N)+M0D(MF2)*(10*M-l0fN) 

OUTPUT  THE  LOCATION  OF  THE  FIRST  DATA  POINT  FOR  THE  USER 
<:  THIS  WAY  HE  CAN  CHECK  THAT  ALL  IS  WELL  BEFORE  TURNING 

C  THE  TRAVERSE  MECHANISM  FROM  MANUAL  TO  AUTO 

URITE<3f807)DAT(1f1) f DAT (2-  1) 
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807  FORMAT ('  FIRST  DATA  POINT:  VERTICAL^' »  F  1  0 . 3  » 

1  '  HORIIONTAL='  ,  F  1  0 . 3  ) 

C  FOR  EACH  LIME  OF  DATA  IN  THE  PLANE  -  PROCESS  THE  POINTS 

DC  20  , 1 -  1  iM 

C  FOR  EACH  POINT  IN  A  LINE  -  PROCESS  THE  POINT 

00  20  1=1 »N 

C  GET  INDEX  OF  THE  NEXT  DATA  POINT 

K  =  I S  E  P  P  (  J  •  i  ; 

C  NONE  THE  TRAVERSE  MECHANISM  TO  THE  DESIRED  LOCATION 

CALL  MOVE < DAT' 2, K) ,  D  A  T  (  1  »  K  > * 

C  WAIT  1  SECOND  FOR  CONDITIONS  TO  STABILIZE 

CALL  WAIT ( 1 , 2, ISP ) 

C  CALL  THE  DATA  ACQUISITION  SUBROUTINE  TO  ACTUALLY  READ  THE  DATA 

CALL  TAKDAT  <  DATA  r ZERO .STRAP  ) 

DO  10  11=1.12 

10  DATA' II)=DATA( II)¥CAL(II) 

C  REPLACE  THE  DESIRED  LOCATION  WITH  THE  ACTUAL  LOCATION 

DAT  (  1  ,K  )  =  D A T A  '  .1  > 

DAT < 2 . K ) =  D  A  T  A ( 1 1 ) 

C  REDUCE  THE  DATA  TAKEN  WITH  THE  COMPRESSIBLE  SEVEN  HOLE  PROBE 

CALL  COMC ALCO.DAT A. KALPHA.KBETA. KG, KZEROf KM fPATM. DAT<4, K) , 

1  DAT (  5  ,  K  ) .DAT (6.K) .DAT (  7  ,  K  ) »DAT (8.K)  .DAT (  9  ,  K  ) .DAT ( 10, K)  , 

2  DAT(ll.K) .DAT (12, K) . D  A  T ( 1 3 , K ) . D  A  T ( 1 4 , K ) ,DAT(15,K) ,  D  A  T  ( 1 6  ,  K  ) , 

3  DAT ( 17 ,K ) , DAT ( IS  »K> , DAT ( 19.K) .DAT (20  ,K)  ) 

20  CONTINUE 

RETURN 

END 

C ********* ******  **************************************** *************** 

c 

c 

C  THIS  ROUTINE  DOES  THE  ACTUAL  DATA  ACQUISITION  THRU  THE  LPS-11 

C  ANALOG  TO  DIC1 TAL  CONVERTORS. 

C 

C 

C**t**t ************ *t* I  *  ********************************************* 

SUBROUTINE  TAKDAT ( X , ZCF , STRAP > 

DIMENSION  X ( 1 2 ) , STRAP ( 12) t ZCF (12) 

C  FOR  EACH  OF  THE  12  CHANNELS  OF  DATA  TO  BE  ACQUIRED 

DO  2000  11=1,12 
1  =  11 

C  IF  THIS  IS  THE  ELEVENTH  CHANNEL  THEN  ACTUALLY  USE  CHANNEL  0 

IF  (II. EQ. 11)1=0 

C  SET  THE  DATA  TO  ZERO  INITIALLY 

RDAT  A  =  0 . 

C  TAKE  10  SAMPLES  AND  AVERAGE  THEN 

DO  1900  J=l,10 

C  GET  THE  VALUE  OF  THE  VOLTAGE  IN  DIGITAL  FORM 

1500  CONTINUE 

CALL  ADC (I»  DATA ) 

C  CONVERT  IT  TO  A  VOLTAGE 

DATOUT=  < ( ( (DATA/64 . ) /2047 . 5) * STRAP ( 1 1 ) ) -STRAP ( 1 1  )  ) 

C  CHECK  FOR  DATA  OUT  OF  RANGE 

DATA=DATA/64. 

IF  ( (DATA. GE. 4095. ) .OR. (DATA. LE.O. ) )G0  TO  1700 
C  ADVISE  THAT  THE  CHANNEL  VOLTAGE  VALUE  WAS  OUT  OF  RANGE 

WRITE(5» 1600)  I 
GO  TO  1500 

1600  FORMAT ( IX TAKDAT !  CHANNEL  '.02,'  OUT  OF  RANGE') 

1700  CONTINUE 

RDATA=RDATA+DATOUT 
1900  CONTINUE 

All 
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2000 


C  *  *  *  *  *  * 
C 

c 

c 

c 

c 

c 

c****** 


DAT0UT=RDATA/10. 

X ( 1 1 )=DAT0UT-ZCF( II > 

CONTINUE 

RETURN 

END 

************************************************************* 
THIS  SUBROUTINE  EXPECTS  THE  HORIZONTAL  TRAVERSE  POSITION 
ON  A/D  INPUT  CHANNEL  0  AND  THE  VERTICAL  TRAVERSE  POSITION 
ON  A/D  INPUT  CHANNEL  1.  THE  *H‘  OUTPUT  CABLE  CONNECTS  TO 
THE  HORIZONTAL  DRIVE  AND  THE  *V  OUTPUT  CABLE  CONNECTS  TO 
THE  VERTICAL  DRIVE.  THIS  CHANGE  AND  COMHENTS  ADDED  BY 
CAPTS  BOLICK  AND  SISSON  ON  27  MAY  30. 
************************************************************* 
SUBROUTINE  MOVE(Z,Y) 

DIMENSION  V<2> »r (3) ,  DEL<2) , IBUF<3> , IRATE  I  2) 

DIMENSION  I3BC2) ,FLC2> ,CA(2) ,CB(2) »VC(2) ,AV<2> 

SET  THE  LIMIT  FOR  USING  MAXIMUM  DRIVING  VOLTAGE  TO  3  INCHES 
FOR  BOTH  HORIZONTAL  AND  VERTICAL 
CA 1 1 >  =3 . 

0  A  <  2  )  =  3  . 

SET  FRACTIONS  FOR  COMPUTING  DRIVING  VOLTAGES  INSIDE  OF  3  INCHES 
CB<1)=9./15. 

CB  <  2  >  =9 . / 15  . 

VCC1)=13./15. 

VC(2)=13./15. 

SET  DESIRED  HORIZONTAL  POSITION 
P  ( 1  >  =  Z 

SET  DESIRED  VERTICAL  POSITION 
P  (  2  )  =  Y 

SET  THE  COMPLETE  FLAG  FOR  BOTH  TO  ZERO 
DO  5  1  =  1,2 
FL ( I ) =0 . 

CONTINUE 

CHECK  LOCATION  AND  SET  DRIVE  VOLTAGE  FOR  BOTH 
DO  35  1=1*2 

THE  ANALOG  CHANNELS  ARE  0  AND  1  FGR  HORIZONTAL  AND  VERTICAL 
AXIS  DRIVES  RESPECTIVELY 
1  =  1-1 

SET  CURRENT  POSITION  TO  ZERO 
V(I)=0. 

TAKE  A  4  TIME  AVERAGE  SAMPLE  OF  THE  LOCATION 
DO  16  I A= 1 f  4 

CALL  ADC (  1 1 ,  A  V  (  I  ) ,1,ISB> 

<I5=V(I)+AV(I)/131072. 

V( I ) =V ( I ) / 4 . 

IF  AN  ERROR  WAS  ENCOUNTERED  ADVISE  THE  USER 
IF ( ISB  < 1 > . NE . 1 >  WRITE ( 5, 908 ) ISB ( 1 > 

COMPUTE  THE  DELTA  DISTANCE  BETWEEN  ACTUAL  AND  DESIRED  LOCATIONS 
DEL ( I ) =- ( (V(t)/l.-l.)*10.-P(I)) 

IF  GREATER  THAN  +3  INCHES  AWAY  -  USE  MAX  POSITIVE  VOLTAGE 
IF(DELCI) 'GT»CA(I) )G0  TO  20 

IF  BETWEEN  0  AND  +  3  USE  PROPORTIONAL  VOLTAGE  BASED  ON  DELTA 
IF ( DEL ( I ) . GT . 0 . ) GO  TO  25 

IF  BETWEEN  0  AND  -  3  USE  PROPORTIONAL  VOLTAGE  BASED  ON  DELTA 
IF ( DEL ( I ) .GT.-CA(I))GO  TO  30 

MUST  BE  GREATER  THAN  -  3  INCHES  AWAY,  30  USE  MAX  NEGATIVE  DRIVE 
IBUF< I+4)=INT(2047.5-VC< I) *2047. 5) 

GO  TO  35 

IBUF ( 1+4 ) = I NT (2047 . 5  + VC  < I ) *2047 .5' 

GO  TO  35 

IBUF < 1  +  4) =INT( < 1 . +CB ( I )  +  <  VC ( I ) -CB ( I ) ) *DEL ( I > /CA ( I ) >*204  7.5) 
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FOSITIVE  VOLTAGE 


IF  BETWEEN  0 
IF ( DEL ( I ) .GT, 
IF  BETWEEN  0 
IF ( DEL  < I ) .GT, 


VOLTAGE  BASED 


VOLTAGE  BASED 


DELTA 


DELTA 
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r 


V 

( 

0 

0 


J 


L 


GO  TO  33 

30  I BUF  < 1+4  >  =  INT<  C 1 .-CB< I >  +  ( VC< I ) -CB( I >  > *DEL < I > /CA ( I >  > *2047 . 5  > 

IF  THE  DELTA  DISTANCE  IS  LESS  THAN  .01  INCHES  THEN  WE 
MUST  SET  THE  COMPLETE  FLAG  AND  STOP  THIS  DRIVE 
33  IF<ABS<DEL<  I )) .GT.  .01 )G0  TO  35 

F  L  (  I  >  =  1  . 

I  BUF  ■:  I+4>=2048 

35  CONTINUE 

C  SET  UP  TO  DRIVE  THE  DIGITAL  TO  ANALOG  CONVERTORS  ON  THE  LPS-11 

I R A T E  < 1 ) =2 
IRATE ( 2  >  =  1 

36  IBUF<7)=IBUF<5) 

IBUF ( 8 ) =IBUF ( 6 ) 

C  OUTPUT  THE  VALUE  TO  THE  D-TC-A  DRIVERS 

CALL  SD AC ( IBUF , 8  » 1 6  t I R ATE , 7 »  0 » 2 » I  SB » 1 ) 

C  IF  AN  ERROR  WAS  ENCOUNTERED  -  TELL  THE  USER 

IF< I  SBC  1 ) .GT . 1) UR I TEC  5*907) ISB( 1  ) 

C  IF  BOTH  COMPLETE  FLAGS  ARE  SET  THEN  CHECK  FOR  ACTUAL  COMPLETION 

IF<FL< 1 >+FL<2) .EQ.2)  GO  TO  55 

C  GO  AND  SEE  IF  WE  ARE  CLOSE  ENOUGH  TO  STOP  YET 

GO  TO  15 

C  CHECK  ON  THE  FINAL  POSITION  AGAIN  AND  CLEAR  COMPLETE  FLAGS 

55  DO  75  1=1,2 

11=1-1 

CALL  ADC ( 1 1 ?  V ( I ) *  1 » I SB ) 

DEL ( I ) =- < ( V( I )/131 072.-1 . ) *  1 0 . -P  < I >  > 

75  FL ( I >  =0 . 

C  IF  BOTH  ARE  WITHIN  .01  THEN  THEY  ARE  WITHIN  LIMITS 

IF< ABSCDELCl ) ) .LT. .01 . AND. ABS<DEL(2) ) .LT. .01 )G0  TO  56 

C  AT  LEAST  ONE  IS  OUT  OF  LIMITS  -  GO  MOVE  IT  AGAIN 

GO  TO  15 
56  CONTINUE 

907  FORMAT (  '  SDAC  ERROR  'rI4) 

908  FORMAT ( '  MOVE  ADC  ERROR  '  »I4) 

C  INSURE  THAT  BOTH  DRIVES  ARE  STOPPED 

IF  (IBUF (5) .NE.2048)G0  TO  910 
IF  ( IBUFC6) .E0.2048)G0  TO  999 
910  CONTINUE 

C  AT  LEAST  ONE  DRIVE  WAS  NOT  STOPPED  -  STOP  THEM  AND 

C  GO  TO  RECHECK  THE  POSITIONS 

TBUF ( 5 ) =2048 
I BUF  <  6 ) =2043 
IBUF ( 7 ) =2048 
IBUF ( 8 ) =2048 

CALL  SDAC ( IBUF  f 8  » 16  f I RATE > 7 1 0 r 2 » I SB » 1 ) 

GO  TO  15 
999  CONTINUE 

C  BOTH  DRIVES  WERE  STOPPED  BUT  WE  WILL  INSURE  IT  BY  SETTING 

C  THE  DRIVE  VOLTAGES  TO  ZERO  ONE  LAST  TIME 

CALL  SDAC ( IBUF  »  3  > 1 6  > I RATE  *  7 »  0  »  2  » I SB  » 1 ) 

RETURN 

END 
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THIS  PROGRAM  IS  USED  TO  BUILD  THE  ’MENU. DAT' 

DATA  FILE  REQUIRED  B  V  'GRAFS’.  IF  A  MENU  FILE 
ALREADY  EXISTS  THEN  IT  READS  THAT  AND  WORKS 
FROM  THERE.  OTHERWISE  IS  ASSUMES  CERTAIN  DEFAULTS. 
THE  USER  IS  ALLOWED  TO  MODIFY  ANY  PARAMETER  SETTINGS 
BEFORE  THE  DATA  IS  REWRITTEN  TO  DISK. 


C******************* *********  ********** ******************** 


DATA  DECLARATIONS  AREA 


C*** *************** ******************************************** 

LOGICAL*l  D  L  A  B  ( 1 0  »  5 ' 

LOGICAL*l  DLABELC 10 >20) > LSCR'200) >STITLE(20) > LTITLE(oO) >  TITLE ( oO 
LOGICAL  NOCHG > EOF , MNMX , DEFALT 
DIMENSION  ID-L I  ST  <  5  )  >  DAT  (5*100)  >  X ( 50 ) 

C** ****** ***************** ************************************* 

C 

c 

C  FORMAT  DECLARATIONS  AREA 

C 

C 

Ct* *****************  *********4*?** ** *************** ***************** 


0  =  C0NTCjUR  f  1=AX0N0METR  I C  > 


FORMAT (  '  ENTER  PLOT  OPTION  I  0  =  C0NT0URf  l=AXONOMETRICr 
'  2=CR0SS  FLOW'  ) 

FORMAT(FIO.O) 

FORMAT (  '  YOU  MAY  PLOT  ANY  OF  THE  FOLLOWING  QUANTITIES: 
FORMAT (5  < IX  >  12  >  IX  > 10A1  )  ) 


FORMAT ( 


OPTION  = 
I. TITLE  = 


r  F  2 . 0  f  /  f 
>  6  0  A  1 t / t 


'  3.X  AXIS  ANGLE  =  '  > F 1 0 . 3 r / t 
'  4.Y  AXIS  ANGLE  =  '>F10.3>/> 

'  5 .  Z  AXIS  ANGLE  =  ;>F10.3>/> 

'  6.X  MIN  =  '  ,  F 1 0 . 3*>  /  t 
'  "  .  X  MAX  =  '  » F10 . 3> / f 
'  3 . Y  MIN  =  ' >  FI  0 . 3  > / > 

'  9  .  '  MAX  =  '  *  F10 . 3) 

FORMAT ( '  10.  I  VARIABLE  =  '>12,'.  M OA 
FORMAT (  '  ENTER  NUMBER  OF  PARAMETER  TO  C! 

'  VALUE  (IF  APPROPRIATE)') 
FORMAT ( I  3  >  F 1 0 . 0  > 

FORMAT ( '  ENTER  NEW  TITLE*  60  CHARACTERS 
FORMAT ( 60 A  1 ) 

FORMAT (  ll.ZMIN  =  '>F10.5>/> 

'  12.ZMAX  =  '  *  F  1  0 . 5  * / * 

'  13. ZINC  =  '  fF10.5) 

FORMAT ( '  DO  YOU  WANT  A  HARDCOPY  PLOT?') 


>  10A1  ) 

TO  CHANGE  AND 


NEW  '  * 


FORMAT ( '  D 
FORM AT ( A  1  ' 
F  0  R  M  A  T  C  '  1 


•  OPTION  =  '  >  F  2 . 0  >  /  > 

2.  TITLE  =  '  * 60A1 >/, 

5 . X  AXIS  ANGLE  =  ' >F10.3>/> 
*. Y  AXIS  ANGLE  =  '  >  F 1 0 . 3  >  /  > 
i.Z  AXIS  ANGLE  =  '  <  F  1  0 . 3  r / * 
>.X  MIN  =  '  *  F  1 0 . 3  *  /  * 

'.X  MAX  =  '  »  F  1  0 . 3  *  /  » 
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'  3. Y  MIN  =  ' .  F 1 0 . 3  .  /  . 

'  Q . Y  MAX  =  ' »F10.3»/, 

'  10. VECTOR  LENGTH  SCALE  FACTOR  =  '.F10.3) 

F ORM AT (  '  DATA  XMIN  -  '.F10.3./. 

'  DATA  XMAX  =  '.F10.3./. 

'  DATA  YMIN  =  '.F10.3./. 

'  DATA  YMAX  =  '.F10.3) 

FORMAT ( '  DATA  ZMIN  =  '.F10.3./. 

' D  DATA  ZM AX  =  ' »F10 .3) 

FORMAT ( '  14. Z  SCALE  FACTOR  -  '.F10.3) 

FORMAT  ('  DO  YOU  WANT  TO  SEE  THE  DATA  MINIMUM /'MAXIMUM'  . 

'  VALUES?  C Y/NO  '  ) 

FORMAT*  A1  ) 


FORMAT ( 

FORMAK 
FORMAT ( 


C**** *******  ************ ****** ***************  *********************** 


DATA  PRESETS  AREA 


C*  ******* ************************************************************ 

DATA  YES/'Y'/ 

DATA  MAXSIZ/lOO/ .N0/5/.LUN/3/ 

C4*t**** ************************************************  ************* 

C 

c 

C  MAIN  PROGRAM  CODE  STARTS  HERE 

C 

c 

C********* ****************** ****************************************** 

CALL  I N I TT ( 960 ) 

DEFALT*. FALSE. 

NOCHG= . TRUE . 

MNMX=. FALSE. 

EOF*. FALSE. 

CALL  ASS  I GN ( 2 . 'MENU.DAT' .8) 

CALL  FDB3EK2.  'READONLY' ) 

READ  <  2  .  END  =  200 ) TITLE . 0PTI0N.HRDCPY.A1.A2.A3.A4.S1.S2.S3.S4. 

1  XMIN . XMAX f X  INC. YMIN . YMAX. Y INC . ZMIN. ZMAX. ZINC. IDLIST 

IOPT=INT (OPTION) 

WR I TE  <  5 . 1 1 3  ) 

READ (3.119) ANS 

IF  (ANS. NE.YES>MNMX=. TRUE. 

CALL  CLOSE ( 2 ) 

IF  (MNMX)GO  TO  1000 
GO  TO  610 

C  EOF  WAS  ENCOUNTERED  ON  MENU  FILE  -  GET  DESIRED  OPTION 

C  AND  USE  THE  APPROPRIATE  DEFAULTS  FOR  THAT  OPTION 

200  CONTINUE 

DEFALT*. TRUE. 

IDLIST ( 1 )=2 
IDLISK2)  =1 
WRITE  <  5 . 10 1 ) 

READ ( 5 . 1 02 ) OPT  ION 
IOPT*INT (OPTION) 

GO  TO  <300. 400. 500)10 PT+l 
GO  TO  200 

C  CONTOUR  PLOT  OPTION  SELECTED 

300  CONTINUE 

A  1  =0  . 

A2=90 . 

A  3  *  0  . 
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S  1  =  1  . 
y  2  =  l . 

1  j  3  -  0  . 

A  A-  0  . 

3  4  -  0  . 

IDlIST < 3  )  =  10 
IDLIST(4)=3 
IDLIST (5)-ll 
GO  TO  600 

i:  AXO NOMETRIC  PROJECTION  PLOT  SELECTED 

400  CONTINUE 
A  1  =0  . 

A2= 1 20 . 

A3  =  90 . 

A4  =  0  . 

Sl  =  l  . 

S  2  =  1  . 

33  =  - 1  . 

S4  =  0  . 

IDLIST ( 3 ) =10 
IDLIST! 4)  =3 
IDLIST  <  5  )  =  1 1 
GO  TO  600 

C  CROSS  FLOU  PLOT  SELECTED 

500  CONTINUE 
A  1  =0  . 

A2=90 . 

A3  =  0  . 

A4  =  0 .  •  ~ 

S 1  =  1  . 

S2=  1  . 

S3=  .  5 
S4  =  0  . 

IDLIST (  3  )  =  1 9 
I  D  L  1 3  T  (  4  )  =  2  0 
ID LIST (5) =3 

GET  DATA  MINIMUM  AND  MAXIMUMS  FROM  THE  DATA 
FILE  AND  DETERMINE  THE  DATA  INCREMENTS  NECESSARY 
TO  PLOT  THE  DATA 
600  CONTINUE 

CALL  CLOSE ( 2 ) 

610  CONTINUE 

CALL  ASSIGNC3, 'DATAIN.DAT' r 10) 

CALL  FDBSET(3» 'OLD' >  'SHARE' ) 

DXM I N  =  999  . 

DXMAX=- 999 . 

DYMIN  =  999  . 

DYMAX=-999. 

DZMIN-999  . 

DZMAX=-999. 

700  CONTINUE 

CALL  GETGRD ( LUN  r  ND  t I DLIST » MAXSI Z » LSCR » X » STITLE » T I TLE » DLAB * 
1  DATrNFr NLf NGfNLATf EOF) 

IF  ( EOF ) GO  TO  800 

CALL  M INMAX  C DAT , DXMIN. DXM AX  »DYMIN» DYMAX , DZMIN • DZMAX  :• 

GO  TO  700 
300  CONTINUE 

CALL  CLOSE  < 3 ) 

IF  (DEFALT. EG. .FALSE. )G0  TO  1000 

XMAX=DXMAX 

YMAX= DYMAX 
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ZMAX=DZMAX 
M  I  M  =  D  X  M  I  N 
YMIN=DYMIN 
ZMIN=DZMIN 
X  M  -  X  M  I  ‘'l 
?  M  =  r  ii  l  N 
Z  M  =  Z  M  I N 

CALL  INCRHT(XM»XMAX»XINC>5.  > 

CALL  INCRMT? YH , YMAX , Y INC  -  5 .  / 

CALL  INCRMT ( ZM > ZMAX , ZINC  .  10  .  ) 

XINC=AMAX1 <  XINC r YINC) 

YINC=X INC 

XMIN=XINC*AINT?XM/XINC> 

YMIN=YINC*AINT  CYM/YINC) 

ZMIN=ZINC*AINT ( ZM/ZINC ) 

IF  (Xh.LT.O.  )XMIN=XMIN-XINC 
If  ( YM.LT.O. ) YMIN=YMIN-YINC 
XMAX=XMIN65. *XINC 
YMAX  =  YMINJ-5  .  *YINC 
IF  (ZhIN.LT .ZM>ZMIN=ZMIN+ZINC 

******** **************************************** i********#i4 


DISPLAY  F'ARAME-ERS  AND  ALLOW  CHANGES 


C********************************************** ************************ 

1000  CONTINUE 

IF  (OPTION.  GT.  1  .  )G0  TO  "2000 
CALL  NEWF'AG 
IF  (MNMX)GO  TO  1005 

WRITE?  5f 1 15) DXMINf  DXMAX  f  DYMINfDYMAX 
WRITE?5f U6)DZMINfDZMAX 
1005  CONTINUE 

WRITE(5f105)0PTI0NfTITLEfA1fA2fA3fXMINfXMAXfYMIN»vMAX 
WRITE?  5  f 106  > IDLIS7?3>  f ( DLAB? I f  3) f 1  =  1 f 10  > 

IF  (0PTI0N.LT.1.)WRITE(5f111)ZMINfZKAXfZINC 
IF  ( OPTION . EQ. 1 . 5  WRITE? 5 f 117) S3 
WR I TE ( 5  t 107) 

READ?  5  f  103)  IPARAM.  VAL.U 
IF  ( IPARAM . EQ . 0 ) GO  TO  3000 

GO  TO  (  IOIOf  1020f  1030f  10->0f  1050f  1060f  1070.  IOSOf  1150)  IPARAM 
GO  TO  ?  1090f  U20r  1130. 1140*  1160)  IPARAM-^ 

GO  TO  1000 
1010  CONTINUE 

IQPT=INT (UALU) 

IF  (  ( IOPT . LT . 0 )  . OR . ? I  OPT . GT . 2 ) ) GO  TO  1000 
OPTION=IOPT 

GO  TO  ? 300 f 400 r 500 ) IOFT+1 
1020  URITE<5f10?) 

READ?5f110)TITLE 
GO  TO  1000 
1030  A  1 =  U  A  L  U 

GO  TO  1000 
1040  A2--VALU 

GO  TO  1000 
1050  A3=UALU 

GO  TO  1000 
1060  XMIN=VALU 

NOCHG=. FALSE. 

GO  TO  1000 
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1070 

<  h  A  X  =  0  A  L  U 

N  0  C  H  G  =  .  F  A  L  S  E  . 
GO  TO  1000 

1030 

Y  fi  I  N  -  0  A  L  U 

NOC AGs. FALSE. 
GO  TO  1000 

1  0  ?  0 

CONTINUE 

I  V  =  I  N  7  <  0  A  L  U  ) 
IF  (  <  I  0 . G T . 0  > 

.AND. 

W  R  I T  E  (  5  »  1  0  3  ) 

W  R I  T  E  (  5  -  1 0  4  >  D 
GO  TO  1000 

LABEL 

10  95 

CONTINUE 

I  D  L  I  S  T  (  3  )  =  I  V 
GO  TO  610 

1120 

Z  H  I  N  -  U  A  L  U 

GG  TQ  1000 

1130 

2  !i  A  X  s  0  A  L  2 

GO  "0  1000 

1140 

ZINOVALL’ 

GO  T 0  1 0 0 0 

1150  continue: 

YnAX=UALU 
NOCHG= . FALSE . 
GO  10  1000 
1160  CONTINUE 
33- VaLO 
GO  TO  1000 
2000  CONTINUE 


CALL  Nt'J^AG 

IF  (  ri N X  :  G 0  10  2005 

WRITE  (  5 . 1 1  5  '  DXri.N  -  DXp'AX  >  DYMIN  -  DYMAX 
2005  CONTINUE 

WRITK<5-  1 14 >CPr I  ON- TITLE-  A1  ,  A2 -  A3 » XM I N - XMAX >  YMIN-YMAX. 
WHITE ( 5  f 107  ' 

R  E  A  D  (  5  -  1  0  S  )  I  P  A  K  A  M  >  0  A  L  J 
IF  ( IFAR'An  .  EG  .  0 )  GO  TO  3000 

GO  TO  ( 2010 - 2020 - 2030  *  2040 - 2050 , 2060 • 2070 - 2080 » 2090  )  I c 
IF  ( I PAR An. EO. 1O)G0  TO  2095 
GO  TO  COCO 

2010  continue; 

iopt=:nt (va_u) 

IF  ( < IOPT.LT. 0) .OR. (I0FT.GT.2) )G0  TO  2000 
OPTION= IOPT 

GO  TO  (300-400-500)  IOF'T  +  1 
2020  URITEC 5-109) 

READ (5-110)  TITLE 
GO  10  1000 
2030  A  1  =  U  A  L  U 

GO  TO  2000 
2040  A2sOAl.ll 

GO  TO  2000 
2050  A  3  =  U  A  L  U 

GU  TO  2000 
2060  XMINsUALU 

NO  ChG- . FALSE . 

GO  TO  2000 
2070  XrtAX=UALU 

N  0  C  H  G  --  .  F  A  L  3  E  . 

GO  T(1  2000 
2  0  80  Inlf.-V  h  L 


w  o 


A  R  A  M 
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NOCHG= . FALSE . 

GO  TO  2000 
2090  CONTINUE 
Y  M  A  X  =  V  A  L  U 
N  0  C  H  G  = .FALSE . 

GO  TO  2000 
2095  CONTINUE 
3  3  =  V  A  L  U 
GO  TO  2000 

C********************************* *********************  ************ 

C 


CLOSE  FILES 


WRITE 


MENU . DAT 


DISK 


3100 


C ********  ************************************************************ 

3000  CONTINUE 

IF  ( NOCHG ) GO  TO  3100 
XM-XMIN 
Y  M  =  Y  M  I  N 

CALL  INCRMT (XM.XMAX.XINC.5.  ) 

CALL  INCRMT (YM*YMAXrYINCr5. ) 

XINC=AMAX1 CXINCrYINC) 

YINC=XINC 

XMIN  =  XINC*AINT  <  XM/XINC) 

YMIN=YINC*AINT(YM/YINC) 

IF  (XM.LT.O. )XMIN=XMIN-XINC 
IF  < YM.LT.O. ) YMIN=YMIN-YINC 
XMAX=XMIN+5 .  *XINC  •  ~ 

YMAX= YMIN+5 . *YINC 
3100  CONTINUE 

CALL  ASSI3M(2» 'MENU. DAT.' »3) 

CALL  FDBSET (  2  > 'NEW' t 'SHARE' ) 

WRITE (5. 112) 

R  E  A  D ( 5 . 113)HRDCFY 

URITE( 2) TITLE » OPTION .HRDCPY.A1 rA2>A3>A4,SlrS2,S3rS4, 

1  XMIN.XMAX.XINC. YMINf YMAX  , YINC » ZM IN  »  ZMAX , Z INC  > I DL I  ST 

CALL  CLOSE ( 2 ) 

CALL  FIN ITT ( 0 » 100) 

END 

C***** ****** ******************************************************** 

c 

c 

C  THIS  ROUTINE  COMPUTES  THE  STANDARD  INCREMENT  FOR  EACH  AXIS 

C 

C 

C*** ******************************************************* ************ 

SUBROUTINE  INCRMT ( XMIN » XMAX » X  INC » COUNT ) 

DIMENSION  FRAC ( 4  ) 

DATA  FRAC/O.  , .301 0299956  » .6989700041 . 1 ./ 

XINC=(XMAX-XM IN) /COUNT 
A-ALOGIO(XINC) 

WHOLE  =  A  I N  T ( A ) 

FRACT=ABS( A-WHOLE) 

IX-1 

IF  ( FRACT . GT . FRAC ( 2 ) ) IX=2 
IF  (FRACT. GT. FRACC3) ) IX=3 
IF  (A.GT.O)  I  X  =  I  X  + 1 
A=WHOLE+SIGN(FRAC( IX) .A) 

XINC=10.**A 

RETURN 
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END 

CI##*M***t*t$****f*************#**#*#*#**##**#####***###**#  ******* 

'  C 

c 

c  THIS  ROUTINE  FINDS  THE  MINIMUM  AND  MAXIMUM  VALUES 

r  C  IN  ARRAY  DAT  FOR  X,  Y,  AND  Z. 

C 

C 

*  **********************************************  ********** 

SUBROUTINE  M I NMAX ( DAT » XM IN  *  XMAX »  YMIN  »  YMAX »  ZMIN »  ZMAX ) 
DIMENSION  DAT ( 5  *  100) 
r  DO  100  1-1*100 

XM I  N  =  AM  INI  (X  MIN  f  DAT  <1  <D) 

YMIN  =  AMIN1 ( YMIN»  DAT(2» I )  ) 

'  ZMIN=AMIN1 < ZMIN» DAT ( 3 . I ) ) 

XMAX= AMAX 1 ( XMAX  t DAT ( 1 » I )  > 

YMAX  =  AMAX1 ( YM AX ,  DAT ( 2  t I  )  ) 

C  ZMAX=AMAXl(ZMAX»DAT<3»I>) 

100  CONTINUE 
RETURN 

C  END 


J 

J 


A21 


*-ir ^j**+*r  ■ 


itvprouucuu  on  f  US  tioverr.n  copier 


m  r>  r?  n  nnnn 

nnn  in  nnnw 

ro  m  mm  r 

«  f*l  F>  f*i  ' 


n  ri  C4  c 
W  m  u“  - 


o  o-  o-  - 
ri  n  ri^ 


03  03  CO  03 
A  03(13  03 


03  (O  03  03 
03  03  03  U3 


a  A  at  tt 
ft  A  A  A 


CCi  CD  I 

A  A  A  A  A  A  I 
AAA  AAAI 


333  33333333333 
333  3333333=33  = 

3  3 
3  3 
3  3 
3  3 
3  3 

333  33333=33333 
333  33333333333 


:z  zxz  zzzzzzzz: 
:z  zzz  zzzzzzzz: 


zzzzzzzzzzz: 

zzzzzzzzzzz: 


nounu 

a  a  a  w 


A  05  05  05 

A  <n  05  cn 


(9QOQ 
©  O©  IS 

o  o 

(S  o 
IS  o 
a  o 

o  oaooooooo 

OOOCDOOCTOO© 


r 


Ct ******** ************************************************** 

c 

c 

C  THIS  ROUTINE  CHECKS  TO  SEE  IF  X  AND  Y  ARE 

C  WITHIN  THE  LIMITS  DEFINED.  IF  NOT  IT 

C  SET  THE  VISIBILITY  FLAG  TO  INVISIBLE. 

C 

C 

C************* ********************************************** 

SUBROUTINE  L I  < I TS ( DAT , NTOT , XM I N , XMAX , YM I N , YM AX , I X , ZM I N , OPT ) 
DIMENSION  DAT (  5  , 100) 

DO  100  1  =  1 »  N  T  0  T 

IF  ( (  D  A  T  (  1  ,  I  )  .LT.XMIN) .OR. (  D  A  T  ( 1 , 1 ) .GT . XMAX)  > DAT ( IX > I ) =1 . 

IF  ( ( DAT ( 2  ,  I ) . LT . YMIN ) .OR. ( BAT ( 2  ? I )  . GT. YMAX ) ) DAT ( IX ,  I )  =  1  . 

DAT (1,1) =D AT (1,1) -XMIN 
DAT ( 2 , I ) =  D  A  T (2,1 ) -YMIN 
IF  < OPT. EG. 1 . )DAT(3, 1 >=DAT<  3, 1 )-ZMIN 
100  CONTINUE 
RETURN 
END 

C*** ************************** ********************************* 

c 

c 

C  THIS  SUBROUTINE  PLOTS  A  SINGLE  CONTOUR  LINE 

C 

c 

c ********************************* ******************************* 

SUBROUTINE  CQNTORCNtXl ,? 1 , Z 1 , V 1 , X2 , Y2 , Z2 , V2 , X3 , Y3 , Z3 , L I NE > 

C 

C 

C  input: 

C  N  NO.  OF  VALUES  IN  COORDINATE  ARRAYS 

C  X 1 , Y 1 , Z 1  INPUT  COORDINATES  OF  ONE  PAIL 

C  VI  VISIBILITY  INDICATOR  OF  ONE  RAIL 

C  X 2 , Y 2 , Z 2  INPUT  COORDINATES  OF  SECOND  RAIL 

C  V2  VISIBILITY  INDICATOR  OF  SECOND  RAIL 

C  X3 ,  Y  3  SCRATCH  ARRAYS 

C  Z3  Z  VALUE  OF  CONTOUR  LINE 

C  LINE  LINE  STYLE  PARAMETER'  FOR  CONTOUR  LINE 

C 

c 

DIMENSION  X 1 <  N ) , Y 1 (N) ,Z1<N),V1(N),X2(N) »  Y  2 ( N ) , Z  2 ( N ) , V  2 ( N ) , X  3 ( N  )  , 
1  Y3(N) 

LOGICAL  FIRST 

Firsts. false. 

I\=  1 
N 1  =  1 
N2  =  1 
IS=1 

IF( <  Z 1 < 1 >-Z3)*<Z2( 1 )-Z3) .GT.O. )G0  TO  10 
CALL  INTER (XI , Y1 ,Z1 ,V1 ,X2, Y2,Z2,V2,X3, Y3,Z3) 

K  =  K  +  1 

FIRST= . TRUE . 

10  CONTINUE 

CALL  CHANGE (N1 ,N»Z1 ,Z3) 

CALL  CHANGE (N2,N,Z2»Z3) 

IF (FIRST )GQ  TO  40 
20  CONTINUE 

I3=MIN0(N1 ,N2) 

I F ( IS. EG. 99?) RETURN 
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isi=is-i 

IF  C IS.EQ.N1 )G0  TO  30 

CALL  INTER(X2( IS1 ) »Y2( IS1 ) »Z2< IS1 ) » V2< IS1 ) »X2( IS) . Y2( IS) »Z2<  IS)  f 

*  V2< IS) »X3(K) t Y3(K) »Z3) 

K  =  K  + 1 

CALL  CHANGE(N2>NfZ2»Z3) 

GO  TO  40 
30  CONTINUE 

CALL  INTER (XI ( IS1 ) * Y1 (IS1) »Z1 ( IS1 ) « VI < IS1 ) .XI ( IS) » Y1 < IS ) » Z1 < IS ) » 

*  VI  (  IS)  »X3(K>  f Y3(K) »Z3) 

K-K+l 

CALL  CHANGE(N1fNfZ1fZ3) 

40  CONTINUE 

I F  =  M I  NO ( N 1 f  N  2 ) 

I F 1= IF- 1 

I F ( IS.GT. I F 1 ) GO  TO  55 

DO  50  I  =  1 3  f I F 1 

IF ( I . GT . N ) GO  TO  70 

CALL  INTER ( XI < I)fY1<I)fZ1<I)fV1<I)fX2<I)fY2<I)fZ2(I)fV2<I>fX3(K)- 
1  Y3<K)fZ3) 

K  =  K  + 1 

50  CONTINUE 

IF( IF.EQ.N1 )G0  TO  60 
55  CONTINUE 

CALL  INTER  (X2<  IF1  )  f  Y2  ( IF1 )  f  Z2  ( IF  1 )  f  V2  (  IF  1  )  f  X2  (  I  F  )  ,  Y  2  (  I F  )  .  Z2  'I  F  )  ». 

*  V2(  IF) f  X3  <  K ) »Y3(K) f  Z3 } 

CALL  CHANGE(N2fNfZ2fZ3) 

GO  TO  80 

60  CONTINUE  •  " 

CALL  INTER  (  XI  (  IF1)fY1<IF1)fZ1<IF1.)fV1<IF1)fX1<IF)fY1<IF)fZ1(IF)f 

*  VI < IF) fX3<K) fY3(K) fZ3) 

CALL  CHANGE(N1»N»Z1»Z3)  . 

GO  TO  80 

70  N=N- 1 

80  CONTINUE 

I  STAR T  =  1 
DO  90  LLL= 1 t K 
KK=LLL 
H=X3(KK) 

V  =  Y  3  (  K  K  ) 

I F ( < I  START. EQ. 1 ) .AND. (X3(KK) .NE.-999 . ) ) CALL  DASHA <Hf V» -1  ) 

IF ( ( ISTART.EQ.O) .AND. ( X3 ( KK ) . NE . -999 . ) ) CALL  DASHA  <  H  f  V  f  L I NE ) 
ISTART=0 

IF(X3<KK) .EQ.-999. ) I START= 1 
90  CONTINUE 

K  =  1 

GO  TO  20 
END 

C***** ***********  ************ ******* ******** **##**  +  ****** ******* 

C 

c 

C  THIS  ROUTINE  LINEARILY  INTERPOLATES  BETWEEN  THE  INPUT  POINTS 
C  FOR  THE  POINT  WITH  THE  SPECIFIED  Z  VALUE 

C 

c 

c*********************** ******************************************* 

SUBROUTINE  INTER < X 1 f Y 1 f Z 1 t V 1 f X2 f Y2 f Z2 t V2 f X3 » Y3 f Z3 ) 

C 

C 

c  input: 

C  XliYlrZl  INPUT  COORDINATES  OF  THE  FIRST  POINT 
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C  V 1  visibility  parameter  of  the  first  point 

C  X2 • Y2  »  Z2  INPUT  COORDINATES  OF  THE  SECOND  POINT 

C  V  2  VISIBILITY  PARAMETER  OF  THE  SECOND  POINT 

C  23  SPECIFIED  Z  VALUE 

C  OUTPUT 

C  X  3  »  Y  3  >  Z  3  OUTPUT  COORDINATES  OF  THE  POINT 

C 

P=t 

IF ( Z2 . EQ . Z3  >  GO  TO  10 
F'  -  (  Z  3  -  Z  1  )/  CZ2-Z1  ) 

10  CONTINUE 

X  3  =  X  1  +  P  *  (  X  2  -  X  1  / 

YJ  =  Y1tF’*(  Y2-Y1  > 

IF ( (VI .NE.O.  )  .OR.  (V2.NE.0. ) ) X  7  ~ - 9  9  9 . 

RETURN 
END 

i:t*4**t*t***t**#Mttt*Mt«******t**M**********MM******M*******M 

c 
c 

c  THIS  ROUTINE  DETERMINES  WHERE  THE  DATA  IN  THE  Z1  ARRAY  CROSSES 
C  THE  Z3  VALUE 
C 
C 

Ct*******t**M*****MM*M*M*M*Mt****MM**MM**MM**MtM*Mt**t 

SUBROUTINE  CHANGE(NlfNfZlfZ3) 

C 

c 

c  input: 

C  N 1  LOCATION  IN  ARRAY  OF  THE  PREVIOUS  CROSSING 

C  N  TOTAL  NO.  OF  VALUES  IN  THE  Z1  ARRAY 

C  Z1  INPUT  ARRAY 

C  Z3  CROSSING  VALUE 

C  output: 

C  N  LOCATION  IN  ARRAY  OF  THE  NEXT  CROSSING 

C 

DIMENSION  Z1(N) 

IF ( N 1 . GE . N  >  GO  TO  20 
N 1 1 =N 1 + 1 
DO  10  I  =  N 1 1 »  N 
N 1  =  I 

IF((Z1(I-1)-Z3)»(Z1(I)-Z3) .T_E . 0 . ) RETURN 
10  CONTINUE 

20  N 1 =999 

RETURN 
END 

C*M***M***************«*********M**M***********M***»***M******M 

C 
C 

C  THIS  ROUTINE  DISPLAYS  A  LEGEND  DEFINING  WHAT  VALUES  THE  DIFFERENT 
C  LINE  STYLES  REPRESENT 

C 
C 

C4M*****4MM*4**4*MM****M*MMM*M*4  ****♦**♦*♦♦#*♦**♦♦**♦♦♦♦#****** 

SUBROUTINE  LEGEND < N f L INE f ZLOU f ZINC f XS f YS f X  INC f Y INC > 

C 

C 

C  input: 

C  N  NO.  OF  LINE  STYLES  TO  DISPLAY 

C  LINE  LINE  STYLE  PARAMETER  ARRAY 

C  ZLOW  Z  VALUE  FIRST  LINE  STYLE  REPRESENTS 

C  ZINC  Z  INCREMENT  BETWEEN  LINES 
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C  X  S  »  Y  S  COORDINATES  OF  LEFT  SIDE  OF  FIRST  LINE 

XINC  COORDINATE  OF  LEFT  SIDE  OF  SIXTH  LINE  IF  NECESSARY 

YINC  COORDINATE  DELTA  BETWEEN  LINES 

DIMENSION  LINE(N) 

LOGICAL*l  NADE(IO) 

Z3=ZLQW 

Y  =  YS 

X  =  XS 

I CNT=0 

DO  30  1  =  1, N 

CALL  DASHA ( X , Y , -1 ) 

CALL  DASHA < X+X INC/2 . , Y , L INE  < I  )  ) 

ENCODE* 10,99, NADE)Z3 

NADE ( 1 ) =32 

CALL  AANSTR* 10, NADE) 

Z3=Z3+ZINC 
Y=Y-YINC 
ICNT= I CNT  +  1 

30  CONTINUE 

RETURN 

9  FORMAT ( F 1 0 ♦ 3 ) 

END 

************************************** ******* ***************** ****** 


THIS  ROUTINE  DISPLAYS  A  GRID  BY  DRAWING  THE  RUNGS  FIRST  AND  THEN 
THE  SECOND  RAIL 


***************** ******* * ****** * ***********  ************** ************** 

SUBROUTINE  D3MESH ( F , N , X 1 , Y 1 , Z 1 , W 1 , X2 , Y2 , Z2 , W2 ) 


INPUT  : 

F  PROJECTION  MATRIX 

N  NO.  OF  VALUES  IN  COORDINATE  ARRAYS 

X 1 , Y 1 , Z 1  COORDINATES  OF  FIRST  RAIL 
W1  VISIBILITY  PARAMETERS  FOR  FIRST  RAIL 


X2,Y2,Z2  COORDINATES  OF  SECOND  RAIL 
W2  VISJ8ILTIY  PARAMETERS  FOR  SECOND  RAIL 


rr  * 


DIMENSION  F ( 4 , 2 ) ,X1(N) , Y 1 <  N ) ,Z1(N> ,W1(N) , X  2 ( N ) , Y  2 ( N  > ,  Z  2  (  N  )  , W  2  <  N ) 
HI  <  I  >  *F  <  1  , 1  >  #X  1  ( I )  +  F  <  2 , 1 )  # Y  1  ( I  >  +  F  <  3 , 1 )  *Z  1  ( I )  +  F  (  4 , 1  ) 
V1(I)=F(1,2)*X1(I)+F(2,2)*Y1(I)+F(3,2)*Z1(I)+F<4,2) 

H2( I )=F< 1 , 1 >#X2< I )+F(2, 1 >*Y2(I)+F(3, 1 >#Z2< I ) +F( 4 , 1 5 
V2< I ) =F( 1 , 2>*X2( I >+F(2,2)#Y2< I )+F( 3*2>#Z2( I )+F<  4 , 2 ) 

DO  50  LLL= 1  ,  N 
I  •=  L  L  L 

IF( ( W1 ( I ) .NE .0. ) .OR. < W2< I ) .NE . 0. > >G0  TO  50 
H  =  H 1 <  I  ) 

V  =  V1 (  I  ) 

CALL  MOVEA  <  H , V ) 

H  =  H2< I  ) 

V  =  V2<  I  ) 

CALL  DRAWA ( H , V ) 

CONTINUE 
I3TART=1 
DO  75  L L L  =  1  ,  N 
I  =  L  L  L 


ivc^i  vuuct-u  01121  iiovLTiiineiil  copier 


n  n  ~  o  n 


~  V  2  <  I  ) 

IF  < ( I  START . EQ. 1 ) . AND. <U2( I ) .EQ.O , ) )CALl  MOVEA(H> V) 
IF((ISTART.EQ.O) .AND . ( W2 < I ) . E Q  .  0 .  >  >  CALL  DR AW A (H» V ) 

1ST  ART  =  0 

IF ( U2 ( I > .NE.O.  )I3TART=1 
75  CONTINUE 

RETURN 
END 

C********************************************************** ******* 


THIS  ROUTINE  DEFINES  AN  AXONOMETRIC  PROJECTION  MATRIX 


C****** ************************************************************ 

SUBROUTINE  FFILL(FfXMIN>XMAX»YMIN>YMAXrZMIN»ZMAXf 
1  SXMIN » SXMAX , SYMIN , SYMAX ) 

C 
C 

c  input; 

C  F 

c 

C  XMINrXMAX 

C  YMIN  f  YMAX 

C  ZMINtZMAX 

C  SXMIN » SYMI N  MINIMUM  SCALE  HORIZONTAL  AND  VERTICAL  VALUES  c 

C  SXMAX  t SYMAX  MAXIMUM  SCALE  HORIZONTAL  AND  VERTICAL  VALUES  " 

C  | 

DIMENSION  F  (  4  >  2  )  .  A<3)-  *  5 

H<  X  r  Y  >Z)=X*F( 1 , 1 ) +Y*F( 2 »1)+Z*F(3»1)+F(4»1)  I 

V(X,Y,Z>=X*F<1 ,2>+Y*F(2>2>+Z*F<3.2)+F<4»2) 

DO  5  1=1 > 3  _  o 

A(I)=F(I»l)/57.2^6  1 

F< I » 1 )=COS( A< I )  )¥F< I »2) 

F  (  I  >  2 ) =S IN ( A ( I ) ) #F  < I >  2  > 

5  CONTINUE 

H1=H(XMIN»YMIN»ZMIN) 

V1=V(XMIN, YMINf ZMIN) 

H2=H(XMIN»YMIN»ZMAX) 

V2=V(XMIN» YMIN»ZMAX) 

H3=H(XMIN»YMAX»ZMIN) 

V3-V(XMIN>YMAX»ZMIN) 

H4=H(XMAXr YMINrZMIN) 

V4  =  V  <  XMAX  t YMIN  t ZMIN ) 

H5=H(XMIN» YMAXfZMAX) 

V5=V(XMIN»YMAX»ZMAX) 

H6=H(XMAX» YMAXf ZMIN) 

V6=V(XMAXf YMAXfZMIN) 

H7=H(XMAX» YMINrZMAX) 

V7=V(XMAX> YMINfZMAX) 

H3=H(XMAXrYMAXrZMAX) 

V8=V < XMAX f YMAX rZMAX) 

SXMIN»AMIN1 ( HI rH2*H3»H4  »H5»H6f H7»H8) 
SXMAX=AMAX1(H1»H2>H3»H4jH5>H6jH7>H8) 

SYMIN=AMIM (VI >V2»V3»V4»V3»V6»V7»V8) 
SYMAX=AMAXl(VlfV2fV3»V4»V5»V6rV7»V8) 

RETURN 

END 

C****** ************************************************************ 


INCLUDES  3  ANGLES  AXES  MANE  WITH  THE  HORIZONTAL  AND  3 
SCALE  FACTORS 

LIMITS  OF  COORDINATE  DATA  THAT  WILL  BE  PLOTTED 
SAME  AS  ABOVE  \ 

SAME  AS  ABOVE  J 
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THIS  ROUTINE  DRAWS  AND  LABELS  A  GRID  SPACE 
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SUBROUTINE  GR I D ( F  *  XM I N *  XM AX  , X I NC  ,  YM I N  r  Y H AX  »  Y I NC ) 

DIMENSION  F  (  4  »  2  ) 

LOG I CAL  # 1  NADE(IO) 

H  <  X  »  Y  ) =  X*F (1,1)  +  Y*F<  2, 1  ) 

V'X,Y)=X*F( 1 , 2 ) + Y#F  <  2 , 2 ) 

HH=H(XMIN,YMIM) 

VV=V(XMIN,YMIN) 

CALL  MOVE A  <  HH , VV ) 

HH  =  H ( XMAX , YMIN ) 

VV=V(XMAX,YMIN) 

CALL  DRAWA ( HH , VV ) 

HH=H ( XMAX , YM AX ) 

V V  =  V  <  XMAX , YMAX ) 

CALL  DRAWA  <  HH , VV  > 

HH=H(XMIN,YMAX) 

VV  =  V  ( XMIN , YMAX ) 

CALL  DRAUA(HH,VV) 

HH=H(XMIN, YMIN) 

VV  =  V<  XMIN , YMIN) 

CALL  DRAWA  <  HH , VV ) 

X  =  XM  IN 

Y  =  YMIN-.  01 5*  (YMAX- YMIN) 

DO  20  1=1,6 

HH=H ( X , Y ) 

VV=V(X,Y)  •  " 

CALL  MOVE A ( HH , VV ) 

HH=H(X,YMIN) 

VV=V ( X , YMIN ) 

CALL  DR AW A  <  HH , VV ) 

X=X+.2*<XMAX-XMIN) 

20  CONTINUE 

Y  =  Y  M  I N 

X=XMAX+.015*(XMAX-XMIN) 

DO  30  1=1,6 
HH=H ( X , Y ) 

VV= V ( X , Y ) 

CALL  MO VE A ( HH , VV ) 

HH=H ( XMAX , Y ) 

VV= V ( XMAX , Y ) 

CALL  DRAWA ( HH , VV ) 

Y=Y+.2*( YMAX-YMIN) 

30  CONTINUE 

ENCODE ( 10,99, NADE)XMIN 
NADE ( 1 ) =32 

X=XMIN-.045#(XMAX-XMIN) 

Y= YMIN-. 05* (YMAX- YMIN) 

HH  =  H  <  X , Y ) 

VV=V ( X , Y ) 

CALL  MOVE A  <  HH , VV ) 

CALL  AANSTR ( 10 , NADE  > 

ENCODE ( 10, 99, NADE) XMAX 
NADE ( 1 ) =32 

X=XMAX-.05#(XMAX-XMIN) 

HH=H ( X , Y ) 

VV=V ( X , Y ) 

CALL  MOVEA<HH,VV) 

CALL  AANSTR ( 1 0 , NADE ) 
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ENCODE ( 10f9?fNADE> YMIN 
1  NADE ( 1  )  =  3  2 

X-XMAX  +  .  05*  (  '<  M  A  X  -  X  M  IN  ) 

HH=H(XfYMIN) 

UV  =  U ( X. YMIN  > 

CALL  M  0  U  E  A  H  H  f  U  > 

CALL  AANSTRf 10. NADE' 

ENCODE  (  10  »  99  f  NADE  )  YMAX 
NADEC 1 ) =32 
HH  =  H  < X  r YMAX  ) 

yy  =  y  < x , ymax  > 

CALL  MOVEA(HHfUU) 

CALL  AANSTR< IOfNADE) 

99  FORMAT ( F 1 0 . 3 ) 

RETURN 

END 

C**************************************************************** 

C 

C 

C  THIS  ROUTINE  PLOTS  A  UECTOR  FIELD 

C 

C 

C****************************** *********************************** 

SUBROUTINE  VF IELD < N f DAT f SCALEF ) 

DIMENSION  DAT  <  5  f 1 00  > 

DO  10  1  =  1  fN 

I F ( D  A  T  <  5  f I ) .NE.O. >G0  TO  10 
CALL  MOVEACDATCIfI) .DAT(2fI) ) 

UU=DAT (3f I ) YSCALEF  -  “ 

VV=DAT ( 4f I >#3CALEF 
CALL  ARROW ( -UUfUU) 

10  CONTINUE 

RETURN 
END 

Cl* ********************************* ******************************* 

c 

c 

C  THIS  ROUTINE  DRAWS  THE  ARROWS  FOR  A  UECTOR  FIELD  PLOT 

C 

c 

c******** *******  ********************  ******************************* 

SUBROUTINE  ARROW(U.U) 

DATA  S2fC2/.707f-.707/ 

CALL  DRAUR(UfV) 

W=SQRT(U*U+U*U) 

IF ( W . LE . .01 )RETURN 

S1=V/W 

C1=U/W 

UU=W/5.*(C1*C2-S1*S2) 

VU=W/5.*<S1*C2tC1*S2) 

CALL  DRAUR(UUfUU) 

RETURN 

END 

C *********  **************************** ********** *************** 

C 

C 

C  THIS  ROUTINE  OUTPUTS  AN  ALPHANUMERIC  STRING  TO  THE  SCREEN 

C 

C 

C**************************************************************** 

SUBROUTINE  AANSTR(NfSTRING) 

A29 
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L0GICAL*1  STRING(N) 

DIMENSION  N A D ( 6 0 ) 

DO  10  1=1 fN 

10  N  A  D  ( I )  =  3  T  R I N  G  ( I  ) 

CALL  A  N  S  T  R  (  N  *  N  A  D  > 

RETURN 

END 

C«*t*** ******* *********************************************  ******** 

C 

c 

c  THIS  ROUTINE  FILLS  ARRAYS  USED  FOR  PLOTTING  FROM  THE 

C  DAT  ARRAY 

C 

c 

c ******************************************************************* 

SUBROUTINE  F I LLI T  < NP  » K 1 » K2 f XI » Y 1 » Z1 » V 1 » X2  » Y2 r Z2  * V2 >  DAT ) 
DIMENSION  Xl(NF')  »  Y 1  (  N  P  )  »  Z 1  (  N  P  )  »  V 1  (  N  P  )  »  X  2  (  N  P  )  rY2(NF')  r 
1  Z2(NP) .V2(NP)»DAi  (5fl00) 

DO  115  J  J  = 1 »  N P 

XI ( JJ ) =DAT ( 1 t K1  + JJ-1 ) 

Y 1 ( JJ) =DAT ( 2 .Kl+JJ-1  ) 

Z1 C JJ)=DAT  <3»N1+JJ-1 ) 

VI  (  JJ  )  =DAT(  4  f  M+JJ-1  ) 

X  2 ( J  J  >  =  D  A  T ( 1 »N2  +  JJ-1  ) 

Y  2 ( J  J ) =  0  A  T (2*K2+JJ-1  ) 

Z2< JJ)=DAT(3fK2+JJ-1 ) 

V2( JJ)=DAT(4»K2  +  JJ-1  ) 

115  CONTINUE 
RETURN 
END 
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THIS  SUBROUTINE  IS  USED  TO  CONVERT  DATA  TAKEN  WITH 
A  SEVEN  HOLE  PROBE  INTO  USEFUL  QUANTITIES.  THIS 
PROGRAM  IS  WRITTEN  TO  BE  USED  WITH  A  PROBE  THAT 
HAS  BEEN  CALIBRATED  FOR  HIGH  SPEED  COMPRESSIBLE 

flow,  ’’he  calibration  coefficients  are  stored  in 

A  FILE  NAMED  •  CPROBE . DAT  * .  IT  IS  THE  USERS  JOB 
TO  ACCESS  THIS  FILE  AND  PASS  THE  CALIBRATION 
COEFFICIENTS  TO  THE  SUBROUTINE.  THE  FILE  MAY  BE 
ACCESSED  AS  FOLLOWS: 


CALL  ASSIGNED.  CPROBE. DAT'>10> 

CALL  FDBSETI 2f  'READONLY ' »  f  f 7) 

DEFINE  PILE  2 ( 7 . 256 . U . I REC ) 

I R  E  C  =  1 

DO  100  1-1.7 

READ<  2 ' IREC  >  <KA< I f J ) . J  =  1 .20 ) . <KB( I . J) . J  =  1 >20  >  » 
1  (KQ( I »  J  >  >J=1f20) f (KZ(IfJ) >J  =  1f20) f 
1  ( KM ( I >  J ) »  J  =  1 » 20  ) 

100  CONTINUE 


C 


SUBROUTINE  CGMC AL ( I  NT  .  X  . KALPH A  > K BET A >  KQ > KZERO >  KM  f  PATH  »  ALPHA  >  BETA 

1  ALPHAT  fBETATf  THETA  >  PHI  f  CTOT  f  C3TAT  >  CA  f  CB  fCQfCZERO  >  CD  YN>  'I 

2  CMfRMACHfUfV)  _  | 

C  INT  IS  NO  LONGER  USER-  *  5 

C  X  IS  THE  INPUT  DATA  ARRAY  WITH  PRESSURES  IN  IT  § 

C  KALPHAf  KBETA.  KQ  f  kZERO.  AND  KM  ARE  THE  PROBE  CALIBRATION  COEr  ; 

‘C  PATM  IS  THE  ATMOSPHERIC  .PRESSURE  IN  PSIA  j 

C  ALPHAf  BETA.  ALF'HATf  BETATf  THETAf  AND  PHI  ARE  THE  ANGLES  5 

C  RETURNED  BY  THE  SUBROUTINE 

C  CTOT  AL i  C3TATIC.  CZERO.  CQ.  AND  CDYN  ARE  THE  PRESSURE  COEFFIC I  E  N  T 

C  RETURNED  BY  THE  SUBROUTINE 

C  CA  AND  CB  ARE  THE  ANGLE  OF  ATTACK  COEFFICIENTS  RETURNED 

C  RMACH  IS  THE  MACH  NUMBER  RETURNED  BY  THE  SUBROUTINE 

C  U  AND  V  ARE  THE  RELATIVE  VELOCITY  VECTORS  RETURNED 

C 
C 

REAL  KALPHAfKBETAfKQfKZEROfKM 

DIMENSION  X(  12)  fKALF'HA(  7>  20)  >KBETA(7f20  )  >KQ(  7  >  20)  f  KZEROv  7  >  20  '  » 

1  I  CHAN  <  9 ) >P( 9) f ICA( 20  >  f ICB(20) > ICM<20) >KM<  7f20> 

C  PRESETS  TO  DETERMINE  WHICH  PRESSURE  IS  WHICH  INDEX  INTO  ARRAY  X 

DATA  ICHAN/2. 3 . 4 . 5. 6 . 7 . 10 > 8 r 9/ 

C  POWERS  FOR  CA  IN  THE  POWER  SERIES  EXPANSION 

DATA  ICA/Of1fOfOf2fOfOf1fOf1f3fOfOf2f2f1fOf1fOf1/ 

C  POWERS  FOR  CB  IN  THE  POWER  SERIES  EXPANSION 

DATA  ICB/OfOf 1 >0f0f2f0f 1f1f0f0f3f0f1 >0f2f2f0f1f 1/ 

C  POWERS  FOR  CM  IN  THE  POUER  SERIES  EXPANSION 

DATA  ICM/0f0.0»1f0.0f2f0f 1f1f0f0f3f0f1f0f1f2f2f1/ 

C  CONVERSION  CONSTANT  FOR  RADIANS/DEGREES 

DATA  RAD/57.296/ 

C  INTEGER  FUNCTION  TO  DETERMINE  SECTOR  NUMBER 

ISECt J)=MOD( J-1f6)+1 

C  FUNCTION  TO  FIX  RANGE  OF  ANGLES  TO  -90  TO  270 

RANGE':A)=AM0D(Af90.  f360.  )-90. 

C  INITIALIZE  THE  HIGH  PRESSURE  VALUE  TO  A  VERY  NEGATIVE  NUMBER 

PH=-1E20 

C  PROCESS  EACH  PRESSURE 

A3  3 


DO  5  1  =  1,7 

CONVERT  THE  PRESSURE  TO  ABSOLUTE  -  F'SIA 

P(I)=X(  I C  H  A  N  ( I  )  )  +  P  A  T  M 

FIND  THE  MAXIMUM  PRESSURE  VALUE 

P  H  =  A  M  A  X  1  (  P  H  •  P  ■:  I  '  ' 

CONVERT  TOTAL  PRESSURE  TO  ABSOLUTE 

(3;  *x;  I c h a n  ■  3 >  :•  +  patm 
CONVERT  STATIC  pressure  to  absolute 
P  (  9  )  =  X  (  ICHAN (?)  ) +  PATM 

IF  THE  HIGH  PRESSURE  IS  NOT  FROM  HOLE  7  SKIP  THIS  SECTION 
IF(PH.NE»P(7) >GQ  TO  20 

THE  HIGH  PRESSURE  IS  FROM  SECTOR  SEVEN  -  COMPUTE  COEFFICIENTS 
I SECT=7 

COMPUTE  AVERAGE  PRESSURE 
P  A  V  G  =  <  P ( 1)+P<2)+P(3)+P(4)+P<5)+P(6) )/6. 

COMPUTE  THE  VALUES  TO  DETERMINE  CA  AND  CB 
CA1  =  <P(4)-P(  1) )/ (FH-PAVG) 

CA2*(P(3)-P<6> ) / ( PH-P AVG  ) 

LA3  = ( P ( 2 ) -P  <  5 )  )/ (PH-P AVG) 

COMPUTE  CA 
CA=CAl+(CA2-CA3>/2. 

COMPUTE  CB 

CB= . 57735* < CA2+CA3 ) 

COMPUTE  CM 
CM=(P(7)-PAVG)/P(7) 

SKIP  THE  NEXT  SECTION  UHICH  IS  FOR  THE  OUTER  SECTORS 
GO  TO  50 
DO  25  I 6 

DETERMINE  WHICH  OUTER-SECTOR  HAS  THE  HIGH  PRESSURE 

IF ( P ( I ) . EQ . PH ) GO  TO  30 

CONTINUE 

SET  THE  SECTOR  NUMBER 
ISECT  =  1 

COMPUTE  AVERAGE  OF  THE  TWO  HOLES  ON  EITHER  SIDE  OF  HIGH  PRES. 
PAVG=<  P < I  SEC ( I  + 1 ) ) +P ( I  SEC  < 1  +  5) >  ) /2  . 

COMPUTE  CA 

C  A  =  (  P  (  I  >  -  P  (  7  )  )  /  (  P  (  I  )  -  F'  A  V  G  > 

COMPUTE  CB 

CB=(P<  ISEC(  1+5  )  )  -P<  ISEC<  I  +  1-)  >  )/(P(  I )  -P  AVG  ) 

COMPUTE  CM 

CM=(P< I)-PAVG)/P<  I  ) 

COMPUTE  ALPHA,  BETA,  CO,  CZERO,  AND  MACH  FROM  CALIBRATION 
ALPHA=0 . 

BETA=0  . 

C0  =  0  . 

CZER0=0  . 

RMACH=0  . 

USE  THE  POWER  SERIES  EXPANSION  AND  THE  COEFFICIENTS  TO 
COMPUTE  THE  ACTUAL  VALUES 
DO  60  1=1,20 

COEFF  =  CA**ICA( I)*CB**ICB(I> *CM**ICM (  I  ) 
ALPHA=ALPHA+CQEFF*KALPHA< ISECT, I ) 

BETA  =  BETA  +  COEFF*KBETA  •'  ISECT  ,  I  ) 

CQ  =  CQ  +  COEF-*!\Q(  ISECT«  I 

CZERO  =  CZEROrCOEc’F  >K  ZERO':  ISECT  ,  I  ) 

RMACH  =  RMACH  +  CQEFF*KM ( ISECT,  I  ) 

CONTINUE 

COMPUTE  TOTAL  PRESSURE  COEFFICIENT 

CTOT  =  -CZERO*(PH-PAVG)/(F (3)-P(?) )+(PH-P(3) >/(P: 3>-F (?)  ) 
COMPUTE  STATIC  PRESSURE  COEFFICIENT  A3. 

CSTAT=-CZERO*(PH-PAVG)/(P(S)-P(?> ) - ( PH-PAVG > / < CCU 
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c 

c 

c 


1 

100 

c 
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c 
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(  P  (  8  )  -  P  <  9  )  )  )  r  (  FH-P  (  9  )  )/(P(8)-P(9)  ) 

COhPUTE  DYNAMIC  PRESSURE  COEFFICIENT 
CDYN=CTOT-CSTAT 

IF  THIS  IS  AN  OUTER  SECTOR  THEN  SKIP  THIS  SECTION 

I  F ( P H  .  N E  .  P  C  7  >  ) GO  TO  30 

COMPUTE  ALF’hAT 

ALPHA f  =  A  L  P  H  A 

COMPUTE  BETAT 

betat=beta 

C  0  M P U  T  E  BETA 

BETA=ATAN < COS ( ALP HA /RAD ) *TAN ( BETAT/RAD ) >  *RAD 
COMPUTE  PHI 

IF CALPHAT.NE.O. )F'HI=AATAN(SIN(BETAT/RAD)*COS(ALPHAT/RAB) * 
COS(BETAT/RAD)*SIN<ALPHAT/RAB>)*RAD 
IF(AL P HAT. LT.O.)PHI=PH 1+180. 
IF(ALPHAT.EQ.0..AND.BETAT.GE.0.)PHI=90. 

IF (ALPHAT.EO.O.  .AND.BETAT.LT  .0.  )PHI=-90. 

COMPUTE  THETA 

THETA=a:aN(5GRTv(TAN(ALPHAT/RAD))**2+iTAN(BETAT/RAH>)**2))*RAD 

SKIP  THE  SECTION  FOR  OUTER  SECTORS 

GO  TO  100 

COMPUTE  THETA 

THET  A  =  ALF'HA 

COnPUTE  PHI 

PHI=BETA 

COhPUTE  ALPHAT 

ALPHAT=AATAN(SIN<THETA/RAD)#COS(PHI/RAB)fCOS<THETA/RAD))*RAD 
COMPUTE  BETAT 

BETAT  =  AATAN <  S I N <  THET A/RfiD ) *SI N  <  PHI /RAD ) f COS < THETA/RAD ) ) *RAD 
COMPUTE  ALPHA 
ALPHA=ALPHAT 
COMPUTE  BETA 

BETA=AATAN< COS  ( ALPHA/RAD ) *S IN < THET A/RAD ) *S IN C PH  I /RAD  )  f 
COS ( THETA/RAD ) ) *RAD 
CONTINUE 

CHECK  THE  RANGE  OF  ANGLES  FOR  ALPHAf  BETA  *  ALPHATf  BETAT*  THETA* 
ALPHA=RANGE< ALPHA) 

BET A  =  R ANGE ( BETA  ) 

ALPHAT=RANGE( ALPHAT) 

BETAT= RANGE (BETAT) 

THETA=RANGE( THETA) 

PHI “RANGE ( PHI ) 

COMPUTE  THE  ANGLE  TANGENTS  FOR  USE  IN  VELOCITY  VECTOR  COMPUTATION 
TAN AL=S I N ( ALPH AT/57 .296 ) /COS ( ALPH AT/57 . 296 ) 

TANBE  =  SIN<  BET AT/37 . 296 ) /COS < BET AT/57 . 296 ) 

INSURE  THAT  CDYN  DOES  NOT  EXCEED  LIMITS 

IF  (CDYN. LT.-.9999)CDYN=-. 9999 

COMPUTE  THE  LENGTH  OF  THE  VELOCITY  VECTOR 

UVR  =  SORT(  (CDYN+1 . )/( l.+TANAL*TANAL+TANBE*TANBE)  ) 

GET  THE  U  COMPONENT  OF  THE  VELOCITY  VECTOR 
U=UVR*T  ANBE 

GET  THE  V  COMPONENT  OF  THE  VELOCITY  VECTOR 

V=UVR*TANAL 

RETUR- 

END 

FUNCTION  A A  T  AN ( T  *  B ) 

ATAN=ATAN2(TfB) 

I F  <  T , GE . 0 . .AND.B.GE.O. )AATAN  =  ATAN 
IF( T.GE . 0 . . AND. B .LT.O. ) A A T AN  =  3 . 1 4 1 5 9  +  AT AN 
IF (T.LT.O. .AND.B.GE.O. ) A A T AN =2 . * 3 . 1 4 1 59+ AT AN 
IF (T. LT.O. .AND. B. LT.O. ) A A  TAN  =  3 . 14159  +  ATAN 

A35 
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RETURN 

END 

FUNCTION  T AN ( A  ) 
TAN=SIN( A)/COS< A) 
RETURN 
END 
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CM******* ******* ************* ******************************************* 

c 

c 

C  ROUTINE  TO  GET  A  GRID  OF  DATA  FROM  A  FILE  IN  THE  STANDARD 

C  DFAN  DATA  FILE  FORMAT. 

C 

C 

C ************************************* ************************  ********* 

SUBROUTINE  GETGRD ( LUN » ND  » I DL I  ST , MAXSI Z f LSCR » RSCR r ST  I TLE » LT I TLE  » 

1  DLABEL  >DAT  ,  NP>NL>NG»  NLAT.EOF  ) 

0 

LUN 
ND 
NP 
NL 
NG 

NL  A  T 
DLABEL 
DAT 

STITLE 
LTITLE 
IDLI3T 
MAXSIZ 
LSCR 
RSCR 
EOF 


LOGICAL  Uty I T  NUMBER  OF  DATA  FILE  TO  BE  WRITTEN 

NUMBER  OF  DIMENSIONS  IN  EACH  DATA  POINT 

NUMBER  OF  POINTS  IN  A  LINE 

NUMBER  OF  LINES  IN  A  GRID 

GRID  NUMBER 

LATTICE  NUMBER 

STORAGE  SPACE  FOR  DIMENSION  LABELS  TO  BE  SAVED  - 
DATA  POINTS  ARRAY  (NDrNL*NL)  i 

STORAGE  SPACE  FOR  20  CHARACTER  SHORT  TITLE  S 

STORAGE  SPACE  FOR  60  CHARACTER  LONG  TITLE  :■ 

ARRAY  OF  INDICES  OF  DATA  TO  BE  SAVED  = 

DIMENSION  OF  SECOND  INDEX  OF  DAT  ARRAY  = 

L0GICAL*1  SCRATCH  ARRAY  (10*MAXSIZ)  c 

REAL  SCRATCH  ARRAY  (MAXSIZ)  ~ 

LOGICAL  END-OF-FILE  FLAG  -  TRUE  WHEN  EOF  WAS  RE Al£ 


LOG  I C AL  # 1  STITLEC20) ?LTITLE(60) >DLABEL(10»ND) »  LSCR ( 1 )  5 

LOGICAL  EOF 

DIMENSION  DAT  (ND»MAXSIZ)_fRSCR(  1  )  r  ID  LIST  (ND)  ,= 

11  FORMAT (515)  1 

READ  LINES  FROM  FILE  UNTJL  THE  GRID  NUMBER  CHANGES 
NF'  =  0 
NL  =  0 
NT0T  =  0 

10  CALL  GETLI N ( LUN » ND *IDLIST» MAXSIZ * LSCR. RSCR* ST ITLE> LTITLE* DLABEL > 

1  DAT ( 1 f NTOT+1 ) rNPl ,  NL1 t NGfNLATtEOF) 

NPsNP+NPI 

NL=NLrNLl 

NT0T=NT0T+NP1*NL1 

READ ( LUN » 11 >END=999)MNB»NP2>NL2rNG2rNLAT2 
BACKSPACE  LUN 
IF(NG2.EQ.NG)G0  TO  10 
RETURN 

999  CONTINUE 
RETURN 
END 

**************M*******M#********M************************************: 


ROUTINE  TO  GET  A  LINE  OF  DATA  FROM  A  FILE  IN  THE  STANDARD 
C  DFAN  DATA  FILE  FORMAT. 

C 

C 

C  MM***************************************************************** 

SUBROUTINE  GETL I N ( LUN , ND r I DL I ST , MAXS I Z > LSCR > RSCR , ST  I TLE f L T I TLE  . 
1  DLABEL  f  DAT , NP  »NL  »  NG  *  NLAT , EOF ) 

C  A.39 

r 


no  onnnnononoonnooooo 


1 1 
12 


50 

100 


999 


LUN 

ND 

IDLIST 

MAX3IZ 

LSCR 


R3CR 


3TITLE 

LTITLE 

DLABEL 

DAT 

NP 

NL 

NG 

NLAT 

EOF 


LOGICAL  UNIT  NUMBER  OF  FILE  TO  BE  READ  ( INPUT ) 
NUMBER  OF  VARIABLES  TO  GET  FROM  DATA  POINT ( INPUT ) 
ARRAY  OF  INDICES  OF  VARIABLES  TO  BE  RETURNED 
MAXIMUM  SIZE  OF  SECOND  DIMENSION  OF  DAT  ARRAY 
LOGICAL  SCRATCH  ARRAY  (10  X  NUMBER  OF  DIMENSIONS 
IN  EACH  DATA  POINT) 

REAL  SCRATCH  ARRAY  (NUMBER  OF  DIMENSIONS*  IN  EACH 
DATA  POINT  LONG) 

SHORT  TITLE  (20  CHARACTERS)  (OUTPUT) 

LONG  TITLE  (60  CHARACTERS)  (OUTPUT) 

DIMENSION  LABELS  ARRAY  (10*ND>  (OUTPUT) 

ARRAY  OF  DATA  RETURNED  (ND*MAXSIZ>  (OUTPUT) 

NUMBER  OF  POINTS  IN  THE  LINE(OUTF'UT) 

NUMBER  OF  LINES  IN  GRID  (OUTPUT) 

GRID  NUMBER  (OUTPUT) 

LATTICE  NUMBER  (OUTPUT) 

LOGICAL  FLAG  -  TRUE  IF  END  OF  FILE  ENCOUNTERED 


NUMBER  OF  POINTS 
AND  LATTICE  NUMBER 


LOG  I  CAL* 1  DLABEL ( 10 *ND> *  ST  I TLE (20) * LTITLE (60) *  LSCR 
LOGICAL  EOF 

DIMENSION  DAT ( ND *  MAX5 IZ ) *RSCR<1 ) * IDLIST (ND) 

FORMAT (515) 

FORMAT ( 10E12 . 5 ) 

EOF=. FALSE. 

GET  THE  NUMBER  OF  DIMENSIONS  PER  POINTr 
PER  LINE*  NUMBER  OF  LINES*  GRID  NUMBER* 

READ (LUN* 11 *END=999)MND*NP*NL*NG*NLAT 
GET  THE  HEADER  FOR  THE  LINE 

CALL  GETHDR(LUN.MND*NB* IDLIST  *LSCK* DLABEL *  ST  I TLE* LTITLE* EOF ) 

I F ( EOF ) RETURN 

COMPUTE  NUMBER  OF  POINTS  IN  LINE/GRID 
NL  =  MAX0 ( NL  *  1 ) 

NTOT=NL*NP 

GET  THE  DATA  POINTS  FOR  -THE  LINE/GRID 
DO  100  1  =  1*  NTOT 

READ (LUN* 12*EN0  =  999) (RSCR( J) * J  =  1 *  1ND) 

IF  THE  DAT  ARRAY  IS  FULL*  DO  NOT  SAVE  THIS  POINT 
IF ( I ,GT .MAXSIZ)GO  TO  100 
DO  50  J= 1  *  ND 

DAT ( J  *  I ) *RSCR ( IDLIST ( J ) ) 

CONTINUE 
CONTINUE 
RETURN 
CONTINUE 
EOF  = . TRUE . 

RETURN 

END 

************************************************** *  ******* **iMttt****** 


ROUTINE  TO  GET  THE  HEADER  RECORDS  FROM  A  FILE  IN  THE  STANDARD 
DFAN  DATA  FILE  FORMAT. 


*********************************************************************** 

SUBROUTINE  GETHDR ( LUN * MND * ND *  I DL I ST  * LSCR * DLABEL *  ST I TLE * L T I TL E  *  EOF 


C 

0 

c 

c 


LUN 

LOGICAL 

UNIT  NUMBER 

OF 

DATA  FILE  T0  BE  READ 

MND 

NUMBER 

OF 

DIMENSIONS 

IN 

A  DATA  POINT 

ND 

NUMBER 

OF 

DIMENSIONS 

TO 

SAVE  FROM  EACH  DATA  POINT 
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IDLIST 

LSCR 

DLABEL 

STITLE 

LTITLE 

EOF 


$RRAY  OF  INDICES  CF  VARIABLES  TO  BE  SAVED 
LQGICAL*1  SCRATCH  ARRAY  FOR  USE  IN  READ 
STORAGE  SPACE  FOR  DIMENSION  LABELS  T0  BE  SAVED 
STORAGE  SPACE  FOR  CO  CHARACTER  SHORT  TITLE 
STORAGE  SPACE  FOR  iO  CHARACTER  LONG  TITLE 
END-OF-FILE  rL AG  -  SET  TO  TRUE  IF  EOF  ENCOUr 


DIMENSION  I DLIST ( ND  ) 

LOG  I  CALK  I  STITLE  C  20  )  .  LTITLE  (  SO )  >  LSCR  (  10  .  MND  >  »  DLABEL  CIO*  ND 
LOGICAL  EOF 

11  FORMAT t 20A1 , cOAl  ) 

12  FORMAT ( 1QA1  ) 

READ  RECORD  2  -  GET  THE  TITLES 

READ ( LUN » 11 . END *999 ) STITLE . LTITLE 

READ  RECORD  3  -  GET  ALL  DIMENSION  LABELS 

READ(LUN.12.£ND=999)  ( CL5CRC I . J ) . I  *  1 .  1 0 )  . J- 1 . MND ) 

SAVE  THE  DESIRED  DIMENSION  LABELS  IN  DLABEL 
DO  100  J  =  1 .NO 
DO  50  1=1.10 

DLABEL  (I.J)=L3CR(I.IDLIST(J)  ) 

50  CONTINUE 

100  CONTINUE 

RETURN 

999  CONTINUE 

EOF= . TRUE . 

RETURN 

END 

******¥************************  *********  **  *********************  **  *** 


ROUTINE  TO  PUT  THE  HEADER  RECORDS  ON  A  FILE  IN  THE  STANDARD 
DFAN  DATA  FILE  FORMAT. 


t********??******;**?**^*****^******  ********  ************************ 

SUBROUTINE  P'JTHDR  (  LUN  .  ND  »  NR  .  NL  .  NG  .  NLAT  r  STITLE. LTITLE  -  DLABEL) 


C 


C 

C 


LUN  LOGICAL  UNIT  NUMBER  FOR  OUTPUT  DATA 

ND  NUMBER  OF  DIMENSIONS  IN  A  DATA  POINT 

NP  NUMBER  OF  POINTS  IN  A  LINE 


NL  NUMBER  OF  LINES  IN  A  GRID 

NG  GRID  NUMBER 

NLAT  LATTICE  NUMBER 

STITLE  20  CHARACTER  SHORT  TITLE 

LTITLE  SO  CHARACTER  LONG  TITLE 

DLABEL  10  CHARACTER  DIMENSION  LABELS 


LQGICALK1  ST  I TLE ( 20 ) . LTITLE ( SO ) . DLABEL  C  10.  ND 
FORMAT( 515) 

FORMAT ( 20A 1 . S0A1 ) 

FORMAT ( 1 0  A  1 ) 

OUTPUT  RECORD  *1 

WRITE  (  LUN  .  1  )  ND  •  NP  »  NL  .  NG  »  NLAT 

OUTPUT  RECORD  *2 

UR  I  TEC LUN. 12) STITLE. LTITLE 

OUTPUT  RECORD  *3 

UR  I TEC  LUN. 13) ( ( DLABEL ( I . J  )  .  1  =  1  .  10  >  .  J*1  .ND  ) 
RETURN 

END  A41 
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c 

c 

C  PLOTTING  PROGRAM  FOR  CONTOUR*  AXONOMETRIC.  AND  VECTOR 

C  FIELD  PLOTS  -  THIS  PROGRAM  ACCESSES  DATA  IN  DATA  FILE 

C  "DATAIN.DAT*  WHICH  IS  IN  THE  STANDARD  DF  AN  DATA  FILE 

C  FORMAT.  THE  FRCGRAM  IS  AN  EXTENSION  OF  '  F I E  L  D  F’  L  T  '  WHICH 

C  WAS  DEVELOPED  AND  WRITTEN  BY  CART  GLYNN  SISSON.  THE 

C  PROGRAM  HAG  BEEN  MODIFIED  AS  FOLLOWS: 

C 

C  MODIFICATION  RECORD: 

C 

C  INITIAL  MODS  AND  COMMENTS  29  JULY  31  CAF'T  BOLICK 

C  MORE  COMMENTS  ADDED  13  SEPT  Si  CAPT  BOLICK 

C 
C 

C  ***************  ********************* ******************************* 

c 

i: 

C  DATA  DECLARATION  AREA 

C 

C 

C******** *************  ********************************************** 

DIMENSION  Xl(10)rYl(10)»Zl(10)»Vl<10)fX2(10)fY2(10)» 

1  Z2 ( 10 )  *  V2 ( 10 ) 

DIMENSION  LINEC20) >  X  3 ( 5  0  )  » Y3 ( 50 ) . FARM ( 3 1 ) 

DIMENSION  F(  4  f  2)  »  DAT  (Sr  100  )  f  II'LIST  (5) 

LOG  I  CAL* 1  LSCR(10»1?) »3TITLE(20> .LTITLE(aO) »  NADE ( 1 0 ) »DLABEL(10»5 
LOGICAL  EOF 
L0GICAL*1  TITLE(AO) 

EQUIVALENCE  ( LSCR < 1 . 1 )  >X3  < 1 >  > 

EQUIVALENCE  (N1 . IDLISTd  ))  ,  :N2» I DLrSTCC)  ) t  (N3r IDLISTC3)  )  , 

1  (N4  » IDLIST ( 4  > ) » ( N5 . IDLIST(S)  ) 

C** *********************************************************** T****» 

C 

c 

C  DATA  PRESETS  AREA 

r 

C 

C**************»****  ************************************** 

DATA  YES/’Y'/ 

C  DISCRIPTOR  FOR  DIFFERENT  LINE  TYPES  TO  BE  DRAWN 

DATA  LINE/'77fl2*32f332f512fS312f72»732>34f334r514» 

*  5314.74. 734. 3212.5212. 7212. 3414, 5414. 7414/ 

DATA  LUN/2/ >  M  A  X  5  I  Z  / 100/ » IBAUD/960/ »  N  D  /  5  / 

d******f*HF***¥***********t******************************************* 

C 

c 

C  MAIN  PROGRAM  CODE  STARTS  HERE 

C 

C 

C**t* ********* t* ***#««¥*** V77********************************* 

C 

c 

C  GET  MENU  AND  INPUT  DATA  FILES  ASSIGNED 

C  SET  UP  TERMINAL  LOGICAL  UNITS  ALSO 

C 
C 

CALL  ASSIGN ( 3 ,  ' MENU  .  DAT  '  , 8 ) 

CALL  F D B S E T ( 3  ,  'READONLY'  ) 

C  READ  THE  INPUT  DATA  FROM  THE  MENU  DATA  FILE 

A43 
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o  noon  n 


REAEK3)  TITLE  *  OPT.  AF'.F.XMINf  XMAX  i X INC »  YM I N »  YMAX » Y I NC  > 
ZM  I N  »  ZMAX  r  Z I NC  » IDLIST 
SAVE  THE  SCALE  FACTOR  FOR  VECTOR  FLOTS 
SC ALEF  =  F ( 3  >  2 ) 

CLEAR  THE  FIELD  IF  THE  OPTION  IS  VECTOR  FIELD  PLOT 
IF  (OPT.GT.l. )F(3»2)=0. 

THE  OUTPUT  CRT  FOR  THE  GRAPH  WILL  BE  ,TT7*.* 

CALL  ASS  I GN ( 5  » ' TT7J ' ) 

IF  A  TK  4662  HARDCOPY  WAS  REQUESTED  CHANGE  UNIT  5  TO 
IF(AP .EQ. YES) CALL  ASS  I GN ( 5 r ' TT 1 J ' ) 

INITIALIZE  THE  GRAPHICS  PACKAGE  -  CLEAR  THE  SCREEN 
CALL  INITT ( IBAUD) 


INPUT 


FILE 


CALL  ASSIGN(2r 'DATAIN.DAT' r 10) 

CALL  FDBSET ( 2  > 'OLD' t  'SHARE' ) 

FILL  THE  F  ARRAY  BASED  ON  THE  MIN/MAX/ ANGLE  DATA 

CALL  FFILL(F»XMIN>XMAX»YMIN*YMAX*ZMIN>ZMAX>SXMIN>SXMAX' 

1  SYMIN  t SYMAX ) 

:  IF  THIS  IS  NOT  A  CONTOUR  PLOT  THEN  SKIP  TO  THE  NEXT  SECTION 

I F  <  OPT . NE . 0  * >  GO  TO  200 

SET  UP  THE  PLOTTING  MINMUM  AND  MAXIMUM  VALUES  AND  DEFINE 
:  THE  PLOTTING  WINDOW  FOR  THE  CONTOUR  PLOT. 

RANGE=SXMAX-SXMIN 
XM I =SXM I N- . 05*RANGE 
XMA=SXMAX+. 45* RANGE 
YMI=SYMIN-. 05* RANGE 
YMA=SYMAX+. 1 *R ANGE 
CALL  DWINDO(XMIfXMA»YMIfYMA) 

CALL  TWIND0<0> 1020»0*7GO> 

ZLOW=ZM I N 

************************************************************* 


CONTOUR  PLOTS  ARE  DONE  HERE 


Cl*ttt*******t*#**t*** ******************** ********************** 

100  CONTINUE 

C  GET  THE  NEXT  GRID/PLANE  OF  DATA 

CALL  GETGRDCLUNr ND» IDLIST . NAXS IZ » LSCR » X3 . ST  I TLE * LT I TLE > DLABEL » 
1  DAT  »NPrNLfNG»NLAT  »  E  0  F  ) 

C  IF  NO  MORE  DATA  THEN  PROCEED  TO  FINISH  UP 

IF  <  EOF ) GO  TO  130 
I CNT  =  0 
Z3=ZL0W 

110  CONTINUE 

ICNT=ICNT+1 
K 1  =  1 
K  2  =  N  P  + 1 
K  =  K  2 

DO  120  J  =  2  f  NL 

C  MOVE  A  LINE  OF  DATA  AT  A  TIME  INTO  THE  PLOTTING  ARRAYS 

CALL  FILLIT<  NPrKl »K2»Xl , Y1 , Z1 » Vl »  X2» Y2,Z2» V2»  DAT) 

C  GO  DO  THIS  CONTOUR  LINE 

CALL  CONTOR  (  NP  »  XI  »  Y 1  fZl  t  V 1  ?X2»Y2»Z2»V2»X3»Y3»Z3*LINE<  ICNT  :■  ' 

K1=K1+NP 
K  2  =  K  2  +  N  P 

120  CONTINUE 

C  SET  UP  FOR  NEXT  CONTOUR  LINE 

Z3=Z3+ZINC 

C  IF  NOT  ALL  CONTOURS  ARE  PLOTTED  THEN  PLOT  THE  NEX"  ONE 
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IF(Z3.LE.ZnAX»G0  TO  110 

C  AFTER  ALL  CONTOURS  ON  THIS  PLANE  ARE  PLOTTED  -  GO  SEE  IF 

C  THERE  IS  ANOTHER  PLANE  OF  DAT~  t 3  f'_0T 

GO  TO  100 
130  CONTINUE 

C  COMPUTE  LOCATIONS  -OR  THE  -EGZND  ON  the  GRAPH 

X-SXMAX+  .  1  *  ■  3XMAX  -SXn  I  M ' 

Y=SYMAX-  .  05*  (  S  (MAX  -  STM  IN  ' 

XI  -  . 25* ' SXMAX-SXMTN i 
ri=.04*(SYMAX-SYMIN) 

C  GO  PUT  LEGEND  ON  THE  GRAPH 

CALL  LEGEND  (  ICNT  *  LINE  »  ZLC-W  -ZINC*  '  •  '  •  <■  I  •  I  1 

C  COMPUTE  LOCATION  FOR  NAME  OP  UhRIAELE  PLOT’ED 

X=.5¥(SXMAX-SXMIN)tSXMIN 
Y=SYMAXr . 0  1  ¥ ( SXMAX-SXMIN ' 

CALL  MOVEA(X»Y) 

CALL  ANMODE 

C  OUTPUT  VARIABLE  NAME  THAT  HAS  -L2"E:' 

CALL  A  A  N  3  T  R  (  1  0  >  D  L  A  ?  t  L  '  1  »  _>  )  ) 

C  GO  TO  FINISH  UP  GRAPH 

GO  TO  SCO 

Clttt******tt********4ti«t**t ********************************* 

C 

c 

C  hXQ NOMETRIC  (3-D)  PLOTS  ARE  DONE  HERE 

C 

C 

C*** ************************** ********************************* 

200  CONTINUE 

C  IF  AN  AXONCMETRIC  PLOT  U AS  NOT  REQUESTED  T  H  E  N  GO  T 0  CROSS  FLOW 

C  PLOT  SECTION 

IF (OPT . N£ . 1 . ) GC  TO  300 

C  COMPUTE  RANGP  AMD  LIMITS  FOR  THE  AXONOMETRIC  PROJECTION  AND 

C  THEN  SET  THE  PROPER  LIMITS  IN  THE  WINDOW  DEFINITIONS  SO  THAT 

C  ALL  OF  THE  GRAPH  WILL  BE  ON  THE  SCREEN 

RANGE  =  AMAX1 ( 3XMAX  f SYMAX) -AMIN1 CSXMINf  SYMIN) 

XM I =SXM IN- . 05*RANGE 
XMA  =  SXMIN+1  , 2TRANGE 
YMI=SYMIN- . OOP RANGE 
YMA  =  SYMIN+1  .  1* RANGE 
CALL  DUINDO ( XMI » XMA •  VMI  . YM A ) 

CALL  TWIND0(125r975r 0,730) 

210  CONTINUE 

C  GET  THE  NEXT  GRID/PLANE  OF  DATA 

CALL  GETGRD(LUNfND»IDLI3TfMAXS:Z,LSCR»X3,3TITLE>LTITLE-DLABEL> 
1  DAT t NP < NL t NG t NLAT » EOF ) 

C  IF  THERE  IS  NO  MORE  DATA  THEN  GO  TO  FINISH  UP  THIS  FLOT 

IF (EOF) GO  TP  200 

C  COMPUTE  TOTAL  NUMBER  OF  DATA  FOINTS  IN  THE  GRID 

NT OT-NP*NL 
K  1  =  1 
K2=NP + 1 
DO  220  J  =  2  >  N  L 

C  MOVE  THE  DATA  INTO  THE  PLOTTING  ARRAYS  ONE  LINE  AT  A  'r  I  ME 

CALL  FILLIT(NP  f  K1  r  R2  »  XI  ,  Y1 . 21  »  VI  f  X2  r  Y2  ►  Z2  »  02  '  DAT  '■ 

C  IF  THIS  IS  THE  FIRST  RAIL  FOR  THIS  PLANE  THEN  CALL  D3’«E3H 

C  TO  PROPERLY  SET  UP  TO  DRAW  THE  MESH 

IF  (  J  .  EQ  .  2  )  CALL  D3MESH  (  >  NP  7  X  1  >  Y  1 , 1  \  r  U  I  ,  X  1  »  Y  1  .  Z  1  •  V  1  '■ 

C  DRAW  THE  NEXT  SET  OF  RAILS  FOR  THIS  PLANE 

CALL  D  3  M  E  S  H  (  F  *  N  P  »  X  1  >  Y  1  7  Z  1  •  V  1  *  X <  Y  2  f  Z  2  f  0  2  ) 

K  1  =  N  1  P  N  P 

A.  >  5 


frfc *  C  l  u  tin  iff  « 
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K2=K2+NP 

220  CONTINUE 

C  GO  SEE  IF  THERE  IS  ANY  MORE  DATA  TO  PE  PLOTTED 

GO  TO  210 
250  CONTINUE 

C  COMPUTE  LOCATION  TO  OUTPUT  NAME  OF  VARIABLE  THAT  WAS  PLOTTED 

X  = . 5* ( 3XMAX-SXMIN ) +SXMIN 
Y=SYMAX+.01*(SXMAX-SXMIN) 

CALL  M 0 V E A ( X  r  Y  > 

CALL  AN MODE 

C  OUTPUT  VARIABLE  NAME 

CALL  AANSTR( 10>DLABEL( 1 »3)  ) 

GO  TO  500 

C#**#* ******************************************************** 

c 


VECTOR  FIELD  PLOTS  ARE  DONE  HERE 

1 


t***4ttt444*4**t4****t*tt****4t**t*4**t*44ttt*444*tt*t**4t44M 

300  CONTINUE 

COMPUTE  RANGE  AND  LIMITS  FOR  PLOTTING  WINDOWS  FOR  THE  CROSS 
FLOW  VECTOR  PLOTS. 

RANGE=3XMAX-SXMIN 
XMI=SXMIN- . 05*RANGE 
XMA=SXMAX+ . 2*RANGE 
YMI=SYMIN- . 05* RANGE 
YMA  =  SYMAX+ . 1  GRANGE 
CALL  DWINDOCXMIf  XMAfY-MH  YMA> 

CALL  TUIND0(125»975»0»780) 

310  CONTINUE 

GET  THE  NEXT  GRID/PLANE  OF  DATA 

CALL  GETGRD ( LUN  >  ND  r ID LIST  »MAXSIZfL3CRfX3*5TITLEfLTITLE*D LABEL » 
1  DAT » NP » NL ? NG » NLAT , EOF ) ^ 

IF  NO  MORE  DATA  THEN  GO^TO  FINISh  UP  THE  GRAPH 
IF ( EOF ) GO  TO  500 
NTOT  =  NP#NL 

CALL  THE  SUBROUTINE  TO  PLOT  EACH  OF  THE  VECTORS 
CALL  VFIELDKNTOTf DATr SCALEF) 

GO  SEE  IF  THERE  IS  MORE  DATA  TO  BE  PLOTTED 
GO  TO  310 
00  CONTINUE 

OUTPUT  GRID  AROUND  THE  PLOT  SPACE 

CALL  GRID ( F. •  XMINrXMAXfXINCf  YMIN , YMAX. Y INC  ) 

CALL  MOVABS < 100 f 760 ) 

TITLE ( 46 ) =32 
TI TLE ( 47)=32 

C  OUTPUT  THE  TITLE  OF  THE  GRAPH 

CALL  AANSTR < 60 f TITLE > 

C  CLOSE  OUT  THE  PLOTTING  PACKAGE 

CALL  FINITT(OfO) 

CALL  EXIT 
END 


A4N 
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APPENDIX  B 


This  appendix  contains  samples  of  the  user  interactions  required  for 
both  the  data  acquisition  and  plotting  programs. 


.touio.t  »  c  )  u  no  ft  mi 


1  _ _ 

■  1  ■:  C  C  IJ  ;-l  T  0  R  N  A  !■:  E  :  fc~  I  i-  r  I  '■< 


=  A  3  3  W  G  F.  L 


R3X-UM  BL2  =  MUL”  I  -USE?:  SYSTEM 


•CD  MORNING 


*  -  n  r  r  -  a 


.'Juut:1  _i  I 


WELCOME  TO  R5X-11M  V3.2  REV  D  fqr  OF  AN 


lE-ClC?-  L^OCi.PSh'Pl 

;  K  *  x  X  M  *  *  t;  *  *  *  ;(c  :*  :f.  A  *  >:  A  A  A  *  A  t  *  A  A  X  fc  A  A  A  *  *  A  A  A  A  A  *  >  A  A  A  A  A  A  A  M  »  I  ‘  A  *  A  *  A  *  *  M  A 


COMPRESSIBLE  SEVEN  HOLE  PROBE  FLOW  FIELD  SUP 
DATA  ACQUISITION  SYSTEM  COMMAND  FILE 


;  K  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  A  n  A  *  A  A  A  A  A  A  A  A  A  A  A  %  %  %  AAA 
PIP  TARE. DAT? A/DE/NM 


PIP  MENU. DAT? A/DE/NM 
PIP  D A  t  p I N . D  A  T  r A / D  E / N H 
:  PIP  DATAOUT . DAT ! A/DE/NM 

t  DO  YOU  WANT  TO  ADD  DATA  TO  AN  EXISTING  FILE"  CY/N3 :£) 


:  PIP  DATAIN . DAT/NV-XWM4 . DAT 
'  PIP  TARE. DAT/NV=C 152>13jNULL.DAT 
ARUM  l1S2»  133CMPSHP 
THE  PRESENT  DATA  TITLE  I S  t 
"WING-BODY,  il  DEG  t  V-lOO  EPS  < XWM4 ) 
DO  YOU  WISH  TO  CHANGE  IT? 


ENTER  DATA  TITLE  (AO  CHARS  MAX) 

777%  PF  -  I  p  H  FIT;  0  s;jp'UEY  Fr-f-PANl 
START  TUNNEL  AND  ENTER  RETURN  WHEN  READ v 


INTER 


M  OSF'HERIC  F’  R  E3  3  U  R  E  (  I  N  -  H  G  ) 


COORDINATE  PAIRS  ARE  (VERT.HORIZ) 

ENTER  LOCATION  OF  LOWER  LEFT  AND  ENTER  RETURN 


ENTER  LOCATION  OF  UPPER  LEFT  AND  ENTER  RETURN 

£=3 

ENTER  LOCATION  OF  LOWER  RIGHT  AND  ENTER  RETURN 


ENTER  LOCATION  OF  UPPER  RIGHT  AND  ENTER  RETURN 

cm 

LOWER  LEFT;  -3.000  6  .  0 0  0 

UPPER  left:  6.000  6.000 

LOWER  RIGHT:  -3.000  3.000 

UPPER  right:  7.000  S.000 

ARE  THESE  LOCATIONS  CORRECT" 


DATA  POINT:  VERTICAL^ 


HOP  ICO  N  T  A L  = 


-IP  DATAIN . D  A  T / N  V  - D  A  T  A  0  U  T . DAT 
■  P I p  DATAOUT  .  DAT . * / D  E . '  M M 
'  PIP  DATAIN  .  DAT/PIJ 

.  M  DO  YOU  WANT  rQ  USE  AN  EXISTING  MENU  DATA  FILE7  L ' 


& 


niiHij u >  Kt'\  w  ii'»  i«.nim>o.rM.‘w 


CO  CO  CO 


AD-A115  971  AIR  FORCE  ACADEMY  CO 

SEVEN-HOLE  PROBE  DATA  ACQUISITION  SYSTEM. (U) 
NOV  81  A  GERNER  *  6  SISSON 
UNCLASSIFIED  USAF A-TN-81-8 _ 

2  2 


13-3EP-91  08  .‘43 


FROM  TT 10 


TO  TTO 


■IE  SRlFFIN/OEWq 

PSX-11M  BL26  MULTI-USER  SYSTEM 
GOOD  MORNING' 

13-SEP-31  08J30  LOGGED  ON  TERMINAL’  TTO t 
WELCOME  TO  RSX-11M  V3.2  REV  D  FOR  DFAN 


T*:M:M****************************:|U***r4********************** 


FLOW  FIELD  SURVEY  PLOTING  COMMAND  FILE 


«7:ln 


s'  ?  * 

> }  .It************************************************************* 

>•*  ENTER  NAME  OF  DATA  FILE  TO  BE  PLOTTED  CS3.* 

>PIP  DAT  A  IN  .  D"AT/NV«XWM4  •  DAT 

>t  DO  YOU  WANT  TO  USE  AN  EXISTING  MENU  DATA  FILE?  CY/N3  4I 
r  P IP  MENU . DAT/NV»C132» 13 3 NULL . DAT 
7RGN  C132» 133MENU 


ENTER  PLOT  OPTION:  0«C0NT0UR»  l*AXONOMETRICr  2»CR0SS  FLOW 

(S 


DATA  XMIN  -  1.000 


DATA 

XMAX  * 

6.000 

DATA 

YMIN  » 

-3.000 

DATA 

YMAX  » 

7.000 

DATA 

ZMIN  a 

-1.913 

DATA 

ZMAX  a 

1.40 

I.  OPTION  *  0. 

2.  TITLE  ■  FWING-BODY t  11  DEGf  V-100  FPS  (XWM4)  13 : 06 : 324-AUG-81 

3. X  AXIS  ANGLE  »  0.000 

4. Y  AXIS  ANGLE  «  90.000 

3 . Z  AXIS  ANGLE  >  0.000 

6 • X  MIN  >  0.000 

7.X  MAX  ■  10.000 

3 . Y  MIN  »  -4.000 

9. Y  MAX  *  6.000 

10.  2  VARIABLE  «  10.  CTOTAL 

II. ZMIN  ■  -1.50000 

12. ZMAX  •  1.40430 

13.  ZINC  «  0.30000 

ENTER^NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 

•  •  _ ...  . 


B4 
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r 


•OPTION  *  0. 

•  TITLE  «  FUING-B0DY  t  11  DEO  »  V*100  FPS  (XUM4)  13 t 06 : 324-AUG-81 

•  X  AXIS  ANGLE  *  0.000* 

.Y  AXIS  ANGLE  »  90.000 

.2  AXIS  ANGLE  *  0.000 

6 • X  NIN  *  -2.000 

7. X  MAX  »  8.000 

8. Y  MIN  «  -2.000 

9 .  Y  MAX  a  8.000 

10.  2  VARIABLE  *  10. 

11. ZMIN  »  -l.SOOOO 

12 .  ZMflX  *  1.40450 

13.  ZINC  «  0.20000 

ENTER  NUMBER  OF  PARAnETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 

S 

DO  YOU  WANT  A  HARDCOPY  PLOT7  EJ 
7  8#d  8 

TTO  —  STOP 
>PIP  MENU.DAT/PU 
>RUN  C1S2» 133GRAFS 


FWlKC-BtJdY,  11  DEC.  U-183  FPS  <KHM>  15: M  24-4UG-81 


-1.S88 
- 1.388 
-1.188 
-4.988 
-8.788 
-8.588 
-8.388 
-8. 188 
8.188 

8.388 
8.588 
8.788 
8.988 
1.188 

1.388 


TTO  —  STOP 

>*  DO  YOU  WANT  ANOTHER  GRAPH?  CY/NIlgJ 

>*  DO  Yoy  WANT  TO  USE  A  DIFFERENT  MENU  FILE?  CY/N3.*® 

7RBN  ClI2> 133MENU 


tO  YOU  WANT  TO  SEE  THE  DATA  MINIMUM/MAXIMUM  VALUES?  CY/N3 

@  85 


Bit 


ui’lJiouuceu  on  a  os  (.uvernincn'  ->pier, 
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.OPTION  *  0. 

2.  TITLE  ■  FWING-BODY.  11  DEG.  V»100  FPS  ( XWM4 ) 

3. X  AXIS  ANGLE  «  0.000 

4 .  Y  AXIS  ANGLE  «  90.000 

5.2  AXIS  ANGLE  «  0.000 

6. X  MIN  «  -2.000 

7. X  MAX  ■  3.000 


15:06:524-AUG-81 


8. Y  MIN  *  -2.000 

9 .  Y  MAX  «  8.000 

10.  2  VARIABLE  «  10. 

11.2MIN  «  -1.50000 

12.2MAX  •  1.40450 

13.2INC  «  0.20000 

ENTER  NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 


.OPTION  *  1. 

2.  TITLE  «  FWING-BODY.  11  DEG.  V»100  FPS  ( XWM4  )  15  J 06 *. 524-AUG-81 

3 .  X  AXIS  ANGLE  *  0.000 

4 .  Y  AXIS  ANGLE  -  120.000 

5. Z  AXIS  ANGLE  •  90.000 

6. X  MIN  «  -2.000 

7. X  MAX  «  8.000 

8. Y  MIN  •  -2.000 

9 .  Y  MAX  *  8.000 

10.  2  VARIABLE  ■  10.  CTOTAL  - 

14.2  SCALE  FACTOR  -  -1.000 

ENTER  NUMBER  OF  PARAMETER  TO  GRANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 


•tin-; 


.OPTION  ■  1. 


2. TITLE  «  FWING 

-BODY  f  11  DEG 

3 .  X 

AXIS  ANGLE 

*  0.000 

4  ,  Y 

AXIS  ANGLE 

«  120.000 

5.2 

AXIS  ANGLE 

a  90.000 

6.X 

MIN  « 

2.000 

7.X 

MAX  » 

8.000 

8.Y 

MIN  » 

2.000 

9 .  Y 

MAX  « 

8.000 

10. 

2  VARIABLE 

a  10.  CTOTAL 

V«100  FPS  ( XWM4.) 


15  5  06  ?  52 4 -AUG- 81 


14. Z  SCALE  FACTOR  - 


-5.000 


ENTER  NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 


DO  YOU  WANT  A  HARDCOPY  PLOT?  fl 

7  t#d  8 

TTO  —  STOP 
>PIP  MENU.DAT/PU 
>RUN  C152. 133GRAFS 


B6 
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15:»<  24-ftUG-ei 


fuimo-sqoy,  u  ofc,  u-nj«  fps  <xnn4> 
CTOTAL 


TTO  --  STOP 

>*  DO  YOU  WANT  ANOTHER  GRAPH?  CY/NI*.Q 


>*  DO  YOU  WANT  TO  USE  A  DIFFERENT  MENU  FILE?  CY/NItfit 
7RBN  C 152  > i 3  3  MENU 

e~ 

DO  YOU  WANT  TO  SEE  THE  DATA  MINIMUM/MAXIMUM  VALUES?  CY/NI 

IS 


.OPTION  *  1. 

.TITLE  *  FUlNG-BODYf  11  DEG  »  V»lOO  FPS  ( X  WM4 )  13  :  06 : 324-AUG-81 

.X  AXIS  ANGLE  *  SO. 000 

.Y  AXIS  ANGLE  ■  90.000 

.2  AXIS  ANGLE  ■  90.000 

6. X  MIN  *  '-2.000 

7. X  MAX  »  8. 000 

8 .  Y  MIN  »  -2.000 

9. Y  MAX  *  8.000 

10.  2  VARIABLE  »  10. 

14.2  SCALE  FACTOR  -  -5.000 

ENTER  NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 

Eg 
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.OPTION  «  2. 

2.  TITLE  »  FWING-BODY »  11  DEG>  V-lOO  FPS  ( XUM4 )  15t06*.524-AUG-81 

3. X  AXIS  AN8LE  ■  0.000 

4 .  Y  AXIS  ANGLE  ■  90.000 

5 .  Z  AXIS  ANGLE  •  0.000 

6. X  MIN  ■  -2.000 

7. X  MAX  «  S.000 

8 .  Y  MIN  «  -2.000 

9 .  Y  MAX  ■  8.000 

10.  VECTOR  LENGTH  SCALE  FACTOR  ■  0.500 

ENTER  NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 


•i-m 


.OPTION  «  2. 

2.  TITLE  *  FWING-BODY t  11  DEG*  V-100  FPS  ( XUM4 )  13.*06:524-AUG-8l 

3. X  AXIS  ANGLE  >  0.000 

4 .  Y  AXIS  ANGLE  *  90.000 

5. Z  AXIS  ANGLE  «  0.000 

6. X  MIN  *  -2.000 

7. X  MAX  «  8.000 

8 .  Y  MIN  *  -2.000 

9 .  Y  MAX  ■  8.000 

10.  VECTOR  LENGTH  SCALE  FACTOR  ■  1.000 

ENTER  NUMBER  OF  PARAMETER  TO  CHANGE  AND  NEW  VALUE  (IF  APPROPRIATE) 

M  YOU  WANT  A  HARDCOPY  PLOT?  Q 

7  BOd  8 

TTO  —  STOP  *■ 

>PIP  MENU.DAT/PU 
>RUN  C 152 > 133GRAFS 

FN1N6-M0Y.  II  MG.  U>1M  PM  (XUM4>  131 M  24-AII6-61 


(.969 


2.M9 


•.MB 


r 


TTO  —  STOP 

>*  DO  YOU  WANT  ANOTHER  GRAPH?  CY/N3:fi) 

>*  DO  YOU  W£NT  TO  SAME  THE  MENU  DATA  FILE?  CY/N3 l Q 
>PIP  MENU «DAT»*/DE/NM 
>PIP  DATAIN.DATJ*/DE/NM 

>*  DO  YOU  WANT  TO  PLOT  ANOTHER  FILE?  CY/N*3:Q 
>e  <EOF> 
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APPENDIX  C 

DEAN  Standard  Data  File  Format 

File  structure:  Formatted,  sequential 
Record  1:  Format  515 

5  Integer  entries  as  follows: 

1  -  ND  -  Number  of  dimensions  in  a  data  point 

2  -  NP  -  Number  of  data  points  per  line 

3  -  NL  -  Number  of  lines  per  grid 

4  -  NG  -  Grid  number 

5  -  NLAT  -  Lattice  number 

Record  2:  Format  20A1 ,  60A1 

2  Logical*l  Arrays  as  follows: 

1  -  STITLE  -  20  character  short  title 

2  -  LTITLE  -  60  character  long  title 

Characters  41-48  are  Time-HH: mm: ss 
Characters  49-60  are  Date-DD-mmm-yy 

Record  3:  Format  10A1 

ND  Logical*l  Dimension  Names/Labels 
1  -  ND  -  DLABEL  -  10  Character  Dimension  Name/Label 


Record  4:  Format  NP*NL+3:  Formac  10E12.5 

ND  Real  numbers  for  a  data  point 

1  -  ND  -  DAT  -  Values  for  each  dimension  in  a  data  point 


This  basic  structure  may  be  repeated  as  many  times  as  necessary  in  a  single 
file. 
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